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ABSTRACT  OF  THESIS 


IDENTIFICATION  OF  POTENTIAL  AIRCRAFT  ICING  REGIONS 
THROUGH  MULHSPECTRAL  ANALYSIS  OF  GOES-8  IMAGERY 

Aviation  weaflier  forecasters  are  still  greatly  challenged  to  forecast  aviation  weather  ha2ar(ls  for 
large  areas  of  atmosphere  with  few  in-situ  observations  and  little  detailed  information  from  which  to  woric. 
Forecasters  generally  lack  sufficient  guidance  to  assist  them  in  the  prediction  of  icing  occurrences  or 
intensity.  In-flight  aircraft  icing  remains  a  significant  aviadon  hazard.  Icing  intensity  is  related  to  the  rate 
of  accretion  as  well  as  die  aircraft  type  and  the  icing’s  effects  on  aircraft  flight  performance.  Icing  regions 
currently  can  only  be  identified  from  subjective  reports  made  by  aircraft  pilots  (PIREPs). 

The  severity  of  aircraft  icing  is  found  to  be  quite  sensitive  to  temperature,  liquid  water  content, 
and  droplet  size  distribution  along  the  flight  path  within  the  icing  cloud.  The  difficulty  of  their  direct 
measurement  and  the  variability  of  these  factors  with  altitude,  position,  and  time,  coupled  with  variable 
aircraft  sensitivities,  make  forecasting  and  identifying  icing  environments  quite  difficult 

Widi  the  advent  of  the  next  generation  of  Geostationary  Operational  Environmental  Satellites 
(GOES),  additional  spectral  channels  with  improved  ground  resolution  are  available  for  use  in  a  wide 
variety  of  research  iqiplications.  One  important  outcome  of  applied  research  based  on  improved  GOES-8 
imager  products  is  analysis  technique  developments  that  will  lead  to  better  forecasts  of  hazardous  aircraft 
flying  conditions. 

In  tiiis  study,  GOES-8  multispectral  satellite  images  and  imagery  data  are  manipulated,  combined, 
and  enhanced  to  retrieve  clues  about  cloud  phase  (ice  or  liquid),  supercooled  cloud  liquid  water 
(SLW),  and  cloud-top  temperature  (CTT).  To  retrieve  information  on  cloud  phase,  SLW,  and  CTT,  tire 
known  spectral  characteristics  of  several  GOES-8  channels  (specifically  1, 2,  and  4)  are  used  to  interrogate 
the  cloud  top’s  radiative  properties.  This  information  is  used  to  determine  whether  aircraft  icing  is  likely 
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widlin  tiiese  dirougb  plication  of  a  simple,  multispectral  analysis  scheme.  The  penod  of  study 
consists  of  known,  significant  in-flight  aircraft  icing  events  ftiat  occurred  in  April  1995. 

TTiis  thesis  includes  a  survey  of  die  literature  on  previous  and  current  icing  studies  and  the  use  of 
PIREPs  as  a  verification  tool,  research  on  new  techniques  using  improved  satellite  imagery  to  better 
identify  regions  of  aircraft  icing  potential,  and  suggestions  for  future  work  or  research  towards  the  goal  of 
the  improvement  of  aircraft  icing  forecasts. 

The  multispectral  channel  difference  method  is  most  useful  when  only  large-scale  stratiform  cloud 
regions  are  considered.  In  these  two  case  studies,  scatter  plots  of  the  GOES-8  channel  2  cloud  reflectance 
component  versus  channel  1  cloud  albedos  show  how  die  more  intense  icing  PIREPs  cluster  within  a  small 
range  according  to  the  clouds'  microphysical  characteristics.  Inclusion  of  convective  clouds  to  this 
multispectral  scheme  leads  to  erroneous  assumptions  of  cloud  phase  and  SLW  based  on  their  colder  CTTs 
when  this  technique  is  applied. 

In  both  cases,  die  majority  of  die  icing  PIREPs  occurred  widiin  areas  with  CTTs  colder  than  die 
typical  icing  CTT  zone  (IZ)  (73%  in  Case  1),  however  many  significant  intensity  icing  PIREPs  did  occur 
within  the  IZ  in  Case  2  (44%)  when  the  cloud  shield  became  more  stratiform.  The  icing  PIREP 
verification  was  fairly  inconclusive  for  stratiform  regions  due  to  a  relative  lack  of  icing  reports  in  some 
portions  of  the  stratiform  cloud  deck. 

The  multispectral  RGB  color  oihancement  method  does  a  good  job  at  distinguishing  clouds  of 
different  heights  with  die  greatest  color  and  shade  contrast  between  the  high  and  low  clouds.  Some 
information  on  cloud  phase  and  SLW  clouds  seems  evident  from  inspection  of  the  middle  cloud  deck  color 
hues,  but  the  result  is  not  yet  reliably  quantitative.  This  information  could  be  used  qualitatively  to  help  in 

initial  icing  assessment  for  larger-scale  cloud  regions. 

Detection  of  potential  icing  environments  using  infiared  and  visible  satellite  imagery  remains  a 
challenge  for  cloud  regions  made  up  of  convective  clouds  or  a  mix  of  both  convective  and  stratiform  types. 

Bradford  D.  Schrumpf 
Department  of  Atmospheric  Science 
Colorado  State  University 
Fort  Collins,  CO  80523 
Sununer 1996 
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Chapter  1  —  Introdoction  /  Motivation 
1.1  Overview  of  the  Aircraft  Icing  Hazard 

In-flight  aircraft  icing  poses  a  significant  hazard  to  die  aviation  conununity,  as  witnessed  by  the 
crash  of  a  airplane  in  northern  Indiana  in  die  fidl  of  1994.  Ihis  crash  had  a  major  impact  on  die 

commuter  airline  safety  regulations  in  the  U.S.,  which  resulted  in  significant  economic  and  other  impacts 
on  die  airline  industty  as  a  whole  (Brown,  1995b). 

The  great  hazards  of  icing,  including  structural  airfiame  icing  and  induction  (carburetor)  icing,  on 
botii  fixed-wing  and  rotary  aircraft  have  long  been  realized.  Widi  over  1000  aviation  accidents  in  die 
continental  United  Stales  during  the  last  twenty  years  (Cole  and  Sand,  1991)  citing  aircraft  icing  as  a  major 
factor,  it  is  a  problem  that  has  not  been  overlooked.  However,  due  to  the  relative  inaccuracy  of  icing 
forecasts  during  thk  period  and  lack  of  detailed  information  fiom  which  to  work,  aviation  weather 
forecasters  have  not  been  very  successful  at  pinpointing  altitudes  and  areas  of  icing  hazard  without 
creating  large  felse-alarm  (over-forecasted)  areas.  When  this  over-forecasting  occurs,  the  pilots’ 
confidence  in  these  forecasts  go  down;  analogous  to  the  litfle  boy  who  cried  wolf  too  often.  Consequently, 
pilots  are  more  apt  to  ignore  these  forecasts  and  sometimes  cause  accidents  as  a  result. 

Cole  and  Sand  (1991)  showed  compelling  statistics  to  support  diis  view  and  their  findings  are 
summarized  in  Chapter  2.  They  found  tiiat  it  is  not  predominantly  die  low-time,  inexperienced  pilots  who 
are  involved  in  icing  accidents.  Apparently,  many  pflots  are  not  taking  appropriate  evasive  or 
precautionary  actions,  when  hazardous  icing  conditions  are  encoimteredtiuring  flight  or  forecasted  in  their 
required  weatiier  briefings  prior  to  takeoff.  Surprisingly,  nearly  60  percent  of  those  pilots  involved  in 
aircraft  icing  accidents  had  been  properly  informed  on  the  icing  tiueat  and  conditions  during  pre-flight 
weather  briefings.  This  suggests  that  there  exists  a  lack  of  understanding  or  concern  in  the  seriousness  of 
icing's  effects  on  the  performance  of  the  aircraft,  a  lack  of  faith  in  tiie  forecast,  and/or  a  general  lack  of 
understanding  of  the  many  factors  contributing  to  icing  conditions. 

Ultimately,  it  is  the  pilof  s  duty  to  know  tiie  strengths  and  weaknesses  of  their  type  of  aircraft  and 
to  evaluate  its  sensitivities  to  icing  conditions  for  different  environments  before  making  the  flight  through  a 
potentially  hazardous  region.  On  the  other  hand,  it  is  the  aviation  forecaster  that  must  be  considerate  of  the 
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flying  community’s  needs  and  by  striving  to  maximize  die  icing  forecast’s  precision,  so  as  to  limit  the 
spatial  extent  (bofli  horizontally  and  vertically)  of  areas  which  would  prohibit  reasonable  access  to  a 
substantial  variety  and  number  of  aircraft.  Thus,  for  aviation  purposes,  die  goal  of  any  icing  forecast 
system  and  for  most  recent  aircraft  icing  studies  should  be  to  minimize  the  forecasted  area  while 
tnaintaming  a  Wgh  dctection/verification  rate.  It  is  undesirable  to  over-forecast  regions  of  icing  and 
unnecessarily  make  them  virtually  inaccessible  to  aircraft.  This  goal  is  also  true  for  AIRMETs,  the  current 
operational  aviation  hazard  forecast,  as  well  as  for  the  forecast  model  algorithms  and  remote  sensing 
techniques  that  are  being  developed  to  improve  them  (Brown,  1996b). 

1.2  Types  and  Causes  of  Icing 

An  abundance  of  literature  has  identified  and  verified  die  general  meteorological  conditions 
required  for  aircraft  in-flight  icing  in  clouds  (Schultz  and  Politovich,  1992;  USAF,  1992;  Hansman,  1989; 
AWS,  1980;  Politovich,  1993;  etc.)  widi  some  new  fiu:tors  presently  being  discovered  and  further 
researched.  Many  of  these  factors  are  discussed  fiirdier  in  Chapter  2,  but  are  summarized  below. 

Two  distinct  processes  control  the  accretion  of  ice  onto  aircraft  surfaces  as  determined  by 
Hansman  (1985).  First,  supercooled  liquid  water  (SLW;  water  existing  colder  than  O^C)  cloud  droplets 
must  impact  the  air  frame.  Secondly,  these  SLW  droplets  must  freeze  onto  the  airframe  as  allowed  by  local 
heat-transfer  processes  along  the  airframe.  As  supercooled  liquid  cloud  droplets  impact  on  the  aircraft, 
they  can  splash  and  freeze  immediately,  or  drey  can  flow  back  along  die  airframe  somewhat  before 
freezing.  If  the  droplet  can  release  its  latent  heat  relatively  quickly  (typically  in  an  icing  environment  of 
colder  temperatures,  smaller  droplets,  and  lower  total  liquid  water  mass),  it  will  freeze  at  the  impingement 
point,  and  rime  ice  wiU  fonn.  This  is  usually  most  dangerous  when  occurring  in  large  amounts  on  smaller 
aircraft.  When  the  impacting  droplets  have  more  time  to  flow  prior  to  freezing  due  to  a  slower  removal  of 
latent  heat  (^ically  in  an  icing  environment  of  warmer  (sub-fieczing)  temperatures,  larger  droplets,  and 
higher  total  water  mass),  clear  ice  or  glaze  wUl  form.  This  type  of  icing  is  much  harder  to  detect  visually 
by  the  pilot,  which  gives  less  time  to  take  precautionary  measures  against  it.  While  clear  icing  is  the  most 
uncommon,  it  is  the  most  dangerous  type  for  aircraft.  Transitions  between  clear  and  rime  icing  can  occur 
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wMun  tile  same  cloud,  and  fligiit  tiirough  tiiese  conditions  can  cause  botii  types  of  ice  to  accumulate.  This 
condition  is  often  reported  as  mixed  ice. 

In  the  past,  tiie  method  of  forecasting  meteorological  conditions  favorable  to  aircraft  icing  was 
largely  a  of  predicting  cloudy  environments  (as  tiie  source  for  tiie  generation  of  SLW)  in  tiie 

appropriate  temperature  range  (between  +4  and  -40  degrees  Celsius).  As  described  in  further  detail  in 
Chapter  2,  the  type  and  intensity  of  tiie  ice  accretion  is  strongly  and  non-linearly  dependent  on  many 
environmental  variables  such  as  temperature,  liquid  water  content  (LWC),  cloud  drpplet  sizes  and 
distributions,  cloud  water  phase,  and  turbulence  or  shear  layers  aloft.  Additionally,  the  intensity  of  tiie 
icing  and  its  effect  on  aircraft  performance  depends  on  both  the  type  and  condition  of  tiie  aircraft.  The 
icing  forecasting  challenge  lies  in  the  fact  that  most  of  these  factors  are  somewhat  hard  to  determine  due  to 
the  natural  variability  of  clouds  (Schultz  and  Politovich,  1992;  Pobanz,  1994),  and  this  characteristic  makes 
it  difficult,  if  not  impossible,  to  completely  isolate  the  Individual  effects  of  tiiese  icing  factors  in  nature. 
One  would  like  to  ultimately  quantify  tiiese  variables  prior  to  flight,  yet  scientific  methods  at  present  do  not 
allow  for  this. 

13  Current  Aircraft  Icing  Forecast  Methods 

Presently,  there  exists  no  formal  training  program  for  icing  forecasting  within  the  NWS. 
Prediction  is  based  on  a  loose  set  of  guidance  took  and  tiieir  emphasis  can  vary  with  each  forecaster.  This 
results  in  a  variety  of  methods  in  use  by  icing  forecasters  at  tiie  Aviation  Weaflier  Center  (AWC,  formerly 
the  National  Aviation  Weather  Advisory  Unit  (NAWAU)),  which  in  turn  leads  to  high  degree  of  variability 
in  forecast  accuracy,  since  the  method  used  for  each  forecast  k  not  recorded  (Politovich  and  Olson,  1991). 
The  standard  aircraft  icing  forecasts,  or  AIRMETs  (Airmen's  Meteorological  Information),  are  created  by 
NWS  forecasters  using  a  wide  variety  of  available  information  and  are  issued  operationally  every  sue  hours 
and  are  amended  as  necessary  by  the  National  Oceanic  and  Atmospheric  Administration’s  National 
Weather  Service  (NOAA/NWS)  Aviation  Weather  Center  (AWC)  in  Kansas  City. 

The  Air  Force  ako  issues  time-phased  icing  forecasts  for  the  Northern  Hemisphere  for  the  layer 
between  10,000  and  55,000  feet  from  tiieir  Air  Force  Global  Weaflier  Center  (AFGWC)  at  Offutt  AFB, 
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Nebraska  (AWS,  1980).  However,  are  gmerally  a  hi^  altitude  siq>plement  to  (and  normally 
^nncictwit  witii)  flie  AIRMETs  given  to  militaiy  aircrews.  Both  forecasting  methods  are  described  further 
in  Chapter  2.  The  icing  AIRMETs  were  also  used  in  tire  icing  forecast,  numerical  model  algoritiun 
comparison  by  Brown  et  aL  (1996b),  as  a  standard-of-reference.  For  tiie  comparison,  it  is  important  to  note 
that  the  NWS  forecasters  issuing  tire  AIRMETs  have  other  supplementary  information,  such  as  satellite 
imagery,  and  current  PIREPs,  available  to  tiiem  to  enhance  tiieir  forecasts  and  this  additional  data  is  not 
taken  into  account  as  of  yet  by  present  automated  icing  algorithms.  In  fact,  AIRMET  forecasters  are  able 
to  use  the  present  AWC/NAWAU  icing  algorithms  as  forecast  guidance. 

1.4  Oveiwiew  of  Research  on  Aircraft  Icing  Variables 

Aircraft  certified  for  flight  in  known  icing  conditions  can  operate  for  only  a  limited  time. 
However,  they  may  not  be  able  to  cope  with  prolonged  exposure  or  with  extreme  encounters.  The  general 
aviation  community  has  been  particularly  vocal  in  calling  for  improved  icing  forecasts  and  for  renewed 
attention  to  this  hazard.  In  recognition  for  this  need,  the  office  of  the  Federal  Coordinator  for 
Meteorological  Services  and  Supporting  Research  (FCM)  issued  a  document  in  1986  addressing  these 
needs,  entitled  "The  National  Aircraft  Icing  Technology  Plan.  It  called  for  an  evaluation  of  existing 
techniques  and  the  development  of  improved  methods  for  detecting  and  forecasting  icing  conditions.  In 
1988,  the  Federal  Aviation  Administration  (FAA)  funded  the  National  Center  for  Atmospheric  Research 
(NCAR)  Research  Applications  Program  (RAP)  to  plan  a  multi-year  program  to  improve  aircraft  icing 
forecasting.  The  FAA  Icing  Forecasting  Improvement  Program  Experimental  Design  (Politovich,  1989a; 
Sand  and  Politovich,  1991)  was  completed  in  tire  Fall  of  1989,  and  the  6  year  FAA  Icing  Forecast 
Improvement  Program  began  in  October  of  that  year. 

In  the  past,  aircraft  icing  studies  have  taken  several  approaches  (Rasmussen  et  al,  1992).  Some 
of  the  aircraft  icing  environment  research  has  been  conducted  in  wind  tunnels,  where  die  nature  of  accreted 
ice  is  examined  on  different  airfoils  (including  Gayet  et  al.,  1984;  Hansman,  1985;  Hansman  and  Kirby, 
1986;  Shin  et  al.,  1991)  and  using  remote  sensing  devices  (Arking  and  Childs,  1985;  Popa  Fotino  et  al., 
1986;  d’Entremont,  1986;  Osborne,  1989).  Numerical  models  have  simulated  accretion,  predicted 
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supercooled  liquid  water,  and  determined  die  effects  of  tempaature,  liquid  water  content,  droplet  size,  and 
oflier  factors  on  icing  ^  and  intensity  (e.g.,  Newton,  1978;  Hansman,  1985;  Smart  et  al.,  1992).  Many 
in-flight  icing  studies  have  concentrated  on  the  analysis  of  environmental  factors  and  dieir  relation  to  ice 
accretion  in  bodi  continental  and  maritime  regions  (Bain  and  Gayet,  1982;  Hansman  and  Kirby,  1986; 
HofBnann,  1988;  Politovich,  1989b;  Hoffinann,  1990;  Politovich,  1993;  Cober  et  al.,  1995;  Fuchs  and 
Schickel,  1995)  or  its  performance  degradation  on  Ae  aircraft  (Cooper  et  al.,  1984;  Sand  et  al.,  1984,  Sand, 
1985).  Relatively  little  research  has  been  conducted  on  Ae  relationship  of  icmg-related  parameters  to  icmg 
observations  (mainly  pilot  reports)  routinely  available  to  operational  forecasters  (Sand  et  al.,  1984; 
Guttman  and  Jeck,  1987).  The  FAA  Icmg  Forecasting  Improvement  Program  chose  to  emphasize  Ae  latter 
studies,  using  Ae  NOAA  Forecast  Systems  Lab  (FSL),  and  Aviation  WeaAer  Center  (AWC;  formerly 
NAWAU)  as  test-beds. 

As  Ais  program  evolved,  it  became  clear  that  Aere  existed  excellent  opportunities  for  a 
collaboration  between  investigators  interested  in  Ae  iq)plied  topic  of  aircraft  icing  and  Aose  interested  in 
Ae  physical  processes  of  winter  storms  and  production  of  supercooled  cloud  water.  It  was  Aen  decided  to 
form  a  larger  field  research  program,  known  as  WISP,  to  incorporate  a  broadening  interest  in  Aese 
variables.  The  Winter  Storms  and  Icing  Project  (WISP;  Rasmussen  et  al.,  1992)  was  designed  to  ftirAer 
Ae  understandmg  of  Ae  dynamical  and  microphysical  processes  leading  to  Ae  production  and  depletion  of 
supercooled  liquid  water  (SLW)  in  winter  storm  clouds  and  to  improve  icing  forecasts.  The  observational 
strategy  utilized  for  WISP  was  based  on  collecting  m-situ  and  remotely  sensed  dato  Aat  could  be  used  to 
mcrease  our  understanding  of  Ae  various  processes  active  during  wmter  storm  events.  Smce  Ae  AWC  was 
Ae  eventoal  recipient  of  Ae  product  unprovements,  WISP  research  focused  on  automatmg  techniques 
currently  used  in  AWC  forecasting,  and  in  developing  improvements  on  Aese  meAods  based  upon  Aeir 
research  results. 

From  WISP,  came  a  new  surge  m  research  in  Ac  area  of  aircraft  icing  and  production  of  SLW. 
Most  of  Aese  stuAes  have  focused  on  Ae  development  of  algorithms  Aat  can  be  applied  to  Ae  output  of 
national  or  regional-scale  numerical  weaAer  prediction  models  in  order  to  Aagnose  and  predict  aircraft 
icmg  conAtions  (e.g.,  Rasmussen  et  al.,  1992,  Schultz  and  Politovich,  1992,  Knapp,  1992,  Forbes  et  al.. 
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1993,  Tliompson  et  al.,  1995  and  1996b,  Mawitz  et  aL,  1996).  Availability  of  such  algorithms  would  be 
very  useful  for  pilots  and  air-trafBc  controllers  if  die  algoridims  can  be  shown  to  be  accurate  and  reliable. 
Currently,  these  algorithms  are  rather  unsophisticated,  being  based  simply  on  whether  modeled 
atmospheric  conditions  meet  specified  temperature  and  humidity  criteria  (Brown  et  al.,  1995a). 
Unfortunately,  determining  the  accuracy  of  aircraft  icing  forecasts  is  not  straightforward,  primarily  due  to 
limitations  in  the  availability  and  quality  of  data  to  verify  such  products.  Currently,  the  best  available  data 
for  verification  of  icing  regions  are  pilot  reports. 

Icing  forecast  verification  via  pUot  reports  (PIREPs)  is  hampered  by  several  deficiencies  that  will 
be  discussed  in  Cluqiter  3.  First,  PIREPs  are  a  function  of  the  number  of  aircraft  flying  through  an  area  at 
any  given  time,  place,  or  altitude.  Icing  conditions  will  not  be  reported,  even  if  they  exist,  if  no  aircraft  fly 
through  the  area.  Also,  PIREPs  are  made  at  the  discretion  of  the  pilot  They  are  made  according  to  the 
pilot's  interpretation  of  the  atmospheric  conditions,  how  it  affects  aircraft  performance,  and  his/her 
experience.  Additionally,  reported  intensities  are  ako  dependent  on  the  individual  aircraft  type,  so  that 
reports  are  not  always  consistent  with  the  existing  meteorological  conditions.  The  result  of  these  defects  in 
the  verification  procedure  is  that  icing  probably  occurs  more  often  and  over  wider  areas  than  is  reported. 
This  lack  of  positive  icing  reports  may  also  resuU  from  die  lack  of  pilots  deciding  to  fly  their  aircraft 
through  regions  for  which  they  have  been  warned  may  be  an  icing  direat. 

Other  recent  research  has  focused  on  the  emergence  of  new  ground-based  remote  sensing  and 
satellite  sensing  technologies  as  a  tool  to  better  detect  icing  conditions  and  predict  icing  environments  (e.g., 
Jones  and  Yonder  Haar,  1990;  Curry  and  Liu,  1992;  Schickel  et  al.,  1994;  Stankov  and  Bedard,  1994;  Lee 
et  al,  1994;  Kleespies,  1995;  Brandes  et  al.,  1995;  Lawyer,  1995;  Politovich  et  al,  1995;  Ellrod,  1996a). 
Several  of  diese  studies  will  be  summarized  in  Ch^ter  2.  Unfortunately,  routine  atmospheric 
measurements  of  relevant  icing  parameters  are  not  nearly  dense  enough  to  accurately  predict  temperatures 
and  LWC,  or  determine  SLW  distributions  within  clouds  on  a  large  scale.  Also,  measurements  or 
Actitnat(»g  of  cloud  droplct  size  are  simply  not  widely  available,  except  for  short,  field  research  flight 
experiments. 
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Research  continues  on  the  use  of  ground-based  microwave  radiometers,  satellite-borne  sensors, 
and  oflier  techniques  for  remote  detection  of  icing  environments.  Satellite  imagery  can  be  useful  for 
diagnosing  icing  /  no-icing  regions  when  combined  widi  a  three-dimensional  temperature  analysis  of  the 
cloudy  environment  Thermal  infrared  (IR)  imagery  can  detect  cloud-top  temperatures  and  imply  cloud 
height  when  not  contaminated  by  overlying  cirrus  clouds.  But  with  the  improved  imaging  package  on  the 
next  generation  of  GOES,  better  resolution  and  additional  spectral  channels  hold  hope  for  an  improvement 
forecasts  of  icing  conditions,  as  sensed  by  satellites. 

1.5  Research  Objective 

Aviation  forecasters  are  still  greatly  challenged  to  forecast  aviation  weather  hazards  for  large  areas 
of  atmosphere  with  few  detaUed  observations  and  little  guidance  mformation  from  which  to  work.  One 
encouraging  tool  for  these  forecasters  are  the  satellites  tiiat  orbit  our  planet.  With  the  advent  of  the  next 
generation  of  Geostationary  Operational  Environmental  Satellites  (GOES-8),  additional  spectral  channels 
with  improved  areal  resolution  are  available  for  use  in  a  wide  variety  of  research  applications.  One 
important  outcome  of  this  applied  research  is  technique  developments  that  will  lead  to  improved  forecasts 
of  hazardous  aircraft  flying  conditions. 

In  this  study,  GOES-8  multispcctral  satellite  imagery  data  (using  channels  1,  2,  and  4)  are 
manipulated,  combined,  and  enhanced  to  retrieve  clues  about  cloud  phase,  supercooled  cloud  liquid  water 
(SLW),  and  cloud-top  temperature  (CTT).  To  retrieve  information  on  doud  phase  and  SLW,  the  known 
characteristics  of  several  GOES-8  spectral  channek  are  used  to  investigate  the  cloud  top's 
emissivity/absorptivity/reflectivily.  Thk  information  k  used  to  determine  whether  aircraft  icing  is  likely 
wifliin  ftiese  clouds  dirough  a  simple,  cloud-top  analysis  scheme.  The  period  of  study  consists  of  known, 
significant  icing  events  that  occurred  in  April  1995.  The  purpose  is  to  infer  through  inspection  of  cloud- 
top  IR  rhannel  brightness  temperatures  and  visible  charmel  albedos  whether  the  cloud  is  an  icing  threat  to 
aircraft  by  comparison  to  known  hazard  regimes  for  icing  variables. 

The  goak  of  this  study  are  to  test  muhkpectral  techniques  based  on  the  improved  GOES-8 
spectral  package  and  satellite  imagery  that  can  be  used  simply  as  an  operational  forecasting  tool  for  the 
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puipose  of  identifying  regions  witii  aircraft  icing  potential,  and  the  analysis  and  short-range  prediction  of 
aircraft  icing.  ITiis  p^wr  will  also  discuss  some  results  after  iqjplying  these  techniques  to  known  icing 
events.  In  particular.  Chapters  2  and  3  will  go  into  more  deptii  on  icing  and  PIREPs.  An  overview  of  die 
iqjplicable  radiative  transfer  tiieory  for  satellite  sensing  /  imaging  of  clouds,  and  tire  GOES-8  satellite 
imager  witii  its  improvements  and  spectral  package  wUl  be  described  in  Ch^ter  4.  Then,  Chapter  5  will 
describe  the  data  collection,  study  metiiodology,  and  data  processing  (computer)  resources  used  to 
accomplish  the  study.  Lastly,  the  results  fi-om  the  PIREP  verification  analyses,  research  conclusions,  and 
suggestions  for  future  work  will  be  discussed  in  Ch^ter  6. 
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Chapter  2  ->  Aircraft  Icing 

Past  methods  of  forecasting  meteorological  conditions  favorable  to  aircraft  icing  depended  mainly 
on  predicting  cloudy  environments  (as  the  source  for  the  generation  of  SLW)  in  the  appropriate 
temperature  range  (between  +4  and  -40  degrees  Celsius).  However,  the  ice  accretion  process  has  been 
found  to  be  also  greatly  affected  by  temperature,  liquid  water  content  (LWC),  die  size  and  concentration  of 
the  cloud  droplets,  the  cloud  phase,  and  the  presence  of  shear  layers  in  the  cloud's  environment,  among 
others.  This  chapter  will  go  into  detaU  on  these  different  fectors  that  affect  the  ice  accretion  process.  The 
accurate  forecasting  of  meteorological  conditions  favorable  to  aircraft  icing  is  difficult  for  several  reasons. 
Most  importantly,  because  these  parameters  are  hard  to  measure  or  estimate  within  the  forecasting 
environment. 

2.1  Aircraft  Ice  Accretion 

The  aircraft  icing  process  is  controlled  by  two  physically  distinct  processes.  The  first  is  the  in¬ 
flight  collision  of  liquid  water  droplets  from  die  ambient  cloud  environment  to  the  surface  of  the  aircraft 
(typically  leading  edge  of  wings  and  fuselage).  Then,  once  the  droplets  have  impacted  the  airframe,  their 
freezing  is  mainly  controUed  by  die  droplet’s  size,  airaraft  velocity  and  sh^,  and  thermodynamic 
processes  (Hansman,  1985). 

2.1.1  Types  ' 

If  die  heat  transfer  from  the  surfece  is  sufficient  to  remove  all  of  the  latent  heat  of  freezing  of  the 

impinging  water,  then  the  droplets  will  freeze  on  impact  resulting  in  a  dry  ice  surfece.  The  ice  shape 
typically  protrudes  forward  into  the  airstream  and  is  commonly  reported  as  "rime  icing".  Rime  ice  is 
opaque  or  whitish  in  color  and  rough  in  texture  due  to  the  quick  freezing  of  small,  supercooled  water 
droplets  trapping  large  amounts  of  air  between  the  impacting  droplets,  making  it  also  quite  brittle  (AWS, 
1980).  Rime  ice  is  most  likely  encountered  in  stratiform  clouds,  but  may  also  occur  in  cumulus-type 
clouds,  with  temperatures  typically  between  0*C  and  -20»C  (but  as  cold  as  -40“C  in  cumulonimbus).  It  is 
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relatively  easy  to  remove  by  conventional  de-icing  methods,  even  thou^  it  distorts  tiie  airfoil  more  dun 
"clear  icing"  (USAP,  1992) 

When  die  heat  transfer  from  die  surfrce  does  not  remove  all  of  die  latent  heat  from  impacting 
droplets,  the  ice  surface  becomes  wet  and  is  reported  as  glaze  or  "clear  icing".  This  clear  or  glossy  glaze  is 
murh  like  that  which  forms  on  trees  as  freezing  rain  fells  to  the  ground.  Clear  ice  is  perhaps  the  least 
common  type,  but  the  most  dangerous  form  of  icing  for  aircraft  because  it  adheres  so  firmly  to  the  airframe 
(Politovich,  1993).  It  is  most  commonly  encountered  in  cumuliform  clouds,  freezing  rain  or  drizzle,  and 
envirotunents  with  high  liquid  water  contents,  large  droplet  sizes  (>40jun),  and  warmer  cloud  temperatures 
between  0®C  and  -10®C  (but  as  cold  as  -25®C  in  deep  convection)  (USAF,  1992).  Clear  ice  is  generally 
smooth,  but  can  be  rough  textured.  It  is  smooth  when  deposited  from  large  SLW  or  drizzle  size  droplets 
that  spread  back  and  freeze  slowly,  taking  die  sh^  of  die  leading  edge  on  which  diey  impinge.  If  the 
impacting  droplets  are  smaller  and/or  mbced  with  ice  crystals  or  graupel  (soft  hail),  die  clear  ice  takes  on  a 
more  rough,  irregular  look  that  builds  hom-like  sh^s  back  into  die  airflow.  This  is  sometimes  referred  as 
"mixed  icing"  and  can  also  occur  during  flight  through  clouds  with  environments  exhibiting  bodi  types  of 
icing.  The  most  severe  aircraft  performance  degradation  is  typically  associated  with  such  homed  ice 
formations  (Hansman,  1989). 

2,1.2  Growth 

The  physical  processes  which  govern  icing  are  distinctly  different  for  wet  and  dry  ice  growth.  In 
dry  growth,  the  ice  accretion  is  controlled  by  die  local  rate  of  Uquid  water  impacting  the  surface.  This 
macc  flux  of  Water  is  determined  by  the  droplet  trajectories  as  they  pass  through  the  flow  field  surrounding 
the  airframe.  Because  the  airfirame  is  not  permeable,  flow  streamlines  do  not  intersect  the  body  of  the 
aircraft  and  die  coUection  efficiency  is  typically  highest  at  the  stagnation  point  for  the  body,  while 
decreasing  downstream.  The  collection  efficiency  of  the  body  (airfrrame  or  wing)  is  a  strong  function  of 
droplet  size  and  body  geometry  (Hansman  and  Kiiby,  1986).  Large  water  droplets  will  tend  to  cross  die 
streamlines  more  readily  and  consequently  impact  more  easily  (leading  to  higher  collection  efficiencies) 
due  to  their  greater  inertia  with  respect  to  smaller  droplets.  Small  droplets  will  follow  the  streamlines 


10 


more,  resulting  in  lower  collection  values  with  most  impacts  at  tiie  stagnation  point  (Figure  2.1).  Also 
important  to  note,  smaller  airfiame  surfeces  (e.g.  leading  edges  of  wings,  tails,  and  propellers)  are  more 
efficient  droplet  collectors,  because  there  is  less  room  for  die  droplets  to  turn  prior  to  impingement.  Thus, 
die  relatively  tiiick  wings  and  canopies  of  larger  commercial  and  cargo  jet  aircraft  have  a  smaller  collection 
potential  than  tiiose  of  tiie  trimmer,  commuter  turbo-prop  or  recreational  aircraft.  Smaller  aircraft  also 
have  less  available  power  to  counteract  the  fli^t  degradation  from  icing  and  typically  little  to  no 
equipment  to  counteract  icing  effects. 


RIME  ICE 


Figure  2. 1  Typical  dry  (rime)  and  wet  (mixed/clear)  ice  shapes  (from  Hansman  1 989) 

During  wet  growtii,  the  ice  collection  is  determined  by  the  rate  at  which  the  latent  heat  of  fusion 
can  be  removed  fi^m  the  surfiice  of  tiie  airframe  or  the  existing  ice  surftice.  This  thermodynamic  heat 
balance  determines  whetiier  tiie  accumulating  ice  surface  grows  tiirough  clear/mwed  or  rime  processes. 
Figure  22  depicts  the  principle  modes  of  energy  transfer  associated  witii  an  icing  surface.  Heat  is  added  to 
the  body  surface  primarily  by  the  latent  heat  of  fusion  released  as  the  droplets  freeze,  but  also  include  the 
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aerodynamic  heating  based  on  the  aircraft's  velocity,  and  droplet  kinetic  energy  impacting  the  surface  to  an 
even  smaller  extent  Heat  is  removed  from  tiie  surface  primarily  by  convection,  and  to  a  lesser  degree  by 
.nibiimation  (for  a  dry  surface)  or  eviration  (for  a  wet  sur&ce)  (Hansman  and  Kirby,  1987).  Also,  heat 
is  absorbed  from  die  surface  as  tiie  supercooled  droplets  impact  and  warm  to  0°C.  The  factors  which 
primarily  influence  die  heat  balance  are  the  temperature  difference  between  the  surface  and  the  droplet 
environment  the  convective  heat  transfer,  and  dre  amount  of  impinging  liquid  water  mass.  For  mixed 
icing,  both  the  droplet  impingement  and  heat  transfer  mechanisms  play  important  roles  in  the  icing  process. 
The  actual  hazard  for  large  propeller-driven  aircraft  comes  from  die  reduction  in  aerodynamic  heating  at 
lower  airspeeds  and  their  longer  flights  at  lower  altitudes  more  conducive  to  icing. 


Figure  2.2  Modes  of  energy  transfer from  an  icing  surface  (from  Hansman  1 989) 
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2.U  Intensity 

Icing  intensity  levels,  as  reported  by  pilots,  are  based  on  die  pilot's  visual  assessment  of  how  well 
his/her  aircraft  can  deal  with  die  ice  accumulating  on  the  airframe.  Ihus,  severity  or  intensity  is  quite 
subjective.  Table  2.1  shows  the  different  icing  intensities  as  reported  in  icing  pilot  reports.  It  depends  not 
only  on  aircraft  but  also  on  die  pilot's  experience  in  icing  conditions.  A  more  serious  problem  for  the 
scientific  researcher  is  that  the  current  index  of  intensities  reflects  little  basis  in  meteorological  parameters 
giving  forecasters  little  scientific  basis  for  including  icing  severity  in  their  forecasts.  Sand  and  Politovich 
(1991)  have  developed  an  icing  severity  index  to  replace  the  current  scheme  of  "trace-light-moderate- 
severe."  During  WISP91,  a  preliminary  version  of  this  icing  severity  index  (ISI)  was  used  to  visually 
display  icing  severity  on  a  pilot's  briefing  display  terminal  to  gain  feed  back  on  the  method.  User  reaction 
was  favorable  and  more  research  has  been  planned 


Table  2.1  Icing  intensities,  types,  and  a  sample  pilot  report  (from  FAA  1 977) 


Intensiiy 

Airframe  ice  accumulation 

Pilot  report 

Trace 

Ice  becomes  perceptible.  Rate  of  accumulation 
dighdy  greater  dian  rate  of  sublimation.  It  is  not 
hazardous  even  diough  deicing/anti-icing  equipment 
is  not  used  unless  encountered  for  an  extended  period 
of  time--over  one  hour. 

Aircraft  identification,  locatii^  time 
(GMT) ,  intensity  and  type  of  icing,* 

Li^t 

The  rate  of  accumulation  may  create  a  problem  if 
flight  is  prolonged  in  this  environmert  (over  <»e 
hour).  Occasional  use  of  ddcing/antidcing  equip¬ 
ment  removes/prevents  accumulation.  ^  It  do«  not 
present  a  problem  if  the  ddcing/anti-icing  equipment 
is  used. 

aimUQe/rJi.a9  airuraiv 

Example  of  pilot* s  transmission: 
Holding  at  Westminister  VOR  1232. 
Li^t  Rime  Icing.  Altitude  six 
thousand,  Jetstar  IAS  200  let 

Moderate 

The  rate  of  accumulation  is  such  that  even  short  m* 
counters  become  potentially  hamrdous  rad  use  ci 
deicing/anti-icing  equipment  or  diversion  is  necessary. 

Severe 

The  rate  of  accumulation  is  such  that  deicing/anti- 
icing  equipment  fails  to  reduce  or  control  the  hazard. 
Immediate  diversion  is  necessary. 

•  Icing  may  be  rime,  dear^  or  mixed. 

Rime  ice:  Rough  milky  opaque  ice  fonned  by  the  instantaneous  freezing  of  small  supercowed 
water  dbroidets. 

Qear  ice:  A  ^ossy,  dear  or  translucent  ice  formed  by  the  relatively  alow  freezing  ot  large 
supercooled  water  droplets. 

Mixed  ice:  A  combination  of  rime  and  dear  ice. 
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12  Meteorological  Factors  CoBtribnting  to  Icing 
22.1  Temperature 

Air  temperature  is  one  of  the  most  important  icing  parameters.  T^picaliy,  flie  ambient  cioud 
temperature  range,  which  is  most  closely  associated  wiA  die  greatest  hazard  for  all  types  of  aircraft,  occurs 
between  O'C  and  -20®C  (or  273-253K).  Meteorologists  normally  work  with  the  ambient  or  environmental 
temperature,  but  pilots  and  aircraft  designers  use  the  concept  of  Total  Air  Temperature  (TAT)  to  include 
aircraft  velocity  (aerodynamic  heating)  effects.  The  TAT  is  the  temperature  at  the  stagnation  point  of  the 
aircraft  and  corresponds  to  the  ambient  temperature  plus  an  additional  temperature  rise  due  to  Ae 
deceleration  of  Ae  oncommg  flow,  which  can  be  significant  at  high  velocities  (as  much  as  30®  at  500  knots 
of  airspeed)  (Hansman,  1989).  This  boost  m  temperature  upon  Aoplet  impact  results  m  a  very  wide  range 
of  temperatures  (wiA  common  occunences  between  +4®C  and  -40»C)  for  which  icmg  can  occur  based 
upon  Ae  type  and  properties  of  Ae  aircraft.  Smee  Ae  airftame  surfece  temperature  can  be  lower  Aan  Ae 
TAT  to  Ae  rear  of  Ae  gtagnarinn  pomt,  Ae  normal  procedure  in  commercial  jet  aircraft  is  to  run  Ae  anti- 
iemg  equipment  at  TAT  values  between  +10®C  and  -10»C  m  Ae  presence  of  visible  moisture.  In  some 
cases,  aerodynamic  heating  is  just  barely  sufficient  to  prevent  ice  accumulation  on  Ae  leading  edge  of  Ae 
wmg,  but  insufficient  to  prevent  icmg  aft  of  Ae  leading  edge  (Aus  to  Ae  rear  of  Ae  anti-  or  de-iemg 
equipment),  if  Ae  mtercepted  droplets  flow  back  over  Ae  wing  surfece. 

An  experiment  by  Olsen  et  al.  (1984)  usmg  Ae  NASA  Uwis  Icmg  Research  Tunnel  demonstrated 
how  small  temperature  changes  can  affect  Ae  ice  Aiqies  formed  and  Ae  resulting  drag  coefficients  (shown 
in  Figure  2.3).  They  found  Aat  at  cold  temperatures  (<  -15®C)  rime  icmg  caused  drag  mcreases  of  2  to  3 
times  Ae  clean  values.  But  as  Ae  temperature  mcreased  above  -10®C,  Ae  drag  mcreased  sharply  wiA 
temperature,  to  a  peak  value  of  nearly  8  times  Ae  clean  drag  (Aown  m  FIGS  Hans89).  This  was  caused  by 
a  transition  from  dry  rime  growA  to  a  wet  mixed  growA  icmg  regime.  The  horns  characteristic  of  mixed 
growA  can  be  observed  m  Ae  high  drag  ice  accumulations.  A  relatively  small  temperature  change  can 
cause  transition  from  relatively  benign  rime  icmg  mto  Ae  dangerous  mixed  icmg  condition,  but  Ae  non¬ 
linear  dependence  of  icmg  mtensity  on  Total  Air  Temperature  make  temperature  alone  a  necessary,  but  not 
sufficient  predictive  variable  for  aircraft  icmg. 
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Figure  2.3  Effects  of  temperature  on  ice  shape  (from  Olsen  et  al  1984) 


22.1  Liquid  Water  Content  (LWQ 

Liquid  Water  Content  (LWC)  affects  die  intensity  of  icing  in  two  primaiy  ways.  First,  increased 
LWC  implies  more  potential  supercooled  water  and  larger  accretions  within  a  given  time.  Thus,  high 
LWC  reflects  a  greater  severity  of  the  icing  encounter  whh  a  shorter  duration  of  flight  needed  to  achieve  a 
higher  icing  intensity.  The  other  effect  of  high  LWC  is  to  cause  the  icing  to  transition  from  rime  to  mixed 
icing  due  to  die  higher  impingement  of  water  mass  on  the  airframe,  since  at  a  higher  LWC  the  threshold  for 
the  transition  will  occur  at  colder  temperatures  (Hansman  and  Kirby,  1987).  Another  factor  which 
complicates  the  icing  process  is  the  large  variability  of  LWC  within  the  cloud  and  for  different  types  of 
clouds.  Various  studies  (Pobanz  et  al.,  1994;  Politovich,  1989b;  Hoffinann,  1988)  have  found  from  aircraft 
flights  through  confirmed  icing  clouds  that  liquid  water  contents  can  vary  according  to  the  size  distribution 
of  die  droplets  (large  droplets  in  low  LWC  regions  and  small  droplets  in  high  LWC),  but  remained 
generally  much  less  dian  1  gm‘*  (usually  between  0.1  -  0.5  gm’*).  Interestingly,  several  wind  tunnel  tests 
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needed  to  greatly  increase  the  LWC  (from  0.5  gm"*  to  over  1  gm*»  )  to  get  greater  desired  icing 

accumulations.  This  leads  to  the  conclusion  that  other  fiictors  are  at  work  within  natural  clouds.  Sadly, 

present  drop  collection  instruments  and  methods  of  quantitatively  detennining  the  observed  LWC  for  large 
cloudy  regions,  needed  to  make  icing  forecasts,  are  not  consistently  accurate  or  currently  available  to  the 
operational  forecaster. 

2.2  J  Droplet  Sizes 

The  size  of  tire  ambient  supercooled  water  droplets  both  in  terms  of  the  Median  Volumetric 
Diameter  (MVD)  and  die  actual  sh^  of  Ae  droplet  size  distribution  can  be  important  to  Ae  icmg  process 
(Hansman,  1989).  MVD  is  defined  as  Ae  Aameter  for  which  half  Ae  total  water  mass  of  a  distribution  of 
droplets  is  contained  m  larger  or  smaller  sizes.  MVD  provides  a  good  characterization  of  Ae  collection 
efficiency  of  Ae  entire  size  spectrum.  Typical  values  of  MVD  for  non-precipitating  Aoplets  m  icmg 
clouds  average  near  20nm,  but  can  range  from  near  lOpm  to  several  hundred  microns  (e.g.  -400  pm  wiA 
drizzle  size  defined  by  rqrprox.  200pm  (0.2  mm)  to  SOOpm  (0.5  mm)).  As  will  be  described  later,  larger 
SLW  droplets  are  more  efficiently  accumulated  by  Ae  airframe  (Politovich,  1993)  (see  Figure  2.4).  Within 
a  typical  cumulus  cloud,  Ae  droplet  size  distribution  is  such  Aat  Ae  bulk  of  Ae  water  mass  for  Ae  cloud 
results  from  only  a  small  number  of  large  droplets.  Ihis  analysis  can  naturally  vary  by  cloud  type,  yet  if 
Aese  larger  Aoplets  are  found  to  reside  togeAer  m  regions  wiAm  Ae  cloud,  it  would  signify  a  much 
greater  hazard  for  aircraft  Aat  should  venture  Arough  ft.  The  MVD  and  Ae  droplet  size  distribution  seem 
to  determme  Ae  effective  collection  efficiency  for  Ae  aircraft.  It  is  not  uncommon  to  have  trace  or  no 
icmg  even  at  high  LWC  for  clouds  wiA  predominantly  high  numbers  of  small  Aoplets.  Likewise,  for 
clouA  wiA  smaller  numbers  of  large  droplets  at  relatively  low  LWC,  ft  is  possible  to  have  moderate  to 
severe  icmg.  Again,  owmg  to  Ae  large  variability  of  clouA,  Ae  uncertamty  m  Ae  limits  of  measurement 
of  MVD  and  droplet  size  distributions  m  clouA  is  still  a  topic  of  controversy  and  information  is  based 
mainly  on  idealized  cloud  models. 
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Figure  2.4  Collection  efficiency  versus  impacting  droplet  diameter  for  different  airfoil  stes  and 
airspeeds  (from  Politovich  1993) 

2.2.3.1  Environmental  Conditions  for  Large  Drop  Regions 

Several  recent  studies  have  concentrated  on  characterizing  die  production  environment  of  large 
supercooled  liquid  water  (SLW)  droplets  in  light  of  its  significant  hazard  potential  for  aircraft  icing 
(Cooper  et  al.,  1984,  Sand,  1985,  Politovich,  1989a,  Rasmussen  et  al.,  1992,  Pobanz  et  al.,  1994,  Marwitz 
et  al.,  1996)  They  found  that  the  meteorological  conditions  accompanying  large  drop  regions  (30-250pm 
diameter)  include  warm  cloud  top  temperatures  (>  -15®C),  weak  thermodynamic  instability,  usually  lower 
droplet  concentrations  in  vicinity  of  the  large  drops,  and  shallow  wind  shear  layers  within  cloud  and  near 
cloud  top.  In  adflitinn,  ample  moisture  and  time  must  be  available  for  their  growth  and  enough  vertical  lift 
must  be  present  to  condense  the  moisture  in  die  updrafts.  Warm  cloud  tops  are  favored  since  they  decrease 
the  probability  of  ice  initiadon,  which  in  turn,  would  deplete  cloud  liquid  water  through  depositional 
growth  and  riming.  A  minimum  of  thermodynamic  instability  leads  to  weaker  convective  updrafts 
providing  time  and  liquid  water  for  hydrometeor  growth  in  the  cloud.  Additionally,  it  tends  to  result  in 
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cloud-top  diy  air  entrainment  and  eviration  of  some  cloud  liquid  Lower  droplet  concentrations  while 
common  are  not  necessary  to  enable  the  fewer  but  larger  droplets  to  collide  and  coalesce  more  readily  than 
smaller  droplets.  Dynamically  unstable  shear  layers  at  cloud  top  may  result  in  Kelvin-Helmholtz  waves, 
die  entrainment  of  subsaturated  air,  and  turbulent  mixing  resulting  in  large  SLW  droplet  regions. 

Pobanz  et  al.  (1994)  have  made  a  compelling  argument  that  dynamically  unstable  wind  shear 
layers  of  small  depth  (generally  less  than  200  meters)  within  the  cloud  environment  (both  inside  and 
outside  of  clouds)  can  contribute  to  the  formation  of  higher  concentrations  of  large  SLW  droplets  (30- 
2S0pm  diameter)  in  regions  near  cloud  top  of  stratiform  clouds.  At  cloud  top,  this  vertical  wind  shear  layer 
(in  direction  and/or  speed)  may  cause  turbulence  and  Kelvin-Helmholtz  waves,  resulting  in  the 
inhomogeneous  mixing  (witii  dry  air  entrained  from  above  cloud)  leading  to  large  SLW  drops.  The 
presence  of  vertical  wind  shear  qipears  to  either  initiate  or  accelerate  the  collision-coalescence  or 
"supercooled  warm  rain"  (i.e.  without  ice  crystals)  droplet  growth  process  in  such  a  way  diat  has  not  been 
sufficiently  physically  explained  and  confirmed.  Unfortunately,  due  to  their  shallow  nature,  it  is  unlikely 
that  ground  based  detection  systems  such  as  wind  profilers  or  Doppler  radars  could  resolve  tiiem. 
Likewise,  they  may  be  beyond  the  resolution  of  numerical  forecast  models. 

2.2.3.2  Hazards  from  Large  Drop  Regions 

The  Winter  Icing  and  Storms  Project  (WISP)  was  designed  specifically  to  study  the  production 
and  depletion  of  supercooled  liquid  water  in  winter  storms  (Rasmussen  et  al.,  1992).  Large  drop  formation 
regions  received  special  attention,  as  flight  plans  were  designed  to  specifically  sample  these  environments 
by  quick  ascents  and  descents  through  the  entire  depth  of  die  storm  clouds.  Ironically,  drops  in  this  range 
of  interest  (30-250nm  diameter)  are  particularly  hard  to  detect  unambiguously  using  instrumentation 
typically  mounted  on  current  research  aircraft  (Pobanz  et  al.,  1994).  Forecasting  large  SLW  droplet 
regions  presents  a  challenge,  since  they  seem  to  be  confined  to  thin  shear  layers  such  that  accurate 

detection  is  beyond  the  capability  of  current  forecast  models. 

The  presence  of  large  droplets  widiin  the  cloud  distribution  can  result  in  additional  hazard  through 
increased  droplet  impingement  over  the  aircraft’s  designed  limits.  Current  FAA  design  guidance 
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Tecommends  fliat  a  diameter  of  40  mioons  be  used  to  detennine  impingement  limits  for  aircraft  safety. 
Since  Ae  droplet  impingement  on  an  airframe  increases  with  droplet  size,  flie  presence  of  significant 
numbers  of  droplets  in  excess  of  40  microns  can  result  in  ice  accretions  occurring  behind  protected  regions 
of  the  aircraft  (aft  of  the  anti-  or  de-icing  equipment)  (Hansman,  1989).  Cooper  et  al.  (1984)  found 
anomalously  high  performance  degradations  when  flight  passed  tiuough  a  cloud  of  droplets  with  diameters 
fi-om  40  to  300  microns.  Marwitz  et  al.  (1996)  states  that  the  mere  presence  of  large  droplets  with  MVD  > 
30pm  was  not  suffreient  to  provide  extreme  decreases  in  flight  capability.  In  fact,  LWC  of  >  -0.2  gm 
were  also  needed  to  produce  severe  flight  degradation,  since  lower  LWC  wifli  MVD  >  30  pm  produced 
performance  loss  similar  to  that  fr-om  smaUer  droplet  encounters.  The  problems  resulting  from  large 
droplets  are  compounded  in  freezing  rain  or  drizzle  in  or  below  clouds,  where  botii  large  droplet  sizes  and 
high  LWC  are  combined.  Freezing  rain  frequentty  causes  clear  ice  accumulations  with  significant  run- 
back  icing  behind  otiierwise  protected  regions  of  the  airframe  (to  flie  rear  of  the  de-icing  boot).  This  is  an 
extremely  hazardous  situation  (to  be  avoided  by  pilots),  since  tiiese  conditions  can  extend  downward  to  flie 
earth’s  surface  and  leave  the  pilot  little  choice  of  evasion  if  at  low  altitude  (FAA,  1977). 


2.2.4  Cloud  Phase 

The  icing  potential  for  a  particular  cloud  is  directly  related  to  the  phase  of  the  hydrometeors  it 
contains.  Icing  normally  results  from  the  collection  of  supercooled  liquid  water  (SLW)  droplets  on  the 
airframe.  GeneraUy,  pristine  ice  crystals  do  not  adhere  to  the  aircraft  surfrices  after  impact  and  are  not 
alone  considered  a  hazard.  But,  if  the  ice  crystals  are  wet  from  vapor  deposition  or  partial  melting,  the 
impacting  crystals  can  stick  as  long  as  the  aircraft  surface  is  also  wet  due  to  recent  de-icing  or  flight 
through  a  region  of  high  LWC  (like  a  rainshower  or  drizzle).  However,  rarely  do  these  conditions  combine 
to  adversely  affect  the  performance  for  an  aircraft.  Though,  it  is  known  that  if  a  significant  number 
concentration  of  ice  crystals  arc  present  within  a  cloud,  the  chance  for  icing  will  be  greatly  diminished 
(Gayet  et  al.,  1984).  This  mainly  occurs  frwm  the  scavenging  of  water  vapor  by  the  ice  crystals  for  riming 
and  depositional  growfli.  Most  icing  forecasting  efforts  are  concentrated  on  identifying  regions  of  SLW  in 
clouds.  Techniques  are  available  to  predict  cloud  phase  in  stratiform  clouds  (AWS,  1980)  and  glaciation 
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usually  results  in  oveiestimation  of  Ihc  icing  severity.  However,  within  cumulifonn  clouds,  the  cloud 
phase  uncertainty  represents  a  potential  source  of  forecast  error,  since  ice  phase  and  SLW  phase 
environments  can  both  be  present  in  significant  amounts  wMiin  updraft  and  downdraft  regions. 

2.2.5  Cloud  Type 

Potential  icing  zones  in  the  atmosphere  are  highly  dependent  on  a  cloud's  structure  and  tire 
synoptic  or  mesoscale  environment  in  which  it  is  formed.  Icing  in  stable,  middle  and  low-level  stratiform 
clouds  is  typically  confined  to  a  layer  less  than  1  km  thick  and  generally  ranges  in  intensity  fiom  Trace  to 
Light-Moderate  witii  the  larger  accretions  occurring  in  tire  upper  portions  of  the  cloud  (near  cloud  top). 
Both  rime  and  mixed  icing  are  common  in  stratiform  clouds  owing  to  their  typical  make-up  of  a  narrower 
size  distribution  of  small  cloud  droplets.  The  main  icing  hazard  lies  in  die  great  horizontal  extent  of  some 
of  these  cloud  decks  and  the  possibility  of  a  long  flight  duration  through  them.  In  addition,  it  has  been 
shown  that  die  presence  of  wmditionaUy  unstable  l{q}se  rates  increases  the  potential  for  significant  icing, 
even  in  stratiform  cloud  systems,  by  triggering  embedded  convection  (Forbes  et  al.,  1993).  High-level 
stratiform  clouds  are  composed  mainly  of  ice  crystals  and  present  little  chance  of  icing. 

The  zone  of  probable  icing  in  cumuliform  clouds  is  smaller  horizontally  but  greater  vertically  than 
in  stratiform  clouds.  Icing  conditions  are  considerably  more  variable  in  these  clouds  depending  on  the 
stage  of  development  of  each  particular  cloud.  Intensities  can  span  die  full  range  of  possibilities  and  icing 
types  are  usually  clear  or  mixed  with  rime  with  the  most  intense/severe  icing  in  the  upper  half  of  the  cloud 
and  near  the  fi^eezing  level  (USAF,  1992).  Aircraft  icing  rarely  occurs  in  cirriform  clouds,  although  some 
contain  small  numbers  of  SLW  droplets.  However,  aircraft  icing  has  been  reported  during  flight  through 
the  cirrus  anvil  tops  of  cumulonimbus  clouds  where  strong  updrafts  are  able  to  maintain 

considerable  SLW  at  rather  low  temperatures. 

2.2.6  Frontal  Zones 

About  85  percent  of  all  icing  conditions  reported  occur  in  the  vicinity  of  frontal  zones  either  in  the 
warm  air  above  or  the  cold  air  beneath  (AWS,  1980).  Usually,  the  greatest  horizontal  extent  of  icing  is 
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assofiatH  wMi  warm  fronts,  and  the  most  intense  icing  wifli  cold  fronts  within  the  unstable  cumuliform 
clouds  (Figure  2.5a  and  2.5b).  Also,  incidences  of  freezing  rain  or  drizzle  are  more  common  around  warm 
fronts.  For  significant  icing  to  occur  above  a  frontal  surfece,  the  warm  air  must  be  fitted  and  cooled  to 
saturation  at  temperatures  below  freezing  producing  the  SLW.  If  the  warm  air  is  unstable,  icing  may  be 
sporadic.  But  if  it  is  stable,  icing  may  be  continuous  over  an  extended  area  with  studies  finding  significant 
mixed  and  clear  icing  up  to  200  miles  ahead  and  rime  icing  in  stratocumulus  up  to  300  miles  ahead  of  the 
surfece  warm  frontal  position.  Moderate  amounts  of  clear  icing  is  usually  limited  to  supercooled 
cumuliform  clouds  within  100  miles  to  the  rear  of  the  surfece  cold  frontal  position  with  increasing  intensity 
towards  the  frontal  zone.  Lighter  icing  is  otten  encountered  wWiin  fee  extensive  layers  of  supercooled 
stratocumulus  clouds  frequently  existing  behind  cold  fronts.  Anofeer  weather  situation  feat  is  sometimes 
conducive  to  icing  conditions  is  fee  deep,  cold-core  upper  level  low  pressure  systems,  when  not 
specifically  associated  wife  a  surface  frontal  system.  Ihis  situation  is  characterized  by  weak  to  moderate 
convective  activity  triggered  mainly  by  fee  greater  instabifi^  from  fee  colder  air  aloft.  Light  to  moderate 
mixed  icing  is  commonly  encountered  in  this  enviromnent.  Forbes  et  al.  (1993)  found  feat  nearly  52 
percent  of  fee  icing  cases  in  their  study  occurred  during  warm  advection  /  overrunning  events  versus  45 
percent  for  cold  advection  /  cold  air  aloft  events.  Interestingly,  they  discovered  cold  advection  /  cold  air 
aloft  regimes  accounted  for  83  percent  of  fee  west  coast  (CA)  cases,  but  only  12  percent  of  fee  central 
plains  (lA/MO)  episodes.  This  presumably  due  to  fee  dominating  local  effects  feat  are  characteristic  of  fee 

weather  systems  moving  through  feese  parts  of  fee  country. 

Marwitz  et  al.  (1996)  has  suggested  anofeer  preferred  location  for  fee  development  of  large 
droplet  SLW  within  a  mature  mid-latitude  cyclone.  It  is  fee  area  known  as  fee  "dry  slot"  in  fee  comma 
cloud  feat  is  often  evident  from  satellite  imagery.  They  hypofeesize  feat  fee  subsidence  or  advection  of  dry 
air  can  remove  fee  predominantly  snow  crystal  laden,  upper  level  clouds  feat  would  otherwise  r^idly 
deplete  fee  SLW  through  its  seeding  from  above.  Under  optimal  conditions,  fee  temperature,  moisture, 
and  wind  profiles  in  this  region  could  readily  combine  to  create  conditions  conducive  to  large  droplet  SLW 
formation  via  fee  vertical  wind  shear  process  described  previously. 
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Figure  2.5a  Cold front  icing zoik  (from  USAF 1992) 
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12.1  Vertical  Updrafts  and  Turbulence 

In  past  forecasting  methods,  vertical  updrafts  within  clouds  have  been  used  to  help  infer 
icing  environments.  It  is  ftiou^t  tiiat  diese  updrafts  would  support  die  formation  of  distributions  of  large 
SLW  droplets  by  holding  them  aloft  and  allowing  them  to  grow  ditough  collision-coalescence  processes. 
The  AWC  (formerly  NAWAU)  has  typically  used  vertical  velocity  information  at  the  700mb  level  fixim  the 
numerical  weather  prediction  models,  in  conjunction  widi  model  estimates  of  temperature  and 
relative  humidity  to  infer  clouds  with  greater  probabilities  of  icing  conditions  (Schultz  and  Politovich, 
1992). 

Small  scale  ambient  turbulence  has  recentfy  been  hypotiiesized  to  potentially  assist  in  the  icing 
process  (Hansman  and  Kirby,  1987).  This  fine  scale  (centimeter  or  smaller)  turbulence  is  believed  to 
strongly  influence  the  convective  heat  transfer  fixim  die  icing  surface.  Flight  test  observations  showed  a 
wide  variability  in  the  transition  threshold  between  rime  and  mbced  icing.  The  variability  was  thought  to  be 
due  to  variations  in  die  heat  transfer  resulting  fitim  die  ambient  turbulence  level.  While  this  effect  has  not 
been  directly  shown,  it  may  be  an  additional  source  of  uncertainty  in  die  process  or  it  may  be  a  function  of 
the  work  done  by  Pobanz  et  al.(1994)  described  above. 

2.2.8  Seasonal  and  Geographical  Influences 

A  ronrinental  U.S.  suTvey  of  USAF  aircraft  accidents  attributed  to  icing  over  the  previous  30 
years  (AWS,  1980),  showed  tiiat  95  percent  of  the  accidents  occurred  during  fall,  winter,  or  spring  with 
nearly  50  percent  during  winter  alone.  Those  accidents  that  occurred  in  the  summer  were  all  due  to 
carburetor  icing  (discussed  below)  or  flight  through  strong  thunderstorms  (which  normally  implies  higher 
intensities  of  icing).  This  should  be  no  surprise  since  these  are  the  seasons  where  most  of  the  previously 
discussed  ftictors  in  aircraft  icing  are  commonly  inesent  (more  cold/warm  fix>nts,  favorable  temperature 
and  humidity  combinations,  larger  cloud  systems,  etc.).  In  the  same  study  (also  confirmed  by  Cole  and 
Sand,  1991),  geographical  maximums  of  aircraft  icing  accidents  occurred  in  the  Great  Lakes  region  and 
similarly  in  the  Northern  Rockies  and  Plains  r^ion  (Figure  2.6).  These  results  tend  to  reinforce  the 
existence  of  previously  mentioned  factors  in  icing,  namely  the  available  moisture  to  make  clouds 
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Airframe  Icing  Accidents  (1975-88) 

(352  Airframe  Icing  Accidents  on  U.S.  Map) 


Figure  2. 6  Plot  of  aircraft  icing  accidents  (1975-88)  on  U.S.  map  (from  Cole  and  Sand  1991) 


(advection  from  die  lakes  or  Gulf  of  Mexico),  colder  temperatures  (higher  latitude,  lower  freezing  level, 
and  more  cold  front  passages),  and  the  availability  of  a  lifting  mechanism  to  form  large  cloud  regions  (be  it 
a  cold/wann  front  or  orogrsqjhic  lifting  from  the  hills  and  mountains).  Particularly  for  orogn^hic 
mountain  clouds,  the  effects  of  orogr^hic  lifting  can  be  evident  up  to  30,000  ft  or  more,  but  the  effects  are 
greatest  at  altitudes  between  5,000  and  18,000  ft  AGL,  since  the  majority  of  aircraft  involved  in  icing 
accidents  fly  at  these  altitudes  (Cole  and  Sand,  1991). 

Over  the  continents,  die  aviation  icing  hazard  generally  extends  frrther  equatorward  that  over  the 
sea  for  two  reasons  (Lee  et  al.,  1994).  First,  low-level  polar  air  masses  can  push  relatively  far  into  mid¬ 
latitude  continental  areas  widiout  wanning  to  temperatures  above  0®C,  which  contrasts  the  oceans  where 
cold  advection  is  quickly  neutralized  by  warming  below.  Second,  certain  orographic  structures  over  land. 
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as  radiatively-cooled  valleys  and  elevated  terrain,  can  be  associated  with  icing  conditions  even  when 
surrounding  areas  have  none.  The  potential  for  aircraft  icing  in  marine  environments  is  generally  confined 
to  high  latitudes  where  SLW  is  based  at  or  near  the  ocean  surface.  This  coincides  with  landing  and  takeoff 
routes  from  coastal  or  island  airports  (Curry  and  Lhi,  1992).  Additionally,  the  U.S.  Navy  is  particularly 
concerned  witii  the  low-level  icing  direat  for  aircraft  landing  or  departing  aircraft  earners. 


23  Aircraft  Icing  Forecasting 

The  AWC  provides  aircraft  icing  forecasts  (AIRMETs)  covering  the  continental  United  States 
every  six  hours  and  are  amended  as  conditions  warrant  (example  in  Figure  2.7).  The  icing  forecasts  are 
updated  as  SIGMETs  (Significant  Meteorological  Infonnation)  when  more  severe  events  (normally  a 
"nowcast"  situation  based  on  PIREP  information)  are  encountered,  or  AIRMET  reports  for  icing  issued. 
AWC  forecasting  procedures  are  generally  based  upon  guidance  from  numerical  weather  prediction  model 
graphics  distributed  by  the  NMC  (now  die  NCEP),  satellite  imagery,  radar  observations,  upper  air 
soundings,  synoptic  analyses,  and  PIREPs  (Tucker,  1983).  The  accuracy  of  this  forecast  is  then  largely 
dependent  upon  the  experience  of  the  forecaster,  his/her  skill  at  quickly  assessing  the  available  data,  the 
receipt  of  PIREPs,  and  the  accuracy  of  the  numerical  prediction  models.  Currently,  the  NMC  models  in 
use  do  not  explicitly  include  eidier  an  icing  hazard  or  a  cloud  liquid  water  forecast  (though  it  is  being  tested 
on  experimental  versions  of  die  Eta  model).  Therefore,  operational  forecasters  must  infer  icing  conditions 
from  other  information,  such  as  temperature,  humidity,  and  vertical  velocity.  Although  specific 
methodologies  differ,  the  AWC  forecasts  focus  mainly  on  larger  synoptic  scale  events.  Ihe  standard 
graphical  display  products  from  the  Eta  model  are  presently  being  used. 


2J.1  AWC  Method 

The  basic  steps  taken  by  die  AWC  to  infer  icing  are  to  forecast  die  fieezing  level  and  the  presence 
of  sufficient  level  at  altitudes  above  that  level.  After  the  areas  of  probable  icing  have  been  defined,  the 
forecaster  carries  out  a  more  detailed  analysis,  employing  empirically  derived  rules  to  establish  the  likely 
intensity  and  type  of  the  icing.  The  synoptic  situation  dictates  the  degree  of  additional  analysis  required  to 
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produce  die  icing  forecast  by  establishing  die  fiontal  zwies,  associated  clouds,  and  die  general  air  mass 
temperatures  for  each  storm  system.  Altitudes  widi  the  desired  icing  temperatures  can  be  inferred  from  die 
Eta  500-1000-mb  diickness  analysis  and  forecasts  with  diidcnesses  between  5220  and  5580  m  considered 
appropriate.  TTie  850-mb  temperature  field  (range  of -15  to  +4»C)  is  also  sometimes  used  as  a  guideline 
for  supercooled  cloud.  The  Eta  model  provides  relative  humidity  (RH)  information  for  the  500-1000-mb 
layer  and  experience  has  shown  that  a  threshold  of  RH  >  55%  for  this  layer  usually  implies  widespread 
cloudiness.  Finally,  the  forecasters  consider  vertical  motion,  which  is  provided  only  at  the  Eta  700-mb 
panel  (Schultz  and  Politovich,  1992).  Ibis  is  used  to  infer  a  condensation  of  the  available  moisture  into 
clouds  and  die  ability  of  the  atmosphere  to  hold  larger  cloud  droplets  aloft.  If  the  preceding  conditions  are 
met,  die  area  is  considered  to  be  an  aircraft  icing  threat 

23.2  AWS  Method 

In  die  Department  of  Defense,  die  U.S.  Air  Force  Air  Weather  Service  (AWS)  and  the  Naval 
Weather  Service  both  forecast  icing.  The  AWS  is  also  responsible  for  forecasting  for  the  Army  aircraft. 
The  Navy’s  techniques  are  very  similar  to  those  used  by  the  Air  Force  which  will  be  discussed  briefly 
below.  For  a  long  time,  die  definitive  document  on  icing  forecasting  in  both  the  military  and  civilian 
communities  was  the  Forecaster's  Guide  on  Aircraft  Icing  (AWS,  1980).  Since  its  publication,  little  has 
changed  except  for  die  general  improvement  of  and  heavier  reliance  on  die  prediction  models  on  which 
some  of  the  guidance  is  based.  It  goes  into  much  depth  on  the  factors  for  icing  listed  above  and  the 
methodology  of  the  icing  forecast  summarized  in  the  previous  paragraph.  The  guide  lists  numerous 
empirical  rules  and  techniques,  many  of  them  requiring  manipulations  of  upper-air  soundings  for  icmg 
threat  regions.  Unfortunately,  die  forecaster  has  precious  little  time  to  use  some  of  the  more  complex 
techniques,  while  trying  to  put  out  the  icing  forecast  in  addition  to  his/her  many  other  duties.  Militaiy 
aviators  obtain  pre-flight  weather  briefings,  including  Air  Force  Global  Weather  Central  (AFGWC)  icing 
products  (example  shown  in  Figure  2.8)  and  AWC  produced  AIRMETs  and  SIGMETs,  from  their  local 
weather  support  units.  These  icing  forecasts  give  an  integrated  picture  of  icing  conditions  in  time  and 
space  (both  horizontal  and  vertical)  along  the  flight  path. 
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Figure  2.7  Sample  AIRMETIcing  Forecast  Product 


Figure  2.8  Sample  AFGWC Icing  Forecast  Product 
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2.4  Aircraft  InHuences  on  Icing 

Several  aircraft  parameters  which  considerably  influence  in-flight  icing  are  discussed  below  in 
tenns  of  their  effect  on  the  physics  and  severity  of  the  ice  accretion.  As  has  been  previously  demonstrated, 
the  aircraft's  airspeed  (velocity)  effects  bofti  die  coUection  of  liquid  water  droplets  and  die  thennodynamics 
of  die  icing  process.  Greater  velocity  can  result  in  more  severe  icing  due  to  a  larger  potential  liquid  water 
exposure  by  increasing  both  the  path  swept  out  over  time  and  the  collection  efficiency  of  the  aircrafts 
surfeces.  But  ft  can  also  affect  the  heat  load  by  first,  bringing  more  liquid  water  into  contact  with  the 
airfiame  (requiring  more  latent  heat  to  be  dispersed),  and  increasing  the  stagnation  point  temperature  for 
higher  airspeeds  (requiring  a  colder  environment  to  off^  for  icing  to  occur  or  decreasing  the  transition 
temperature  between  rime  and  mixed  icing)  (Hansman,  1989). 

2.4.1  Structural  Airframe  Icing 

The  aircraft’s  shsqie  has  a  significant  effect  on  die  local  collection  efficiency  around  the  airfiame. 
Generally,  smaller  bodies  (as  measured  by  dieir  cross-section  into  die  airstream)  are  better  collectors  than 
large  bodies.  Therefore,  slender  components  of  the  aircraft  (e.g.  propellers,  fan  blades,  antennae,  and 
narrow  wing  and  taU  sections)  wUl  tend  to  be  most  sensitive  to  icing  accumulation.  In  fact,  because  of 
their  high  collection  efficiencies,  windshield  wqiers,  outside  air  temperature  probes,  pitot  tubes  or  other 
static  pressure  ports  are  frequently  used  to  visually  identify  and  estimate  icing  conditions  in  flight  (USAF, 
1992).  Ice  accretion  on  wing  and  tail  surfaces  disrupt  the  flow  of  rur  around  diese  airfoils.  This  results  in  a 
loss  of  lift,  an  increase  in  drag,  and  possibly  most  important,  causes  the  aircraft  to  stall  at  a  higher  airspeed 
than  normal.  The  weight  of  die  ice  deposit  presents  less  danger,  but  may  become  important  when  too 
much  lift  and  thrust  are  lost 

When  ice  forms  on  the  wings  or  the  fuselage  of  an  aircraft,  or  the  rotor  of  a  helicopter,  it  may 
seriously  endanger  die  aircraft.  Within  a  short  time  after  passing  dirough  a  region  of  SLW,  it  can  become 
f^AytahiUTPjt  by  ice  loading  At  times,  an  aircraft  can  climb  or  descend  to  exit  a  region  of  icing,  but  these 
strategies  are  not  always  effective.  An  ice-loaded  aircraft  may  be  too  heavy  to  climb  or  descent  may  be 
hazardous  if  the  supercooled  water  layer  extends  to  near  die  surface.  Experiments  have  shown  that  ice 
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accretions  of  only  one-half  indi  on  tiie  leading  edge  of  airfoils  on  some  airoaft  reduce  their  lifting  power 
as  much  as  50  percent  and  increase  drag  by  an  equal  amount,  thereby  increasing  tire  stalling  speed.  The 
serious  consequences  of  this  are  obvious,  since  one-half  inch  or  more  ice  accumulations  can  occur  in  only  a 
minute  or  two,  in  some  cases. 

Icing  on  the  propeUer  hub  and  blades  reduce  the  efficiency  of  tire  propeUer,  which  in  turn  reduces 
airspeed.  Increased  power  settings  may  then  fail  to  produce  sufficient  thrust  to  maintain  flying  speed,  all 
the  while  consuming  more  fuel.  An  even  grater  hazard  is  tire  vibration  of  the  delicately  balanced  propeller, 
caused  by  the  uneven  distribution  of  ice  on  the  propeller  blades.  Ibis  places  dangerous  stress  levels  on  the 
propeller  and  tire  engine  mounts  (USAF,  1992).  Tbus,  propeUers  with  low  revolutions-per-minute  (RPM) 
typical  of  smaller  recreational  and  business  aircraft  are  more  susceptible  to  icing  degradation  than  those 
witij  higher  RPM.  Tring  of  die  pitot  tube  or  other  static  pressure  pmts  is  dangerous,  because  it  causes  an 
inaccurate  indication  of  the  airspeed,  altimeter,  and  vertical  velocity  instruments.  Also,  ice  accumulations 

on  radio  anteimae  may  result  in  a  loss  of  radio  communication,  therefore  not  allowing  the  crew  to  request  a 

I 

change  of  altitude  or  course  to  fly  out  of  tiie  icing  region. 

2.4.2  Type  of  Aircraft 

The  effect  of  icing  varies  significantly  with  the  individual  aircraft  design.  Jet  aircraft  (including 
turbojet  and  turbofan  types)  arc  considered  to  be  tiie  least  susceptible  to  icing,  since  they  can  operate  with 
greater  quantities  of  excess  tiirust  which  can  be  used  to  offset  performance  degradation  from  icing. 
Additionally,  the  jefs  characteristic  flight  profile  is  to  quickly  climb,  cruise  at  higher  altitudes,  and  descend 
rapidly  through  the  lower  troposphere  where  icing  potential  is  greatest  (Hansman,  1989).  The  primary 
hazard  for  jet  aircraft  is  engine  frilure  due  to  Foreign  Object  Damage  (FOD)  usually  from  chunks  of  ice 
which  are  shed  from  otiier  sections  of  the  airframe  and  ingested  into  the  engine. 

Sensitivity  to  icing  conditions  is  quite  variable  for  propeller-driven,  general  aviation  type  aircraft 
(including  turboprops  and  smaller  reciprocating  engines).  Most  large  turboprops  have  been  designed  to  be 
fairly  ice  protected  and  can  operate  in  most  icing  conditions.  On  tiie  other  hand,  small  reciprocating  engine 
aircraft  are  generally  not  equipped  with  icing  protection  and  even  short  flights  tiirough  a  light  icing 
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environment  poses  a  great  hazard  to  the  aircraft  Most  propeUer-driven  aircraft  operate  at  lower  altitudes 
which  make  them  vulnerable  to  a  longer  duration  of  flight  through  regions  of  higher  icing  potential.  Icing 
tends  to  conventional  aircraft  most  Ihrou^  accumulations  on  die  propeller  hub  and  blades,  since  it 
directiy  leads  to  a  loss  of  thrust  (USAF,  1992).  Because  these  aircraft  are  already  using  the  bulk  of  their 
thrust  in  normal  flight  conditions,  there  is  little  margin  to  overcome  large  degradations  in  performance  due 

to  icing. 

Helicopter  icing  has  become  more  important  widiin  the  last  decade  where  flight  in  Instrument 
Meteorological  Conditions  (IMC)  has  become  more  fiequent  Helicopter  operations  occur  almost 
exclusively  at  low  altitudes  where  there  is  significant  icing  potential  for  both  in-cloud  and  sub-cloud 
regions.  AWiough  the  slow  forward  speed  of  die  helicopter  reduces  ice  build-up  on  the  fuselage,  the 
rotational  speed  of  narrow,  main  and  tail  rotor  blades  can  produce  a  r^id  growth  rate  on  certain  surfoce 
areas  (particularly  the  iimer  two-thirds  of  die  rotor  blades).  Rotor  icing  simultaneously  degrades  the  lift 
and  thrust  efficiency  of  the  heUcopter.  Furdiermore,  they  ^ically  operate  with  very,  slim  power  margins 
and  can  therefore  tolerate  only  minimal  ice  accretion.  Other  perilous  results  from  flight  through  icing 
conditions  include  the  blockage  of  engine  inlets  and  air  intake  screens,  and  vibrations  due  to  asymmetrical 
ice  loading  or  shedding  of  ice  from  the  rotors  (Fuchs  and  Schickel,  1995). 

2.43  Induction  (Carburetor)  Icing 

In  addition  to  hazards  to  structural  icing,  an  aircraft  frequently  is  subjected  to  icing  of  the  power 
plant  itself.  Ice  formation  can  be  found  in  the  air  induction  system  or  the  fuel  system.  While  beyond  the 
scope  of  this  study,  it  is  important  to  mention  these  types  of  icing  for  completeness  due  to  their  significant 
impacts  on  aircraft.  Carburetor  icing  is  most  treacherous  and  fi«quentiy  causes  complete  engine  failure.  It 
may  form  under  conditions  that  structural  icing  cannot.  If  the  relative  humidity  of  the  outside  air  being 
drawn  into  the  carburetor  is  high,  ice  can  form  inside  the  carburetor  in  cloudless  skies  with  the  temperature 
as  high  as  22“C  or  as  low  as  -10“C.  Carburetor  ice  forms  during  vqwrization  of  fuel,  combined  with  the 
expansion  of  air  as  it  passes  tiirough  the  carburetor.  The  temperature  drop  in  the  carburetor  due  to 
thermodynamic  cooling  can  be  as  much  as  40^,  but  is  usually  less  than  20'*C.  Ice  will  form  in  the 
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carburetor  passages,  If  cooling  is  sufficient  to  bring  the  temperature  inside  below  fteezing  and  with 
sufficient  moisture.  Most  aircraft  have  a  carburetor  heater,  an  anti-icing  device  which  preheats  die  air 
before  it  reaches  die  carburetor  keeping  the  air  above  fteezing,  which  can  be  turned  on  manually  if  flight  in 
humid  conditions  is  anticipated.  However,  the  heater  will  not  melt  any  ice  build-ups  if  they  have 
previously  formed.  Since  the  use  of  the  heater  affects  the  engine  performance,  it  must  be  used  sparingly 
arroHing  to  the  aircraft's  flight  manual.  Because  water  mixes  easily  with  jet  fuel,  the  fuel  absorbs 
considerable  water  when  the  humidity  is  high.  This  can  cause  fuel  system  icing  problems  when  the  fuel 
temperature  is  at  or  below  freezing.  In  addition,  structural  icing  on  smaller  air  intake  ducts  and  inlet  guide 
vanes  can  build  up  enough  to  restrict  airflow  to  the  engine  or  pieces  of  die  ice  can  break  away  and  be 
ingested  by  die  engine  causing  severe  damage. 

2.4.4  Results  of  Aircraft  Icing  Accident  Study 

The  National  Transportation  Safety  Board  (NTSB)  determined  that  icing  was  a  cause  or  factor  in 
803  aircraft  accidente  between  1975  and  1988.  Cole  and  Sand  (1991)  conducted  a  statistical  study  of  these 
icing  accidents  and  the  orographic  and  meteorological  parameters  that  could  have  had  a  systematic  effect 
on  them.  Their  findings  reveal  some  obvious  trends  and  some  surprising  factors  which  help  to  emphasize 
the  potential  dangers  of  aircraft  icing.  A  summary  of  some  of  tiieir  results  is  listed  below. 

Larger  concentrations  of  accidents  were  found  in  most  of  the  major  mountain  ranges  (53  percent) 
and  near  the  Great  Lakes  region  (14  percent)  (see  Figure  2.6).  These  findings  coincide  with  the 
geographical  factors  described  above.  It  seems  that  conditions  in  these  areas  can  better  increase  the 
moisture  content  of  the  atmosphere  through  lifting  of  air  masses  and  more  efficient  condensation  of  larger 
increased  amounts  of  moisture.  In  addition,  the  relatively  cooler  environment  can  lead  to  a  better 
production  of  supercooled  cloud  droplets. 

Accidents  were  equally  distributed  among  takeoff,  in-flight,  and  landing  of  the  aircraft  with  an 
average  of  nearly  65  per  year  during  tiie  study.  The  takeoff  cases  are  most  concerning,  since  they  mostly 
occurred  while  the  pilots  attempted  to  take  off  with  known  frost,  snow,  or  ice  on  the  airframe,  and  could 
have  been  prevented.  The  specific  effects  that  icing  have  on  the  type  of  aircraft  have  been  described 
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above.  Most  of  flie  aircraft  involved  in  aircraft  icing  accidents  in  die  study  were  the  smaller  general 
aviation  type  aircraft.  Few  of  which  were  equipped  for  flight  into  known  icing  conditions,  nor  did  they 

have  excess  power  to  deal  with  die  airfiame  ice  loading. 

Examination  of  die  pilot-in-command's  experience  indicated  that  icing  accidents  were  not  clearly 
attributable  to  the  less  experienced,  low-time  pUots.  The  average  pilot  had  over  2000  hour  of  command 
time.  Despite  this,  pUots  did  not  always  take  evasive  action  when  icing  conditions  were  encountered  in 
flight.  Also,  they  did  not  always  remove  ice  prior  to  takeoff  instead  sometimes  believing  the  snow  or  ice 
would  blow  off  or  not  affect  flight  performance. 

Weather  briefings  are  important  since  it  is  a  violation  of  FAA  regulations  to  fly  into  known  or 
forecasted  icing  conditions,  based  on  icing  intensity  and  type  of  aircraft.  Though  57  percent  of  the  pilots  in 
the  acridfpr  study  had  received  weather  briefings  and  icing  had  been  forecast  Forecasts  were  substantially 
correct  in  all  but  a  few  cases  according  to  the  NTSB  findings.  Shockingly,  43  percent  of  the  pilots  had  no 
record  of  receiving  a  pre-flight  weatiier  briefing,  which  is  also  a  violation  of  FAA  rules. 

One  surprise  was  that  fully  one-half  of  the  aircraft  accidents  were  attributed  to  carburetor  icing 
with  remaining  amount  structural  airfiame  icing  accidents.  Carburetor  and  induction  system  icing  as 
discussed  previously  can  occur  under  a  much  wider  range  of  conditions  than  for  structural  icing.  Although 
its  effects  are  fer  less  ^t  to  result  in  fetalities,  carburetor  icing  causes  dozens  of  accidents  each  year.  It  is 
mainly  an  engineering  and  operating  problem  ratiter  than  a  meteorological  one,  but  a  concern  to  aviation 

all  the  same. 

2.5  Past  and  Current  Studies  on  Icing  Detection  and  Prediction 
2.5.1  Field  Program  (WISP) 

One  goal  of  the  Winter  Icing  and  Storms  Program  (WISP)  was  for  the  detection  and  forecast  of 
regions  of  supercooled  liquid  water  (SLW)  cloud  droplets  which  would  lead  to  the  improvement  of 
forecasts  of  aircraft  icing  (Rasmussen  et  al.,  1992).  Comparative  verification  of  forecasting  systems  and 
algorithms  clearly  is  an  integral  part  of  die  process  towards  their  improvement.  Such  a  comparative 
evaluation  was  one  of  the  primary  goals  of  the  WISP  Real-Time  Icing  Prediction  and  Evaluation  Program 
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(WRIPEP),  wiuch  was  conducted  as  a  part  of  a  winter  field  program,  WISP94.  Additionally,  anodier  goal 
of  WRIPEP  was  to  develop  an  in-flight  icing  algoriflim  that  can  be  used  as  an  operational  forecast  tool 
(Brown,  1995a,  1996b).  WRIPEP  involved  a  real-time  display  of  icing  forecasts  based  on  several 
algorithms  and  numerical  weadier  prediction  models.  Model  data  and  associated  verification  data  were 
collected  and  archived  to  allow  an  in-depdi  comparison  of  Ae  models  and  algorithms  following  the  real¬ 
time  field  project.  More  details  regarding  WRIPEP,  the  models,  and  algorithms  can  found  in  Thompson  et 

al.  (1995). 

2SJ2  Icing  Algorithm  Research 

Many  studies  using  vaiying  methods  have  been  accomplished  to  try  and  address  this  aviation 
forecasting  problem  of  aircraft  in-flight  icing.  Some  previous  evaluations  of  algorithms  designed  to 
identify  and  forecast  aircraft  icing  regions  using  numerical  model  data  have  included  a  study  by  Schultz 
and  Politovich  (1992).  They  used  data  from  the  Nested  Grid  Model  (NGM)  to  formulate  and  evaluate  an 
algorithm  with  two  “threat  levels”  of  icing.  Forbes  et  al.  (1993)  expanded  on  this  work  through  the 
development  of  an  enhanced  algorithm  which  used  verification  results  to  optimize  the  Schultz-Politovich 
scheme  with  data  from  the  80-km  version  of  the  Eta  model.  Another  algorithm  was  developed  and  tested 
by  the  United  States  Air  Force  (Kn^p,  1992)  using  rawinsonde  data.  The  main  purpose  of  the  current 
algorithm  comparisons  was  to  aid  in  the  develoimient  of  a  better  automated  method  to  improve  on  the 
current  state  of  icing  environment  forecasting  for  commercial  and  military  aircraft  represented  by  Airman  s 
Meteorological  Information  (or  AIRMETs)  (Brown  et  al.,  1995a). 

2.5.2.1  Icing  Algorithms 

As  part  of  die  WRIPEP  for  Ihe  WISP94  period  (25  Jan  -  25  Mar  1994),  several  icing  algorithms 
were  evaluated  using  two  different  weaflier  prediction  models.  One  was  an  experimental  version  of  the 
National  fVntw  for  Environmental  Prediction's  (NCEP)  Eta  model  with  40-km  horizontal  resolution  and 
38  vertical  levels.  The  otiier  was  a  version  of  die  Mesoscale  Analysis  and  Prediction  System  (MAPS), 
which  is  also  run  by  die  NCEP  as  the  Rapid  Update  Cycle  (RUC)  model.  The  primary  icing  algorithms 
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evaluated  by  Brown  et  al.  (1996b)  are  listed  below.  In  addition,  cloud  liquid  water  (CLW)  estimates  by  Ae 
Eta  model  were  evaluated  to  sec  how  it  would  perform  as  a  prediction  variable  for  icing  environments  and 
supercooled  liquid  cloud  water  (SLW)  forecasts. 

The  RAP  algorithm  is  an  enhanced  version  of  Ae  algoriAm  developed  by  Forbes  et  al.  (1993) 
(also  included  in  Ae  comparison)  which  was  an  extension  of  Ae  Schultz  and  Politovich  (1992)  algorithm. 
The  RAP  algorithm  includes  four  different  categories  of  icing  (stable,  unstable,  freezing  rain,  and  freezing 
drizzle),  but  Ae  primary  criteria  for  Ae  yes/no  forecasts  are  based  on  temperature  and  relative  humidity 
thresholds.  The  AWC  /  NAWAU  algorithm  also  was  derived  from  Ae  Schultz-Politovich  scheme,  wiA 
improvements  and  adjustments  based  on  Aeir  operational  experience.  Like  Ae  Schultz-Politovich 
algorithm,  it  predicts  two  (threat)  levels  of  icing  wiA  Ae  second  level  having  more  restrictive  temperature 
and  humidity  criteria,  corresponding  to  a  greater  icing  Areat.  The  Air  Force  algonAm  was  developed 
using  rawinsonde  data  and  is  designed  to  predict  type  and  severity  categories  of  icing  (Knapp,  1992)  using 
pre-determined  temperature  and  relative  humidity  thresholds.  However,  for  WRIPEP  Ae  algorithm  was 
evaluated  only  in  terms  of  its  ability  to  predict  yesfoo  icing  conditions. 

2.5.2.2  Algorithm  Verification 

Murphy  and  Winkler  (1987)  describe  a  general  framework  for  forecast  verification  which 
encompasses  Ae  characteristics  of  Ae  forecasts,  Ae  corresponding  observations,  and  Aeir  relationship.  In 
particular,  Ais  framework  is  based  on  Ae  joint  Astribution  of  Ae  forecasts  and  observations,  and  on  two 
conditional  distributions  and  marginal  distributions  underlying  Ae  joint  distribution.  Then,  to  folly 
characterize  Ae  quality  of  a  forecasting  system,  all  of  Aese  distributions  must  be  known. 

Brown  et  al.  (1995a,  1996b)  goes  into  forAw  detail  regarding  icing  algorithm  forecasts  by 
describing  wheAer  a  binary  yes/no  icing  forecast  and  PIREP  verification  will  fit  Ais  fiamework  and  how 
Ae  ProbabiUty  of  Detection  (POD)  and  False  Alarm  Ratio  (FAR)  can  be  properly  used  wiA  PIREPs  in  Ais 
case.  Due  to  Ae  nature  of  PIREPs,  forecast  area  and  volume  were  used  as  measures  of  Ae  forecasted 
attribute,  since  it  is  difficult  to  estimate  overforecasting  parameters.  Because  Ais  framework  assumes  Aat 
forecasts  and  observations  are  collected  in  a  way  Aat  is  representative  and  consistent,  problems  arise  in  Ae 
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use  of  PIREPs  in  Ae  icing  observation  and  individual  verification  process  due  to  tiieir  relative  unsystematic 
nature  in  time  and  space  (much  unlike  tiie  gridded  data  fields  from  tiie  models).  (See  Chapter  3  for  more 
information  on  PIREPs)  However,  they  found  tiiat  a  POD  value  fiom  traditional  icing  verification  could 
be  used  to  advantage  in  a  relative  sense,  such  as  comparing  the  detection  rates  of  those  different  algoritiuns 
listed  above. 

2.5.2.3  Model  Comparisons 

Brown  et  al.  (1996b)  did  a  thorough  comparison  of  the  present  algorithm  methods  as  part  of  the 
Winter  Icing  and  Storms  Program  (WISP)  Real-Time  Icing  Prediction  and  Evaluation  Program  (WRIPEP) 
and  found  most  algorithms  to  be  adequate  and  promising,  but  none  superior  to  tiie  AIRMETs,  in 
forecasting  these  factors  or  atmospheric  conditions  and  tiie  implication  of  icing  regions.  The  verification 
results  for  the  explicit  CLW  output  from  the  Eta  model  suggest  tiiat  it  has  particular  potential  usefulness  as 
an  icing  predictor.  Differences  between  the  two  models  cm  which  tiie  algoritiuns  were  run,  showed  only  a 
slight  improvement  in  Probability  of  Detection  (POD)  witii  the  Eta  forecasts  over  tiie  MAPS  model.  The 
results  point  to  an  obvious  need  for  a  greater  physical  basis  in  the  icing  algorithms.  Other  comparisons 
have  been  accomplished  on  a  smaller  scale  by  Kietii  and  Cornell  (1995). 

2.53  Other  Recent  Important  Icing  Studies 

2.5.3.1  Using  An  Array  of  Ground-Based  Remote  Sensing  Systems 

Other  recent  approaches  to  determine  the  meteorological  factors  related  to  aircraft  icing  have 
focused  on  remotely  sensing  tiie  environment  (whether  from  ground-based  or  satellite  based  platforms). 
Stankov  and  Bedard  (1994)  used  an  array  of  ground-based  remote  sensors  to  successfully  predict  aircraft 
icing  conditions  aloft.  The  combination  of  profiling  radars,  radio-acoustic  sounding  systems  (RASS), 
multi-channel  microwave  radiometers,  and  lidar  ceilometers  were  effectively  used  to  define  tiie  spatial  and 
temporal  distributions  of  supercooled  liquid  water  (SLW)  which  would  be  used  to  infer  icing  conditions. 
The  remote  sensing  techniques  were  rqiplied  to  obtain  vertical  wind,  temperature,  and  humidity  profiles,  as 


35 


well  as  oti»er  physical  parameters  relating  to  clouds  (LWC,  water  v^r  content,  cloud-base  height).  Also, 
using  the  measurements  and  the  moist  adiabatic  tqiproximation,  a  cloud-top  height  was  inferred. 

While  successful  in  some  cases,  diis  system  of  remote  sensors  has  some  drawbacks.  Liquid  water 
falling  below  cloud  base  or  its  presence  on  the  radiometer  reflector  surfaces  causes  errors  in  its  retrieval  of 
moisture  or  temperature  profiles  and  precipitation  can  also  adversely  affect  the  lidar  ceilometers  by 
degrading  cloud  base  information.  In  addition,  die  vertical  limit  of  these  sensors  as  a  whole  is  around  3  km 
or  10,000  ft  AGL.  While  this  may  provide  valuable  icing  environmental  data  at  lower  levels  where  a 
majority  of  aircraft  fly,  it  cannot  be  very  successful  at  higher  altitudes.  Moreover,  the  cost  to  fully 
implement  their  system  of  sensors  would  be  enormous  if  it  were  required  to  cover  large  spatial  extents. 
This  horizontal  limit  (<  100  km^)  of  usefulness  makes  it  ideal  though  for  placement  around  the  major 
airports  with  increased  flight  Uaffic  requiring  ascent/descent  through  possible  icing  conditions  and  larger 
budgets  to  cover  the  cost  of  safety. 

Ideally,  Stankov  and  Bedard's  (1994)  suite  of  ground-based  sensors  could  be  used  in  conjunction 
with  satellite-sensed  information  in  order  to  provide  top-to-bottom  detection  for  the  suspected  icing 
environment.  The  fact  that  upper-level  cloud  cover  can  reduce  the  usefulness  of  the  satellite  to  detect 
cloud-top  temperatures  and  infer  SLW  of  lower  cloud  layers,  emphasizes  the  need  for  combining  a  range 
of  remote  sensor  information  (including  die  Doppler  radar)  to  receive  optimum  results.  However,  die 
emphasis  will  need  to  be  automation  of  each  of  die  analyses,  since  it  would  take  a  forecaster  quite  awhile  to 
timely  assimilate  all  of  die  information  into  a  usable  icing  forecast  product. 

2.S3.2  Satellite  Sensing  of  Cloud  Liquid  Water  Over  Oceans 

Satellite  remote  sensing  of  integrated  cloud  liquid  water  provided  the  basis  for  a  study  to  estimate 
aircraft  icing  regions  over  the  oceans  (Lee  et  al.,  1994).  They  used  the  Special  Sensor  Microwave  Imager 
(SSM/I)  85.5  GHz  channels  aboard  the  polar-orbiting  Defense  Meteorological  Satellite  Program  (DMSP) 
satellite  to  skillfully  redieve  images  of  low-level  water  clouds  over  the  northern  Pacific  Ocean.  While  an 
accurate  satellite,  integrated  cloud  liquid  water  parameter  has  yet  to  be  successfully  applied  over  large 
continental  regions,  die  microwave  sensor  has  proven  more  reliable  over  oceanic  environments  in  sensing  a 


36 


column  of  cloud  liquid  water  dirougb  Ihe  atmosphere,  owing  to  the  6ct  that  die  surfiice  emittance  is  feirly 
constant  for  water  surfaces  in  this  spectral  regime.  Using  die  integrated  cloud  liquid  water  (CLW) 
retrieval  images,  additionally  with  infiared  sensed  cloud-top  temperatures  and  numerical  model  output  for 
vertical  temperature  and  geopotential  height  structure  from  the  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS),  diey  developed  a  prototype  aircraft  icing  forecasting  system  that  could  be 
used  quickly  and  simply  by  4i  forecaster  to  assess  potential  icing  environments  over  oceanic  regions. 
Scientifically-based  icing  zone  temperature  (firom  O'C  to  -20®C)  and  CLW  (0.2  kg  m'*  or  greater) 
thresholds  were  specified  and  evaluated  as  potential  forecasted  watch  areas,  and  further  hone  the  accuracy 
of  the  retrieval.  Their  inclusion  of  satellite  data  helped  to  considerably  refine  icing  watch  areas,  as 
compared  to  those  fixim  the  Schultz-Politovich  forecast  algorithm. 

Unlike  ground-based  microwave  radiometry,  which  estimates  integrated  cloud  liquid  water 
continuously  at  a  single  geographical  point  (or  over  a  small  area),  the  SSM/I  gives  a  day  or  night, 
"snapshot"  view  of  integrated  CLW  over  a  much  wider  area  (-1400  km  swath),  but  unfortunately  only 
when  it  passes  over  die  area.  Though  ftie  method  can  produce  overestimates  of  icing  potential  when  ice 
clouds  cover  the  lower  water  clouds,  thus  blocking  the  infrared  (IR)  view  of  the  water  cloud-top 
temperature.  Even  known  errors  in  die  DMSP  IR  cloud-top  temperatures  of  1-2  K  introduce  an  additional 
source  of  error  into  icing  top  estimates.  However,  the  resulting  100-200m  height  assessment  errors  are 
small  compared  to  beneficial  improvements  in  icing  top  measurements  (usually  500  m  or  more)  when  the 
IR  method  is  incorporated.  Scattered  convective  clouds  can  cause  an  under  estimation  of  icing  potential 

owing  from  their  cloud-top  temperatures  below  -20°C. 

Precipitation-sized  droplets  (>  500pm)  can  also  cause  erroneous  retrievals  due  to  their  scattering 
properties  in  die  85  GHz  channel.  Since  die  presence  of  significant  precipitation  can  infer  some  amount  of 
supercooled  liquid  cloud  water  at  that  point  in  the  cloud  system,  Lee  et  al.  (1994)  tried  to  use  this 
information  to  dieir  advantage  by  spatially  averaging  the  integrated  CLW  through  the  precipitation  areas 
and  setting  a  direshold  of  greater  than  about  1  kg  m’*  to  indicate  precipitation.  Unfortunately,  due  to  the 
relative  scarcity  of  icing  reports  over  die  oceans,  they  were  unable  to  perform  a  validation  of  their  system. 
Their  method  is  most  appropriate  for  use  in  poleward  regions,  where  subfiwzing  temperatures  and  liquid 
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cloud  water  are  likely  to  exist  near  tiie  surfece.  Fortunately,  the  polar-orbiting  DMSP  satellite  provides 
relatively  frequent  coverage  in  these  regions.  Also,  at  high  latitudes  where  traditional  daytime  visible 
imagery  has  a  hard  time  distinguishing  low  clouds,  the  CLW  images  are  available  day  and  night  and  are 
free  from  sunglint  and  cirrus  that  tend  to  obscure  low  clouds  in  die  visible  channel. 


Satellite  Sensing  of  Cloud  Liquid  Water  Over  Land 

Cloud  liquid  water  (CLW)  has  been  shown  to  be  a  very  important  variable  to  the  study  of  aircraft 
icing.  With  the  abundance  of  commercial  and  military  air  traffic  flying  widiin  typical  altitudes  for  icing 
hazards,  mainly  below  15,000  ft  AGL,  over  the  continental  U.S.,  the  importance  of  icing  environment 
detection  and  forecasting  over  land  areas  is  much  more  important  than  over  most  oceanic  regions.  With 
most  research  concerning  satellite  microwave  remote  sensing  of  cloud  liquid  water  restricted  to  oceanic 
regions  (as  evidenced  by  Lee  et  al.  (1994)  above),  the  need  for  a  viable  method  to  detect  cloud  liquid  water 
over  large  regions  of  land  remained.  One  advantage  of  using  the  microwave  frequencies  lies  in  the  feet 
that  they  can  penetrate  thin  ice  clouds  or  cloud  tops  enabling  better  sensing  of  non-precipitating  convective 
clouds.  Again  however,  the  effective  use  of  satellite  microwave  remote  sensing  breaks  down,  when 

precipitation  size  droplets  are  encountered  within  or  below  die  cloud. 

Lawyer  (1995),  in  his  M.S.  thesis  written  at  Colorado  State  University,  developed  a  satellite 
derived  index  based  mainly  on  passive  microwave  imaging  data  from  Defense  Meteorological 
Satellite  Program  (DMSP)  polar  orbiting  satellites.  His  methods  followed  closely  those  of  Jones  and 
Yonder  Haar  (1990)  and  helped  to  retrieve  accurate  vertically  integrated  cloud  liquid  water  estimates  over 
continental  areas.  Similar  to  Lee  et  al.  (1994),  Lawyer  (1995)  also  used  the  high  frequency  85.5  GHz 
channel  of  the  SSMfl  instrument  aboard  the  DMSP  polar  oibiter.  Additionally,  visible  and  infrared  (IR) 
satellite  imagery  from  the  GOES-7  and  rawinsonde  information  were  used  to  supplement  the  microwave 
image  data. 

The  main  problem  for  accurate  microwave  CLW  retrievals  over  land  areas  is  resolving  their 
surface  emhtance  characteristics.  Over  water  surfaces,  the  emittance  from  the  surface  in  the  microwave 
portion  of  Ihe  spectrum  is  fairly  uniform  and  smaU.  This  feet  makes  the  CLW  calculations  based  upon  the 
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microwave  radiative  transfer  equation  much  simpler.  Retrieval  of  CLW  over  land  is  more  difficult  due  to 
die  greater  emittance  magnitude  and  variability  (mainly  from  soU  moisture  differences)  from  the  land 
surfrices,  which  causes  the  microwave  radiometer  to  be  insensitive  to  CLW  content  variations.  In  order  to 
retrieve  CLW  estimates,  both  SSM/I  and  GOES-7  IR  data  under  clear  sky  conditions  over  land  only  were 
used  to  make  die  surfiu»  emittance  calculations,  which  acted  as  a  surface  boundary  condition  for  the  total 
emittance  sensed  by  the  microwave  radiometer.  Lawyer's  resulting  CLW  values  were  consistent  with  those 
rnaHft  during  research  aircraft  experiments  under  similar  conditions. 

Then,  by  combining  seemingly  accurate  CLW  measurements  with  estimated  cloud-top 
temperatures,  inferred  cloud  depths,  and  cloud  texture  information  from  GOES-7  imagery.  Lawyer 
developed  an  icing  index  nomogram  as  a  proposed  icing  forecasting  tool.  A  positive  comparison  to  a  small 
set  of  icing  PIREPs  for  verification  showed  that  die  proqiects  for  a  microwave  CLW  retrieval  and 
multispectral  satellite  imagery  technique  for  icing  piediction  is  encouraging.  However,  further  testing  and 
automation  of  the  various  computational  tasks  widiin  the  mediod  must  be  completed  before  its  operational 
use  would  be  possible. 

More  importantly.  Lawyer 's  (1995)  research  has  shown  that  with  proper  knowledge  of  the  land 
surface  properties  (specifically  ground  moisture  and  surfrice  emissivity),  the  satellite-retrieved  estimates  of 
cloud  liquid  water  will  give  a  fairly  good  depiction  of  the  cloud  regions'  icing  potential  for  non¬ 
precipitating  clouds  (i.e.  stratiform  cloud  areas  and  some  cumuliform  clouds).  Especially  when  used  in 
conjunction  with  odier  remote  sensing  methods  in  the  future,  passive  microwave  imaging  can  be  a 
powerful  tool  for  the  researcher  or  forecaster.  However,  for  improved  spatial  coverage  and  more  timely 
forecasting,  a  dedicated  microwave  radiometer  placed  aboard  a  geostationary  satellite  would  allow  for  a 
more  complete  observational  system  (better  temporal  resolution)  and  increase  the  amount  of  information 
that  could  be  used  to  further  improve  icing  environment  detection  and  prediction. 


2ii3.4  Using  Multispectral  GOES-8  Imagery 

Anotiier  more  recent  sqrplication  of  satellite  remote  sensing  towards  detection  of  aircraft  icing 
regions  has  been  undertaken  by  Ellrod  (1996a).  Using  the  availabiUty  of  geostationary  satellite  imagery 


39 


from  Ae  advanced  GOES  I-M  series,  he  was  able  to  infer  toe  presence  and  spatial  extent  of  some  of  toe 
meteorological  fectors  important  in  aircraft  icing.  His  technique  uses  characteristic  information  garnered 
from  each  of  three  separate  GOES-S  spectral  channels  (CHI-Visible.  0.7pm;  CH2-Shomvave  IR,  3.9pm; 
CH4-Thermal  IR,  10.7pm)  to  diagnose  potential  aircraft  icing  zones.  As  they  pertain  to  this  study,  these 
channel  characteristics  will  be  further  mqjlored  in  Ch^ter  4.  EUrod  mainly  explored  stratiform  cloud 
systems  over  toe  northeastern  U.S.  for  known  significant  icing  events  (from  pilot  reports),  where  the 
reported  icing  intensities  ranged  from  light  to  severe  (although  most  were  of  moderate  or  greater  intensity). 

Digital  infrared  (IR)  brightoesk.  temperatures  and  visible  brightoess  counts  were  collected. 
Imagery  in  each  channel  was  also  studied  for  any  qualitative  information  that  would  help  to  detect  toe 
heavier  icing  regions.  As  noted  in  the  study,  when  displayed  in  toe  normal  IR  convention  (cold  shown  as 
white),  CH2  IR  often  appeared  significantly  warmer  (darker  gray;  higher  brightoess  temperature)  than  CH4 
due  mainly  to  the  respective  emission  and  reflection  characteristics  of  cloud  liquid  water  droplets  and  ice 
crystals  during  daytime  for  each  channel.  This  temperature  difference  between  toe  two  channels  (AT  = 
CH2  -  CH4)  was  utilized  by  use  of  an  empirical  scheme  toat  would  combine  toe  AT  and  the  visible  (CHI) 
channel  brightoess  counts  of  known  icing  cases  into  a  decision  tree  that  would  assess  icing  risk  in 
stratiform  clouds.  Ellrod  also  suggests  toat  using  a  visually  displayed  and  enhanced  combination  of 
channels  2  and  4,  it  may  be  possible  to  <^ct  icing  potential  for  clouds  in  toe  absence  of  obscuring  cirrus 
cloud  layers.  His  preliminary  results  show  some  skill  in  assessing  aircraft  icing  potential  by  combining 
data  from  the  three  imager  channels  listed  above.  WhUe  there  was  little  skill  in  icing  detection  or 
prediction  using  the  channels  individually,  the  CH2  imagery  revealed  some  success  at  highlighting  cloud- 
top  regions  where  heavier  icing  was  possible,  owing  from  its  increased  sensitivity  to  larger  liquid  cloud 
droplets.  While  tois  method  is  not  currently  a  stand-alone  solution  to  toe  discovery  of  aircraft  icing 
environments,  it  shows  promise  for  use  in  combination  with  improved  numerical  modeling  data  and/or 
other  ground-based  or  satellite-based  remote  sensing  techniques. 
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Chapters  — PUot Reports  (PIREPs) 

3.1  Description  of  PIREPs 

of  the  lack  of  physical  measurements  of  icing  conditions  in  the  atmosphere,  icing 
forecasts  have  primarily  been  verified  tom  subjective  observations  made  by  aircraft  pilots.  Pilot  reports 
(PIREPs)  can  provide  valuable  information  about  numerous  weather  conditions  observed  by  weraft  pilots 
while  in  flight,  during  landing  (descent)  or  taking  off  (ascent).  Some  examples  can  include  cloud  cover, 
fi-eezing  level,  air  temperature,  and  turbulence  and  icing  with  altitudes  and  intensities  as  determined  by  the 
pilot.  Unfortunately,  most  pilots  flying  conunercial  aircraft  do  not  include  die  ambient  air  temperature  in 
dieir  pilot  reports.  In  die  U.S.,  all  pilots  are  required  by  FAA  rules  to  report  hazardous  aviation  conditions, 
including  the  existence  of  icing  (Brown,  1995b).  The  pUots  voluntarily  issue  these  PIREPs,  as  soon  as 
aircraft  safety  aUows,  by  voice  over  the  radio  to  the  nearest  FAA  Flight  Service  Station,  though  sometimes 
at  die  request  of  an  air-traffic  controller.  The  reports  are  manually  entered  into  a  computer  as  part  of  die 
air  traffic  control  system. 

PIREPs  of  «<^tng  indicate  both  the  intoisity  and  type  of  aircraft  icing  (see  Figure  3.1). 
Severity/intensity  is  rated  on  a  scale  from  "Trace"  to  "Severe"  and  types  of  icing  include  Rime,  Clear,  and 
Mixed.  The  reported  severity  depends  on  a  subjective  assessment  of  how  well  the  aircraft  handles  die 
ice,  which,  in  turn,  depends  on  die  aircraflfs  deicing  or  and-icing  c^ability  while  in  flight 
This  variability  widi  aircraft  type  somewhat  diminishes  die  utility  of  die  icing  intensity  forecast 
information.  PIREPs  also  may  be  "negative",  indicating  diat  no  icing  was  observed.  However,  negative 
reports  are  rare,  since  the  pilot  will  usually  not  take  die  time  to  send  a  PIREP  stating  that  there  is  no 
significant  weather  hazard  unless  prompted  by  the  air-traffic  controller  (usually  during  takeoff  or  landing). 

PIREPs  have  some  drawbacks  diough.  Appropriate  to  the  voluntary  nature  of  PIREPs,  they  are 
biased  bodi  in  time  and  space.  This  represents  the  most  significant  problem  when  applying  them  to  a 
forecast  verification  analysis.  For  example,  very  few  PIREPs  are  recorded  at  night  due  to  the  relative  lack 
of  aircraft  Hnring  diis  time.  Moreover,  large  regions  of  the  U.S.  have  very  few  or  sparse  reports  (see 
Figure  3.2).  Larger  U.S.  cities  (defined  as  population  centers  of  250,000  or  larger  widi  a  major  airport  in 
the  vicinity)  with  correspondingly  greater  amounts  of  air  traffic  are  found  to  be  well  covered  with  PIREPs. 
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Figure  3 . 1  Example  pilot  report  (PIREP)  plot  with  icing  intensity  codes,  icing  type  symbols, 
and  icing  base  /  top  altitudes  (in  hundreds  of  feet  AGL) 


Figure  3.2  Spatial  distributions  ofPIREPs  in  1  “latitude/longitude  regions  (each  X represents  at  least  40 
PIREPs  over  a  two  year  period  (1990-1992)  (from  Brown  1993) 

Although  tiicse  large  cities  have  more  qiportunities  for  PIREPs,  many  parts  of  tite  U,S.,  particularly  the 
West,  have  few  large  cities.  Brown  (1995b)  found  tiiat  of  58,000  PIREPs  in  her  study,  33  percent  of  these 
reports  were  located  within  80-km  regions  around  large  cities.  These  fects  give  rise  to  a  bias  in  PIREP 
icing  reporting  mainly  during  daytime  and  around  large  cities  or  along  the  heavier  commercial  air  traffic 
corridors.  Additionally,  because  PIREPs,  overall,  are  not  systematic  in  time  or  space,  they  represent  a 
difficult  data  set  on  which  to  build  an  icing  climatology  (Brown,  1995b).  Other  problems  with  PIREPs 
result  from  errors  in  location  and  time.  The  pilots  often  can  pass  incorrect  information,  when  due  to 
aircraft  safety  concerns,  diey  cannot  take  note  of  the  proper  time  or  place  of  the  aviation  hazard,  and  were 
too  busy  to  pass  on  die  information  immediately.  Thus,  times  and  places  get  roughly  estimated,  and 
weadier  information  gets  entered  into  the  PIREP  system  archive  incorrectly. 

3.2  Methods  of  Scientific  Research  Verification  Using  PIREPs 

Since  PIREPs  do  not  give  a  representative  spatial  or  temporal  identification  of  the  icing  regions 
outside  areas  of  higher  air  traffic  volume,  they  are  to  be  used  very  carefully  for  icing  study  verification 
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purposes.  Tliese  aspects  of  the  statistical  use  of  PIREPs  arc  covered  weU  in  Brown  (1996a,  1996b)  and  are 
discussed  below.  Unfortunately,  at  this  time  tiiere  is  no  better  available  or  reliable  piece  of  data,  or  less 
expensive  method  tiiat  enables  die  verification  of  this  variable  to  such  a  large  spatial  extent  as  needed  by 
tiie  operational  forecaster. 

In  weather  forecast  verification,  the  verifying  data  typically  are  a  direct  measure  of  the  quantity 
that  is  being  forecast  with  the -verification  data  measiu^d  systematically  (with  one  obsen'ation,  yes/no,  to 
match  each  forecasted  occurrence).  It  is  important  to  note  that  the  non-occurrence  of  an  event  is  as 
important  to  measure  as  the  occurrence  of  an  event  For  such  weather  phenomenon  as  precipitation,  the 
non-occurring  events  are  systematically  recorded  in  time  and  space,  but  for  icing  events  this  is  generally 
not  the  case.  Although  negative  (no  icing)  reports  are  sometimes  recorded,  fliey  are  far  from  systematic, 
and  die  lack  of  reported  icing  in  an  area  cannot  usually  be  interpreted  to  mean  that  there  was  no  icing 
present.  In  feet,  the  small  relative  fiequencies  of  negative  PIREPs  do  not  accurately  reflect  the  large 
proportion  of  the  country  that  is  not  experiencing  icing  conditions  at  any  particular  time. 

Instead  of  being  recorded  at  regular  time  intervals  and  over  locations  of  interest  spread  equally 
over  large  spatial  domains,  pilot  reports  (PIREPs)  are  recorded  in  a  seemingly  haphazardly  fashion. 
Ironically,  the  distribution  of  PIREPs  is  actuaUy  quite  non-random.  Patterns  of  icing  reports  generally 
represent  flight  patterns  for  day-of-die-week,  time-of-day,  and  altitude  (Schultz  and  Politovich,  1992).  For 
instance,  air  traffic  is  normally  heavier  during  mid-week  than  on  weekends,  such  that  more  PIREPs  are 
received  on  weekdays.  Approximately  93  percent  of  the  1990-1992  PIREPs  occurred  between  1200  and 
0400  UTC,  which  corresponds  roughly  to  the  daylight  hours  over  the  continental  United  States.  Icing 
report  altit»4«i<!  for  that  same  period  were  mainly  concentrated  below  5  km  (-15,000  ft)  above  mean  sea 
level,  where  most  aircraft  are  likely  to  fly;  only  7  percent  of  the  PIREPs  occurred  at  altitudes  higher  than  5 
km  and  70  percent  were  between  1  km  and  3.5  km  (3000  and  10,000  ft)  (Brown  et  al.,  1993).  Tliese 
findings  suggest  that  it  would  be  best  to  restrict  icing  verification  analyses  to  tiie  daylight  hours  on 
weekdays  to  ensure  maximum  PIREP  representation.  However,  PIREPs  are  also  a  fimetion  of  the  severity 
of  the  weatiier  present  during  flight,  so  PIREPs  would  also  be  more  abundant  any  time  there  is  flight 
tiirough  significant  weather  systems. 
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Tbe  verification  of  icing  and  turtmlencc  finecasts  has  a  number  of  serious  pitfalls,  some  of  which 
were  noted  in  Brown  et  al.  (1995a,  1996b).  These  problems  are  primarily  related  to  the  nature  of  die  only 
data  diat  is  presently  widely  available  for  verification  of  aircraft  icing  events,  namely  the  pilot  reports 
(PIREPs).  Oiaracteristics  of  die  PIREPs  and  dieir  collection  have  important  implications  for  the  types  of 
verification  analyses  that  can  be  undertaken.  Moreover,  they  make  it  difficult  to  apply  techniques  based  on 
a  more  general  and  recognized  forecast  verification  scheme  based  on  work  by  Murphy  and  Winkler  (1987) 
and  further  applied  by  Brown  et  al.  (1995a,  1996a,  1996b).  Despite  their  drawbacks,  PIREPs  have  been 
used  to  estimate  standard  verification  measures  such  as  the  Probability  of  Detection  (POD)  and  False 
Alarm  Ratio  (FAR),  but  their  interpretation  is  not  straightforward  owing  to  die  nature  of  PIREP  collection. 
Brown  et  al.  (1995a)  concluded  their  verification  study  diat  estimates  of  POD  obtained  do  not  provide  an 
absolute  measure  of  the  ctqiability  of  die  algoridims  to  properly  forecast  icing  events,  but  POD  can  be 
useful  as  a  comparison  to  the  aqrabUities  of  different  icing  algoriduns.  Also,  estimates  of  FAR  were  found 
to  be  inappropriate  when  used  with  observation-based  measurements  (PIREPs)  ,  because  of  the  lack  of 
negative  reports  and  die  PIREPs  poor  coverage  of  the  forecast  grid.  The  inability  to  measure  FAR  in  this 

fashion  means  that  the  extent  of  overwaming  cannot  be  measured. 

In  Brown  and  Murphy's  (1996a)  pjqier,  a  procedure  is  described  that  may  possibly  overcome  the 
problem  of  ^grimating  the  icing  forecast  FAR  dirough  application  of  die  concept  of  "extending  die 
argument".  This  mediod  is  based  on  die  use  of  a  third  variable  (covariate)  on  which  to  condition  the  joint 
probabilities  of  forecasts  and  observations.  To  be  useful  in  this  context,  the  covariate  values  should  be 
strongly  related  to  the  observations.  For  an  icing  event,  such  a  covariate  might  be  cloud  coverage, 
supercooled  liquid  water  (SLW)  clouds,  or  integrated  liquid  water  as  estimated  fi-om  satellite  data  or  icing 
inferred  from  surface  remote  sensing  or  surfece  observations  of  precipitation.  One  distinct  advantage  of 
this  approach  is  that  the  covariate  can  be  a  variable  that  is  measured  at  the  forecast  valid  time  rather  than 
die  forecast  issue  time.  The  covariate  iqiproach  shows  promise  in  better  estimating  FAR  and 
overforecasting  of  icing  regions.  However,  the  actual  icing  region  is  unknown  in  real-world  cases,  which 
limits  the  abdity  to  realistically  determine  the  reliability  of  die  statistics  (Brown  and  Murphy,  1996a). 
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Chapter  4  —  Satellite  Background 

Satellite  imagery  seems  to  be  a  relatively  under-utilized  source  of  information  for  potential  icing 
detection  and  prediction.  When  compared  to  ground-based  sensors,  satellites  have  a  greater  capability  to 
accurately  measure  various  meteorological  parameters  over  larger  areas,  but  yet  with  sufficient  detaU  for 
jffpaiw  scales.  Some  work  on  detecting  icing  conditions  has  been  done  using  microwave  imagery  over 
oceans  from  polar-orbiting  satellites  (Curry  and  Lhi,  1992;  Lee  ct  al.,  1994;  Thompson  et  al.,  1996a). 
While  microwave  retrieval  of  integrated  cloud  liquid  water  from  satellite  has  been  relatively  restricted  to 
ocean  regions,  some  skiU  has  been  shown  in  summer  over  die  central  U.S.  (Jones  and  Vonder  Haar,  1990; 
Lawyer,  1995).  Widi  the  more  frequent  and  higher  resolution  imaging  data  attainable  from  die  next- 
generation  GOES-8/9  satellites,  die  potential  is  there  to  improve  die  icing  hazard  assessment  over  large 
continental  as  weU  as  oceanic  regions.  This  chapter  will  focus  on  die  radiative  transfer  theory  and 
desirable  characteristics  of  the  infrared  (IR)  and  visible  (VIS)  spectral  regimes  important  to  this  study.  It 
will  also  describe  the  GOES-8  imager  platform,  «q>plicable  characteristics  of  common  spectral  channels, 
and  the  improvements  made  overpast  meteorological  satellites. 

4.1  Passive  Visible  and  Infrared  (IR)  Remote  Sensing 

4.1.1  Atmospheric  Windows 

The  propagation  of  radiation  dirough  die  atmosphere  is  affected  mainly  by  two  processes: 
absorption  and  scattering.  During  absorption,  a  fraction  of  the  energy  passing  through  a  volume  element  of 
the  atmosphere  is  absorbed  by  atmospheric  constituents  and  re-emitted  at  different  wavelengdis.  Most 
absorption  is  spectrally  selective  by  major  gases,  such  as  COj  ,02,0,,  and  11,0  vapor,  with  lesser 
amounts  due  to  atmospheric  contaminants.  AdditionaUy,  a  fraction  of  die  radiant  energy  passing  through 
the  atmosphere  is  scattered.  The  amount  and  direction  of  this  scattering  is  largely  a  function  of  the 
scattering  element's  size  as  compared  to  die  wavelength  of  the  incident  radiation  being  scattered. 

Passive  imaging  sensors,  responding  to  radiation  in  atmospheric  windows,  are  very  important  in 
the  study  of  cloud  properties,  and  their  variability  and  distribution  in  space  and  time.  Atmospheric 
windows  in  the  infrared  and  visible  spectrum,  so-called  from  their  relative  insensitivity  to  certain 
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atmospheric  molecules  (water  vapor,  m  this  case),  are  used  so  that  the  tops  and  edges  of  low  clouds  as  well 
as  high  clouds  can  be  seen  (or  their  respective  cloud  particles  detected).  Moreover,  the  lack  of  atmospheric 
interference  in  these  wavelengths  will  ensure  that  the  emitted  radiation  being  sensed  will  be  dependent  on 
only  the  cloud's  microphysical  properties  or  the  Earth's  surfece.  For  example,  the  top  of  a  low  cloud  might 
be  only  a  few  hundred  meters  above  the  ground  leaving  about  99%  of  the  mass  of  the  atmosphere  between 
the  cloud  and  the  sateUite  sensor.  Thus,  energy  sent  from  the  cloud  towards  the  sateUite  would  be  very 
sensitive  to  attenuation  (reduction  in  returned  energy;  i.e.  absorption  or  scattering)  by  atmospheric 
molecules.  To  get  an  accurate  indication  of  the  cloud  properties,  die  amount  of  attenuation  by  the 
atmosphere  must  be  minimized,  and  transmittance  must  be  maximized. 

Atmospheric  windows  which  have  been  commonly  used  for  cloud  detection  and  are  relatively 
unaffected  by  water  vsqior  are  in  die  approximate  wavelengdi  ranges  of  0.4  -  1.3, 1.6  -  1.7, 2.1  -  2.3, 3.5  - 
4.0,  and  10  -  13  microns  (pm)  (shown  in  Figure  4.1).  The  spectral  response  of  satellite  imagery  is  usually 
narrower  than  diese  bandwiddis  since  it  is  desirable  to  avoid  some  molecular  absorption  bands  within  die 
water  vapor  window  and  also  to  take  advantage  of  variations  in  the  reflectivity  of  Earth  surfaces  (Bunting 
and  Hardy,  1984).  The  reflectivity  or  albedo  of  most  land  surfaces  in  the  visible  spectrum  tends  to  increase 
as  a  fimetion  of  wavelength  between  0.5  -  1.0  pm,  so  a  spectral  response  in  the  shorter  wavelengths  from 
0.55  -  0.75  pm  may  give  a  sharper  contrast  for  clouds  over  land  than  a  response  from  0.5  -  1.0  pm.  For  the 
shorter  wavelength  windows,  sunlight  (solar  radiation)  is  mostly  reflected,  so  tiiat  tiiey  are  mainly  useful 
during  the  daytime.  The  water  vapor  windows  at  longer  wavelengths  sense  thermal  energy  (temperature) 
originating  from  the  clouds  or  tiieir  backgrounds  with  only  small  contributions  firom  the  atmosphere  itself 
and  can  be  used  day  or  night.  The  3.5  -  4.0  pm  window  senses  radiative  contributions  from  both  reflected 
solar  and  emitted  tiiermal  radiation.  It  marks  the  transition  between  tiie  dominant  effects  of  reflected  solar 
(visible  /  near-infrared)  and  tiiermal  infrared  (IR)  regions  of  the  radiation  spectrum  on  Earth.  Thermal 
radiation  windows  are  also  found  at  longer  wavelengths  known  as  millimeter  or  microwave  fr-equencies. 
The  microwave  spectrum  roughly  extends  firom  1  mm  (1000  pm)  to  1  m  in  wavelength.  These  windows 
are  generally  more  sensitive  to  larger,  raindrop-size  water  particles  (>  500  pm)  within  or  below  clouds 
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raflicr  dian  tiie  entire  cloud  itself.  While  not  die  purpose  of  this  stu^,  these  fiequencies  are  commonly 
to  gadier  information  on  precipitation  and  cloud  liquid  water. 


Figure  4. 1  Visible  and  infrared  transmittance  versus  wavelength  showing  atmospheric  windows  (from 
Kidder  and  Yonder  Hoar  1 995) 


4.1J2  Radiative  Transfer  Theory 

Since  a  cloud's  characteristics  are  very  important  to  satellite  sensing,  the  interplay  between  the 
clouds,  and  solar  and  terrestrial  radiation  must  be  understood.  A  cloud's  optical  properties  such  as 
reflectivity  or  emissivity  depend  on  the  cloud's  physical  properties  as  well  as  the  wavelength  of  radiation 
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observed  by  a  satellite.  For  many  i^jplications,  it  is  important  to  distinguish  cloud  particles  from  the  larger 
precipitation  droplete  and  water  particles  from  ice  particles.  Though  greatly  variant,  typical  continental 
clouds  tend  to  consist  of  particles  on  the  range  of  5  to  300  microns  in  diameter  with  most  between  15-40 
pm.  Usually,  water-cloud  droplets  tend  to  be  much  smaller  tiian  prec^itation  and  ice-cloud  particles, 
although  sizes  vary  greatly  depending  on  the  cloud  type  and  the  particular  environment  in  which  the  cloud 
is  formed.  In  general,  pure  water  clouds  have  many  small  particles,  while  ice  clouds  have  fewer,  but  larger 

particles. 

Given  particle  sizes  and  number  distributions,  it  is  possible  to  calculate  the  radiation  sent  from  die 
cloud  towards  the  satelUte  (Bunting  and  Hardy,  1984).  The  calculations  start  with  the  amount  of  energy 
scattered,  absorbed,  and  emitted  by  individual  particles.  The  results  depend  on  the  refractive  index  and 
scattering  efficiency,  which  varies  wMi  the  wavelengdi  of  the  radiation  and  die  water  or  ice  phase  of  die 
cloud  particle.  Also,  the  results  hinge  on  the  radius  of  the  particle  (r)  and  the  wavelength  (k)  consistent 
wife  fee  relation, 

X  =  2nrlX 


where  x  is  known  as  fee  size  parameter. 

For  smaller  non-prec^hation  sized,  spherical  liquid  cloud  droplets  in  the  visible  and  infrared 
regimes  with  x  between  about  0.1  and  50  (where  droplet  size  is  nearly  equal  to  the  incident  wavelength), 
Mie  theory  determines  the  scattering  and  absorbing  characteristics  for  the  particles,  while  the  calculations 
are  less  straightforward  for  ice  particles  and  their  irregular  shapes  (Figure  42).  In  the  Mie  regime,  the 
angular  distribution  of  scattered  radiation  is  quite  complicated  and  varies  rapidly  with  %,  but  at  larger  x  (> 
10)  forward  dominates  over  back  scattering  with  litde  radiation  scattered  back  towards  the 

source  (Kidder  and  Yonder  Haar,  1995).  In  the  case  of  infrared  radiation  whose  source  is  frie  Earth,  this 
forward  scattering  preference  for  cloud  droplets  in  the  Mie  regime  results  in  more  terrestrial  radiation 
sensed  by  sensors.  Since  a  satellite  can  detect  billions  of  cloud  particles  in  every  pixel,  the 

radiative  intw^rrinn  between  diem  (multiple  scattering  across  the  droplet  size  spectrum)  must  be  included 
in  the  calculations  of  cloud-sensed  energy.  Because  all  visible  wavelengdis  are  scattered  nearly  as  well,  the 
cloud  appears  white  in  this  wavelength.  For  larger  precipitation-size  particles  and  drizzle  droplets  in  the 
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visible  and  IR  regime,  and  most  cloud  droplets  in  die  visible  with  x  >  50,  geometric  optics  are  used  to 
the  reflecting  and  refiacting  radiation  at  the  scatterer's  surfece.  Again,  at  larger  size  parameters,  die 
phntft  fimction  shows  a  greater  inclination  for  forward  scattering. 


Figure  4.2  Scattering  regimes  (from  Kidder  and  Yonder  Hoar  1 995) 


4.13  Spectral  Properties  of  Clouds 

A  detailed,  quantitative  theoretical  study  of  die  radiative  properties  of  water  droplet  and  ice  clouds 
at  visible  and  infrared  wavelengths  in  water  v^r  windows  is  presented  by  Hunt  (1973).  The  results  of  his 
ralrniatinns  of  cloud  optical  properties  are  presented  below.  Ibere  are  three  processes  which  control  the 
incident  radiation  upon  a  body  (cloud):  Absorptivity  (a).  Reflectivity  (p),  and  Transmissivity  (t),  and  they 
can  vary  according  to  cloud  density,  diickness,  height,  and  phase  of  particles.  For  energy  conservation  of 
die  inrident  beam,  die  three  quantities  are  related  by  a  +  p  +  t  a  1.  Assuming  the  cloud  tops  are  in  local 
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Ihennodynamic  equilibrium,  KirchofPs  law  states  that  at  each  wavelength,  a  body  that  has  good 
absorptivity  (a)  will  have  an  equally  good  emissivity  (e)  or  a  =  e.  For  a  Planckian  radiator  or  perfect 
blackbody,  aH6  =  l,andp  =  x*0  (Schott  and  Henderson-Sellers,  1984).  For  opaque  objects  (optically 
diick  clouds),  the  transmission  goes  to  zero  and  o  +  p  *  1  for  objects  in  die  visible  wavelengflis,  where  die 
primary  concern  is  with  reflected  energy.  Then,  in  the  longwave  inftared  wavelengths,  e  +  p  =  1,  where 
self-emission  is  dominant.  For  water  clouds,  visMe  wavelengths  are  more  sensitive  to  optical  depth,  while 
shortwave  infrared  wavelengths  are  more  sensitive  to  the  cloud  particle  size.  A  good  synopsis  of  IR 

radiance  retrieval  theory  is  contained  in  d'Entremont  (1986). 

Assuming  die  water  clouds  to  be  horizontally  homogeneous,  and  optically  diick  (cloud  optical 
depdi  (8)  >  10  or  clouds  at  least  several  hundred  meters  thick)  or  poorly  transmitting  at  shorter  (visible) 
wavelengtiis,  diey  have  high  reflectivities  due  to  low  emission  /  absinption  in  die  visible  regime.  A  high 
albedo  from  die  clouds  due  to  more  solar  radiation  reflected  and  scattered  back  to  the  satellite  by  smaller 
(ice  or  liquid)  cloud  particles,  results  in  a  higher  brightness  count  from  the  satellite  sensor.  Larger  cloud 
particles  have  the  effect  of  lowering  of  the  sensed  albedo  sUghtly  in  the  visible  band  due  to  lesser  amounts 
of  multiple  scattering  interactions. 

In  the  near-re.  (12  -  1.3, 1.6  -  1.7,  and  2.1  -  2.3  pm)  and  diortwave  IR  (3.5  -  4.0  pm)  windows, 
water  clouds  do  not  behave  as  perfect  radiators.  The  shortwave  IR  band  marks  the  transition  of  dominant 
radiant  energy  from  the  reflected  solar  regime  to  the  longwave  thermal  infrared  region  of  the  Earfli's 
electromagnetic  (EM)  spectrum.  Although  the  absorption  by  Uquid  water  clouds  increases  in  these 
windows,  their  emission  is  much  less  than  in  the  thermal  IR  window.  However,  there  is  a  large 
contribution  to  the  emergent  radiation  field  due  to  reflection  by  the  cloud  top  at  the  shorter  IR  wavelengths 
that  is  absent  in  the  longer  IR  wavelengths.  Ihe  reflectivity  in  the  shortwave  IR  window  is  0.10-0.17  for 
optically  tiiick  water  clouds,  but  only  0.02-0.05  for  tiiin  clouds  according  to  Hunt  (1973),  but  more 
importantly,  the  reflectance  is  highly  dependent  upon  particle  size  and  phase  with  smaller  particles  (radii  < 
15  pm)  being  better  reflectors  than  larger  ones  (Arking  and  Childs,  1985).  When  the  sun  is  rising  or 
cAtHng  tire  optical  depth  for  tiie  incident  solar  radiation  is  relatively  high.  The  water  cloud  reflectivity  is 
thus  inrreasing  witii  decreasing  sun  elevation  angle. 
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The  emissivities  of  ice  clouds  are  less  than  diat  of  water  drq)lct  clouds,  as  a  consequence  of  die 
ice  cloud's  lower  water  content  and  much  larger  particle  sizes.  Water  clouds  also  will  reflect  more  energy 
than  ice  clouds  in  this  region  of  the  spectrum.  Thus,  more  radiant  energy  will  be  received  at  the  satellite 
sensor  ftom  water  clouds  with  respect  to  ice  clouds  during  die  daytime.  Also,  small  cloud  water  droplets 
are  more  reflective  dian  a  like  number  of  larger  droplets  (Kleespies,  1995).  Plus,  number  distributions  of 
small  droplets  within  a  cloud  are  generaUy  much  greater  than  that  of  larger  droplets,  which  would  tend  to 
multiply  this  effect  Kleespies  (1995)  also  found  that  for  water  cloud  droplets  with  radius  >  6  pm,  the 
optical  depth  of  die  cloud  was  not  a  significant  factor  in  die  cloud  reflectance  returned  to  a  satellite  sensor. 
These  conclusions  reveal  some  of  the  useful  qualities  of  the  shortwave  inflared  band  for  distinguishing 
water  and  ice  clouds  during  the  daytime. 

In  the  longwave  or  tiiermal  IR  (10  -  13  pm)  window,  both  water  and  ice  clouds  are  very  poor 
reflectors  or  scatterers,  but  excellent  emitters  /  absorbers.  The  emergent  cloud  top  radiance  is  dominated 
by  the  emission  from  the  cloud  layer  (p  «  0,  so  e  «  1).  Therefore,  because  the  cloud  behaves  nearly  as  a 
blackbody  in  this  wavelength,  die  radiance  emitted  by  the  cloud  tops  can  be  used  to  estimate  die 
temperature  of  the  clouds  by  means  of  the  Planck  function  (which  relates  temperature  and  wavelength  to 
emitted  radiance).  Some  ice  clouds  can  be  opticaUy  thin,  therefore  they  transmit  a  significant  amount  of 
energy  fr-om  beneath  the  cloud  and  emit  correspondingly  less  energy  from  die  cirrus  cloud.  In  turn,  they 
usually  appear  warmer  than  dieir  true  temperature,  due  to  "warmer"  radiances  being  transmitted  through 
the  thin  clouds.  The  energy  sensed  by  the  satellite  radiometer  is  not  a  reliable  indicator  of  cloud 
temperature  for  these  thin  cirrus  clouds,  and  dieir  effects  on  the  returned  radiance  are  not  negligible.  The 
radiative  properties  of  clouds  in  die  diermal  IR  are  feirly  sensitive  to  die  size  spectra  of  the  cloud  particles. 
Increasing  the  size  of  the  particles  has  the  effect  of  increasing  the  transmissivity  (fr-om  die  Earth  source 
below),  decreasing  the  reflectivity  (scattering),  and  increasing  the  absorptivity  /  emissivity  of  the  cloud 
layer. 
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A2  Satellite  Sensing 

4^.1  Passive  Satellite  Imaging 

SatcUite-based  imaging  radiometers  (imagers)  are  commonly  used  for  cloud  identification  and 
characterization  purposes.  These  sensors  are  designed  to  provide  broad  horizontal  coverage  of  areas 
beneadi  tbe  at  die  finest  resolution  possible  while  enabling  the  satellite  to  sense,  store,  or  transmit 

the  Hata  back  to  Earth.  Most  meteorological  satellite  imaging  devices  scan  the  Earth  back  and  forth,  taking 
radiance  measurements  at  discrete  time  steps  and  producing  a  series  of  digitized  elements  along  each  scan 
line.  This  process  is  repeated  until  the  entire  Eardi  or  a  portion  thereof  has  been  sensed.  Each  rectangular 
element  is  assigned  a  specific  location  wWiin  die  image  by  referring  to  a  line  (or  scan)  number  and  an 
element  (or  sample)  number  (Kidder  and  Yonder  Haar,  1995).  When  displayed  on  a  monitor,  each  element 
is  referred  to  as  a  picture  elemoit  or  pfacel.  Using  the  sensed  digital  value  (Le.  brightness  temperature  or 
count)  for  each  pixel,  a  color  or  gray-shade  can  be  assigned  and  displayed  on  die  image. 

A  satellite  instrument’s  resolution  is  defined  as  die  size  of  the  smallest  element  of  a  scene  that  can 
be  resolved  by  the  instrument  It  is  a  fimction  of  the  field  of  view  of  the  sensor  optics  and  its  distance  finm 
the  scene.  Also  known  as  the  imager's  footprint,  the  resolution  can  be  thought  of  as  the  dimensions  of  the 
rectangular  box  on  the  scene  from  which  radiation  is  gadiered  at  each  time  step  during  die  scan.  The 
resolution  is  usually  dependent  on  the  wavelengdi  (spectral)  band  or  channel  being  interrogated.  Present 
geostationary  sateUites  are  capable  of  1  km  visible  and  4  km  IR  ground  resolution  directly  under  the 
satellite  (at  the  satellite  subpoint):  Polar-orbiting  satellites  are  capable  of  1  km  or  better  resolution  in  both 
visible  and  IR  channels,  but  diey  are  often  limited  by  die  data  storage  capability  of  the  spacecraft.  It  should 
be  noted  that  for  individual  cloud  structures,  cloud  breaks,  or  other  bodies  smaller  than  the  pixel  size,  the 
sensed  radiance  is  integrated  from  die  object  and  its  background  into  a  single  response.  Hence,  the 
radiometer  is  unable  to  satisfactorily  detect  in  isolation  any  source  diat  is  less  than  the  resolution  of  the 
channel  (Brimacombe,  1981).  Each  pixel  sampled  by  die  radiometer  is  recorded  as  a  voltage,  which  is  then 
converted  to  a  digital  brightness  count  onboard  and  transmitted  to  die  ground  on  a  scale  from  0  to  255. 
Calibration  tables  are  used  to  convert  diese  numbers  into  radiance  values,  whereby  temperatures  can  be 
inferred  and/or  images  can  be  produced  using  a  gray  scale  from  white  to  black. 
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The  feet  feat  fee  Earth's  surfece  is  curved  leads  to  some  importaiit  effects  when  its  surfece  is 
viewed  from  satellites.  First,  fee  degradation  of  fee  horizontal  resolution  of  fee  as  a  result  of  changing  fee 
nadir  angle  imagery  (known  as  foreshOTtening)  is  accentuated,  since  most  sensors  scan  out  for  some 
distance  from  fee  satellite  subpoint  Also,  fee  data  require  corrections  at  oblique  viewing  angles  for 
changing  viewing  geometry  between  fee  sun,  fee  cloud,  and  fee  satellite  for  solar  channels,  and  for  the 
increasing  path  length  through  fee  atmosphere  for  IR  channels.  Additionally,  fee  satellite  tends  to  see  more 
clouds  at  oblique  viewing  angles,  since  it  looks  ferough  more  atmosphere  and  has  a  higher  probability  of 
encountering  a  cloud.  This  also  leads  the  sensor  to  confuse  sides  and  lower  cloud  layers  for  fee  tops  of 
clouds.  Lastly,  higher  clouds  when  viewed  at  oblique  angles  may  not  be  correctly  assigned  spatially  due  to 
parallax  errors.  Thus,  pixels  may  be  misrepresented  horizontally  and  in  height  Other  potential  errors  in 
gftnsing,  such  as  sensor  lag,  signal  interference,  and  normal  atmospheric  attenuation  or 
contamination,  may  also  introduce  non-trivial  errors  into  the  remotely  sensed  product 


4.2.2  Satellite  Spectral  Band  Characteristics 

At  first  the  spectral  bands  used  on  weather  satellites  included  only  the  visible  and  thermal 
(longwave)  infrared  (IR)  regions  of  fee  electromagnetic  spectrum,  since  these  windows  held  most  of  the 
obvious  cloud  information.  But  over  fee  past  few  decades,  the  operational  and  scientific  use  of  other 
bands,  such  as  the  near-IR  window,  shortwave  IR  window,  and  water  v^r  band,  has  become  more 
commonplace.  While  most  meteorological  satellites  now  make  measurements  throughout  the  visible  and 
infinred  regions  of  fee  spectrum,  others  have  additional  capabilities  in  fee  ultraviolet  and  microwave 
realms.  A  brief  description  follows  below  regarding  fee  more  commonly  used  satellite  charmels  and  some 
of  their  uses  that  apply  to  cloud  sensing.  The  reader  is  referred  to  other  detailed  guides  on  satellite  imagery 
interpretation  and  detection  of  atmospheric  phenomena  that  go  much  further  beyond  the  scope  of  this  study 
(Weber  and  Wilderotter,  1981;  Brimacombe,  1981;  Scorer,  1986;  Rao  et  al.,  1990;  Bader  et  al.,  1994; 
Kidder  and  Vender  Haar,  1995). 
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42J1A  Visible  Band 

Probdjly  the  most  commonly  and  longest  used  q)ectral  band  from  satellite  imagers  is  in  the 
visWe  portion  of  the  spectrum.  The  visible  channel  typicafry  includes  wavelengths  from  0.4  to  0.9  pm  and 
exists  within  a  water  v^  window  and  at  the  peak  of  die  solar  EM  irtadiance  spectrum.  This  channel 
gives  us  imagery  in  colors  (or  giay-shades),  contrast,  and  textures  that  are  most  recognizable  to  the  human 
eye  and  for  a  quick  assimilation  of  land,  water,  and  cloud  features.  Also,  the  visible  channel 

typically  offers  the  greatest  spatial  resolution  from  the  satellite  sensor.  The  main  drawback  is  that  its 
imagery  is  only  useful  during  the  dayli^t  hours  unless  the  sensor  is  specially  equipped  widi  a  low-light 
(moonUght  Ulumination)  capability  (Kidder  and  Yonder  Haar,  1995).  The  visible  brightness  or  albedo  of  a 
target  depends  on  three  factors  -  the  intensity  of  die  solar  radiation,  the  sun's  elevation  angle  or 
illumination,  and  the  reflectivity  properties  of  the  body  itself.  The  reflectivity  of  a  cloud  depends  primarily 
on  its  tiiickness,  but  also  on  the  nature  of  die  cloud  particles,  whether  ice  or  water  and  dieir  sizes. 
Sometimes  a  higher  diin  cloud's  brightness  mtty  be  enhanced  when  the  underlying  surface  has  a  high 
albedo  and  the  radiation  from  diis  surfiice  increases  the  net  radiance,  making  it  hard  to  detect  thin  cirrus 
layers.  However,  under  similar  illumination,  a  water  droplet  cloud  will  appear  brighter  than  an  ice  crystal 
cloud  of  comparable  thickness  due  to  the  greater  number  of  multiple  reflections  that  are  possible  in  a  liquid 
cloud,  especially  for  smaller  water  droplets  (Brimacombe,  1981). 

During  visible  data  processing,  die  image  is  displayed  such  diat  die  greater  the  returned  radiance, 
the  brighter  die  pixel.  Visible  imagery  can  assist  in  die  interpretation  of  high  and  low  cloud  layers,  since 
generally,  thicker  or  higher  clouds  will  appear  whiter  or  brighter  due  to  higher  reflectivities  (albedo). 
Visible  data  is  often  used  to  discern  texture,  and  dius  cloud  type,  frxim  the  tops  or  edges  of  clouds.  It  can 
also  help  to  detect  high  clouds  (cumulonimbus)  embedded  within  lower  cloud  decks  dnough  recognition  of 
their  shadows  during  low  sun  (elevation)  angle  periods.  Visible  imagery  is  good  for  locating  low  cloud  or 
fog  covered  areas  diat  don't  show  up  well  in  the  infrmed,  but  it  can  cause  problems  discerning  snow  or  ice 
cover  from  clouds  witii  little  contrast 
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4.2M  Longwave  IR  (LIR)  Band 

The  thermal  or  longwave  infiaied  (LIR)  window  channel  fiom  10.0  -  12.5  fun  is  generally  die 
most  commonly  used  portion  within  die  middle  IR  realm  (about  3  to  15  pm).  In  diis  band,  where  die 
terrestrial  radiation  spectrum  peaks,  reflectivity  is  negligible  and  the  emissivity  of  land,  water,  and  thick 
clouds  is  very  near  one.  What  little  atmospheric  absorption  occurs  in  this  channel  is  principally  due  to 
water  vapor  and  may  cause  die  radiating  surface  to  appear  sUghtly  colder  than  it  actually  is.  This  factor 
becomes  most  pronounced  when  the  path  length  is  long  or  where  there  is  a  high  water  vapor  content  along 
the  path  (Brimacombe,  1981).  After  a  sUght  correction  for  these  effects,  these  surfaces  can  be  assumed  to 
be  blackbody  radiators.  Brighmess  temperature,  also  known  as  the  equivalent  blackbody  temperature,  is 
often  used  in  die  inftared  regime  to  express  die  amount  of  radiant  energy  impinging  upon  the  satellite 
sensor.  However,  nntike  in  the  microwave  spectrum  where  radiance  is  simpfy  proportional  to  temperature, 
the  infrared  brightness  temperature  must  be  determined  by  inverting  die  Planck  function  rather  dian  by 
simple  division  (Kidder  and  Yonder  Haar,  1995).  Inversion  of  remotely  sensed  radiances  into  atmospheric 
or  cloud  parameters  is  never  a  simple  problem  when  accuracy  is  a  must  Calculations  must  be  made  of 
atmospheric  transmittance,  Mie  scattering  parameters,  and  of  the  radiative  transfer  for  a  variety  of 
conditions,  and  potential  instrument  calibration  errors  and  assumptions  in  the  methodology  must  be  taken 
into  account. 

Since  IR  radiation  can  be  related  to  die  emitting  temperature  of  a  body,  and  because  die 
troposphere  generally  cools  with  height,  using  the  thermal  IR  channel  can  be  used  to  give  us  valuable  data 
about  the  temperature,  thickness,  and  height  of  clouds.  Thus,  the  images  constimte  a  thermal  map  of  the 
eardi’s  surface  and  cloud  tops.  In  visible  imagery  the  contrast  between  clouds  can  be  poor,  whereas  IR 
imagery  can  easily  discern  higher  clouds  and  lower  clouds  ftom  their  temperature  differences.  Another 
important  aspect  of  infitired  imagery  is  die  ability  to  provide  information  on  clouds  during  the  night,  when 
visible  imagery  is  useless. 

While  received  fi’om  visible  wavelengths  are  produced  like  ordinary  photographic  prints, 

LIR  images,  as  well  as  aU  other  IR  channel  imagery,  are  generally  produced  and  displayed  in  negative 
(Scorer,  1986)  (see  Figures  4.3  and  4.4  for  a  comparison  of  univerted  and  inverted  IR  images). 
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Uninverted  GOES-8  CH4 
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Figure  4.3  Uninverted  IR  (CH4)  image  for  10  Apr  95  1745Z 


Inverted  (Normal)  GOES-8  CH4 
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Figure  4.4  Inverted  IR  (CH4)  image  for  10  Apr  95  1745Z 


These  "inverted"  jmagtts  result  in  highw  (warmer)  radiance  values  being  represented  by  darker  pixels, 
instead  of  by  brighter  pfacels  as  they  would  in  die  visible.  Used  for  die  sake  of  comparison  widi  visible 
imagery,  diis  convention  means  diat  clouds  or  surfece  features  ^ipear  daiker  where  diey  are  warm  and 
whiter  where  diey  are  cool.  Ihus,  the  hipest,  coldest  (mosdy  ice)  cloud  tops  and  snow  cover  (widi  low  IR 
radiances)  ^ipear  white  in  this  channel,  and  die  land,  water  bodies,  and  lower  Uquid  water  clouds  ^pear 
darker  varying  with  their  warmer  temperature  (higher  radiances).  This  presents  a  problem  for  detecting 
low  clouds  at  night,  because  their  radiating  temperatures  can  be  nearly  the  same  as  the  background  (land  or 
sea). 

Since  in  this  channel  die  radiometer  is  effectively  measuring  emitted  temperatures  in  die  field  of 
view  and  not  reflected  energy,  some  poor  reflectors  display  different  qualities  that  are  not  apparent  in  a 
visible  image  (Brimacombe,  1981).  For  ecample,  thin  high  cirrus  clouds  consisting  mainly  of  ice  crystals 
at  low  temperatures,  while  occasionally  transparent  in  die  visible,  are  well  distinguished  in  the  IR. 
Although,  they  may  rqipear  warmer  fipom  higher  sub-cirrus  radiances  penetrating  through.  Also,  sea 
surfaces  in  the  visible,  except  for  sun  glint,  are  consistentiy  very  dark,  but  show  up  as  various  dark  gray 
shades  due  to  temperature  variations.  Temperature  variations  that  occur  naturally  fi-om  day  to  night  and 
with  seasonal  and  latitiiriinal  changes  wiU  affect  die  overall  appearance  of  infiared  images.  When  breaks 
between  clouds  in  die  field  of  view  are  smaller  dian  the  resolution  of  die  sensing  instrument,  higher 
radiances  fiom  underlying  suifeces  will  be  combined  widi  the  cloud  top  resulting  in  an  incorrectly 
integrated  radiance  for  die  pixel.  These  errors  tend  to  favor  a  warmer  pixel  radiance  and  might  reduce  the 
contrast  between  the  surfece  and  cloud  or  infer  die  cloud  to  exist  at  a  lower  altitude. 

4.2.2  J  Shortwave  IR  (SIR)  Band 

More  recently,  the  use  of  the  shortwave  infiared  (SIR)  window  channel  (3.5  -  4.0  pm)  has  become 
an  important  tool  for  researchers  and  forecasters,  though  it  must  be  used  with  some  caution.  The 
jqipearance  of  die  shortwave  IR  image  is  complex  during  the  day,  because  it  includes  both  radiative 
contributions  from  reflected  solar  radiation  and  emitted  terrestrial  radiation  in  varying  fiactions  dependent 
upon  the  solar  illumination  and  die  surface's  (cloud,  land,  or  water)  spectral  properties  (Figure  4.5). 
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Surfeces  with  high  reflectivities  in  tins  band  (liquid  water  clouds)  are  observed  as  surfeces  vwfli  higher 
brightness  temperatures  (satellite  images  iqipear  darker)  than  tiie  actual  temperature,  because  a  portion  of 
the  sun's  radiation  is  added  (reflected)  to  the  near-blaclcbody  radiation  emitted  from  the  cloud  surfiice 

(Liljas,  1986)  (see  Figure  4.6). 

A  very  important  property  of  the  SIR  spectral  band  is  that  ice  is  very  absorptive  of  these 
wavelengths.  Snow  covered  ground  and  sea  ice  are  highly  absorptive  and  scatter  very  little  of  the  radiation 
in  this  band.  Clouds  composed  mainly  of  ice  crystals  could  absorb  more  radiation  than  those  made  up  of 
liquid  water  droplets  of  like  numbers  and  sizes.  But  more  importantly,  most  water  clouds  are  composed  of 
a  larger  distribution  of  much  smaller  particles  on  the  order  of  tiie  incident  wavelength  in  this  channel 
(Scorer,  1986).  Not  only  is  there  much  less  absorption  from  ^ical  ice  clouds  (cirrus),  but  the  additional 
solar  reflection  component  from  tiie  water  clouds  results  in  much  greater  radiances  received  by  satelhte 
sensors.  The  shortwave  IR  image's  interpretation  is  not  as  straightforward  as  each  of  the  previous 
channels,  since  it  cannot  quantitatively  establish  cloud-top  temperatures  like  tiiermal  IR,  nor  can  it  reveal 
all  cloud  features  as  easily  represented  by  visible  imagery.  This  band  tends  to  display  more  charactenstics 
typical  of  the  visible  regime  during  the  day  and  the  thermal  infrared  band  at  night. 

While  the  shortwave  infrared  channel  shares  much  in  common  with  tiie  two  previously  described 
spectral  bands,  its  imagery  is  somewhat  hard  to  interpret  witiiout  a  review  of  its  channel  characteristics  that 
are  listed  above  in  the  summary  of  Hunfs  (1973)  study.  Some  other  aspects  of  clouds  in  the  shortwave  IR 
band  follow  from  research  by  Aricing  and  Childs  (1985),  Scorer  (1986),  Liljas  (1986),  and  others.  Daytime 
fog  and  low  stratus  have  a  high  reflectivity  due  tiieir  large  distribution  of  small  water  droplets  and  appear 
darker  (warmer)  to  the  satellite  sensor  from  tiie  increased  radiance.  But  during  the  night,  they  appear 
slightly  cooler  (brighter)  in  tiiis  wavelengtii  versus  tiie  thermal  IR  channel  (d'Entremont  and  Thomason, 
1987;  Ellrod,  1995;  Nebon  and  Ellrod,  1996).  Thin  cirrus  (ice)  clouds  tiiough  cold  appear  relatively 
warmer  at  night,  especially  when  over  low  clouds,  due  to  their  high  transmissivity  and  resultant  upwelling 
radiance  contamination  from  below  (Ellrod,  1992).  Low  clouds  over  snow  cover  or  ice  on  water  bodies 
appear  daik  owing  to  tiieir  higher  reflection  and  absorption  combination  (Kidder  and  Wu,  1984;  Allen  et 
al.,  1990).  Developing  cumulonimbus  and  towering  cumulus  clouds  ^pear  warmer  (darker)  than  do 
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Figure  4.5  Schematic  showing  the  radiation  components  reaching  channel  2  of  GOES-8  imager  (adapted 
from  Liljas  1986) 


Figure  4.6  Brightness  temperature  difference  versus  time  for  SIR  and  UR  channels  on  GOES-VAS  (from 
Kleespies  1995) 
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decaying  cumulonimbus  (brighter),  because  they  are  mainly  liquid  water  topped  during  formation  and 
glaciated  when  mature  or  dissipating.  Sometimes  warm,  dark  cores  wiU  ^)pcar  in  a  developing  or  mature 
cumulonimbus  cirrus  shield  (Setvak  and  Doswell,  1991).  Ibese  are  likely  the  result  of  greater  absorption 
by  larger,  supercooled  Uquid  water  drops  and  reflection  by  larger  numbers  of  smaller  liquid  droplets  which 
have  been  pushed  aloft  to  extremely  high  and  cold  levels  by  the  strong  updrafts  within  the  storm. 
Detection  of  clouds  during  the  night  using  the  shortwave  IR  band  is  much  like  that  of  the  thermal  IR  band, 
but  with  not  nearly  as  strong  of  a  blackbody  return,  since  the  dominant  radiant  energy  feature  is  only  the 
emission  from  each  surface.  While  this  channel  can  be  used  to  identify  clouds  at  night  (Olesen  and  Grassl, 

1985;  d'Entremont,  1986),  it  is  not  used  within  the  scope  of  this  study. 

In  research  studies,  the  shortwave  IR  channel  has  been  used  alone  and  in  conjunction  with  other 
satellite  rhanimk  to  interrogate  many  atmospheric  and  Earfli  surface  phenomenon  (Ellrod,  1992;  Kleespies, 
1989).  Kidder  and  Wu  (1987)  and  Prins  and  Menzel  (1992)  have  used  toe  channel  to  detect  very  hot 
surfece  temperature  regions  such  as  industrial  hot  spots,  volcanoes,  mass  burning,  or  forest  fires.  Use  of 
toe  SIR  and  LIR  channels  through  a  channel  differencing  scheme  led  to  an  improved  ability  to  detect 
nighttime  rainfall  (Ellrod,  1994).  Arking  and  Childs  (1985)  used  toe  visible,  shortwave  IR,  and  thermal  IR 
bands  on  toe  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  to  determine  cloud  fiaction,  cloud- 
top  temperature,  optical  thickness,  and  cloud  microphysical  properties.  Some  researchers  have  devised 
schemes  to  remove  toe  emitted  portion  of  toe  radiation  in  tois  channel  by  digitally  subtracting  toe  emitted 
radiance  of  anotoer  channel,  assumed  to  be  nearly  toe  same  in  toe  thermal  infrared,  in  order  to  isolate  toe 
solar  component  This  process  has  been  used  to  support  toe  study  of  global  feedback  mechanisms  caused 
by  clouds  in  a  global  change  scenario.  Other  researchers  (Allen  et  al.,  1990;  Lee  et  al.,  1996;  Kleespies, 
1995)  using  polar-orbiting  and  geostationary  satellite  instruments  have  used  tois  process  to  quantitatively 
investigate  continental  and  marine  cloud  properties  during  toe  daytime.  More  recently,  several  qualitative 
applications  of  tois  channel  have  appeared.  They  can  identify  low  clouds  over  snow  cover  and  identify 
toin  stratus  clouds.  The  unices  have  been  used  to  identify  shiptrack  effects  on  marine  stratus  caused  by 
changes  in  reflectivity  due  to  cloud  particle  size  changes  fixim  ship  exhaust.  It  has  also  been  suggested  that 
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the  SIR  band  can  be  applied  wMi  oflier  channels  to  detect  large  sqiercooled  water  droplet  regions  that 
present  icing  hazards  to  aircraft  (Menzel  and  Purdom,  1994;  Ellrod,  1996a). 


4.2^.4  Water  Vapor  Band 

Another  imager  channel  coming  into  more  frequent  usage,  known  as  the  water  vapor  channel, 
takes  advantage  of  upper  atmospheric  water  vapor  absorption  in  the  spectral  band  between  5.4  and  7.1  pm, 
and  most  commonly  centered  at  6.7  pm  in  the  IR  regime.  At  this  wavelength,  most  of  the  sensed  terrestrial 
energy  comes  from  the  middle  to  top  of  die  troposphere  roughly  between  300mb  and  600mb  (Kidder  and 
Yonder  Haar,  1995).  This  band  is  used  to  infer  middle  to  upper  tropospheric  water  vapor  or  humidities. 
Since  this  is  not  a  water  vapor  window  channel,  low  level  features  can  not  be  resolved  unless  the  sensed 
column  is  eictremely  dry  where  diere  would  be  a  lack  of  attenuation  of  surface  radiance  by  water  vapor.  In 
this  case,  the  water  vapor  image  would  appear  dark  as  is  the  convention  for  high  radiances  in  infrared 
imagery.  For  abundant  atmospheric  moisture  or  diick  high  or  middle  clouds,  the  imagery  appears  white 
owing  from  die  lack  of  IR  radiance  being  sensed.  While  diis  imagery  can  detect  clouds  that  are  deep  in 
moisture,  it  is  most  commonly  used  to  monitor  the  water  vapor  movement  for  cloud-free  regions  through  a 
satellite  image  looping  procedure.  Sinking  motions  witiiin  die  troposphere  can  be  implied  due  to  adiabatic 
warming  and  drying  processes  resulting  in  daric  regions  on  the  imagery.  The  opposite  is  true  for  rising 
motions  widi  adiaharif,  cooling  and  condensation  processes  resulting  in  whiter  areas  (looking  much  like 
clouds). 

4.2,2.5  Microwave  Region 

The  passive  imaging  of  die  microwave  portion  of  the  Earth's  EM  spectrum  has  been  a  rich  source 
of  data  on  precipitation,  atmospheric  and  sur&ce  temperatures,  ocean  surface  characteristics,  cloud  water, 
water  droplet  phase,  and  soil  moisture.  In  addition,  the  microwave  instruments  have  the  capability  to 
measure  the  polarization  of  the  sensed  radiation  from  die  land,  ocean,  and  clouds.  This  also  assists  in  the 
detection  of  the  above  phenomena.  While  not  the  purpose  of  this  study,  much  research  has  been  devoted 
towards  further  understanding  and  utilizing  the  breaddi  of  information  from  microwave  imagery.  Lawyer 
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(1995)  developed  a  satellite  derived  icing  index  based  mainly  on  passive  microwave  imaging  data  from 
polar  orbiting  satellites.  His  mediods  paralleled  those  of  Jones  and  Vonder  Haar  (1990)  and  helped  to 
retrieve  cloud  liquid  water  estimates  over  continental  areas.  Their  research  has  shown  that  with  proper 
knowledge  of  die  land  surfece  properties  (specifically  ground  moisture  and  surface  emissivity),  die 
satellite-retrieved  estimates  of  cloud  liquid  water  can  give  a  fiuriy  good  idea  of  the  cloud  regions’  icing 
potential  for  non-precipitating  clouds.  When  used  in  conjunction  with  other  remote  sensing  methods, 
passive  microwave  imaging  can  be  a  powerful  tool  for  the  researcher  or  forecaster. 

43  Analysis  and  Display  Methods  for  Satellite  Imagery 

Several  satellite  imagery  display  methods  diat  have  been  commonly  used  in  the  past  will  next  be 
Also,  some  new  display  and  analysis  techniques  involving  multispectral  combinations  of 
satellite  image  channels  will  be  described  prior  to  dieir  use  in  this  study. 

4.3.1  Black  and  White  /  Grayscale  /  Color  Enhancements 

The  typical  channels  of  early  GOES  imagery  were  primarily  displayed  individually  as  a  black  and 
white  (BAV)  linear  scale  image  widi  intermediate  gray  shades  denoting  varying  values  of  radiances  or 
albedos,  infrared  chaimel  images  were  usually  inverted  for  display  in  order  that  die  coldest  cloud  tops 
would  match  that  of  die  brightest  cloud  tops  in  die  visible  channel  imagery.  Image  processors  displayed  up 
to  256  (8-bit)  gray  shades  on  dieir  monitors,  ranging  linearly  from  black  (zero)  to  white  (255).  Although  a 
processor  was  capable  of  producing  256  tones  of  gray,  only  a  limited  range  of  them  can  be  discerned  by 
the  image  interpreter’s  eye.  Thus,  a  wide  assortment  of  gray  shading  enhancement  scales  to  produce  a 
contrast  for  the  observer  were  developed  to  take  advantage  of  the  diermal  characteristics  of  the  LIR 
imagery.  This  led  to  better  ways  to  assess  cloud  hei^t  and  type  dian  for  unenhanced  linear  scale  LIR 
images.  In  addition,  techniques  of  false  color  enhancement  (creating  color  shadings  unlike  that  which  is 
seen  by  the  human  eye)  have  greatly  improved  die  prospects  for  visual  interpretation  of  meteorological  and 
land  surfece  phenomenon  in  satellite  imagery.  It  has  been  found  that  such  colorized  imagery  is  easier  to 
interpret  dian  the  more  conventional  single-channel  BAV  imagery. 
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43  J  MuUispectnil  RGB  Color  Enhancement  Method 

Wiflj  better  imageiy  display  tedmology  and  improvements  of  die  diannels  on  more  recent  polar 
orbiting  and  GOES  satellites,  methods  incorporating  combinations  of  satellite  channels  (known  as 
multispectral  combinations)  have  come  into  greater  usage.  MuWspectral  analysis  involves  the  use  of 
different  spectral  channels  togedier  in  order  to  highlight  a  desired  parameter.  This  technique 

enables  a  fast  application  of  all  avaUable  channels,  and  avoids  the  visual  inspection  of  a  number  of  separate 
B/W  or  color  enhanced  images.  Some  atmospheric  phenomena  can  be  more  easily  identified  by  using 
satellite  Hata  fixim  more  than  one  channel  in  a  multispectral  image,  such  as  fog  or  thin  cirrus  over  low 
clouds.  Color  composite  imagery  is  not  a  new  concept  and  die  use  of  different  multispectral  techniques  to 
detect  atmospheric  or  surface  phenomenon  has  been  commonly  used  in  die  research  community  in  the  past 
two  decades.  Studies  1^  Colwell  (1983),  d'Entremont  and  ThomasOT  (1987),  Kleespies  (1989),  and  more 
recently  BeUec  and  Le  Gleau  (1992)  have  all  explored  ways  of  creating  new  multispectral  color  composite 
images  for  display  and  iqiplying  dieir  use  to  meteorological  features.  Andreadis  et  al.  (1995)  also  have 
gone  on  to  investigate  odier  means  for  displaying  a  color  composited  image  using  the  three  components  of 
color  as  perceived  by  the  human  eye,  namely  the  hue,  saturation,  and  intensity. 

One  metiiod  which  will  be  explored  in  diis  study  involves  simultaneously  combining  three  B/W 
satellite  images  into  a  color  enhanced  product  This  process  involves  the  assignment  of  each  of  the  three 
channels  each  to  one  of  the  color  "guns"  (red  (R),  green  (G),  and  blue  (B))  used  in  typical  television  image 
production.  The  resulting  composite  image  is  then  referred  to  as  being  RGB  enhanced.  For  each  color 
image  plane,  there  exists  a  range  of  resulting  colors  tiiat  tiie  color  guns  virill  produce  varying  with  tiie 
contribution  from  tiie  channel.  For  example,  tiie  red  gun  ranges  from  dark  red  to  pink  to  white  depending 
on  the  display  scale  for  the  image  pixels  according  to  each  channel's  display  convention  (white 
representing  hot  brightness  temperatures  for  uninverted  IR  or  cold  for  inverted  IR).  If  for  any  given  image 
pixel,  tiie  measured  brightness  temperatures  are  tiie  same,  then  the  contributions  to  the  color  guns  will  be 
tiie  same  (d’Entremont  and  Thomason,  1987).  When  any  tiiree  channels  are  combined,  equal  contributions 
of  red,  green,  and  blue  wUl  yield  a  shade  of  gray  in  tiie  color  image.  But  more  commonly,  the  varying 
contributions  of  each  channel  create  a  rainbow  of  colors  varying  in  hue  with  the  resulting  colors  a 
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<»ftmhiniitiftn  of  Ae  tcd,  green,  and  blue  guns.  Tbe  intermediate  colors  for  die  RGB  image  are  cyan  (green- 
blue),  yellow  (red-green),  and  magenta  or  purple  (red-blue).  To  quantify  the  RGB  color  combination  of 
the  diree  fiiannak  may  be  easier  when  all  diree  channels  are  in  die  IR  regime,  since  it  would  be  more  likely 
diat  the  brightness  temperatures  would  conelate  for  similar  RGB  contributions.  However,  the  resulting 
image  would  not  take  advantage  of  the  different  spectral  characteristics  for  visible,  SIR,  and  LIR  regimes. 
So,  in  the  case  of  RGB  combinations  of  different  spectral  channels,  the  inteipretation  is  best  left  to  a 
qualitative  understanding  of  the  phenomena  highhgbted  through  its  use. 


4  J J  Multispectral  Channel  Differencing  (AT,)  Method 

SateUite  channel  differencing  methods  have  been  around  for  quite  some  time  and  have  been 
frequently  used  in  cloud  classification  schemes.  Channel  differencing  commonly  involves  the  manual  or 
automated  digital  subtraction  of  brightness  temperatures  (T,)  for  each  image  pixel  of  two  different  satellite 
imagery  channeb  in  die  IR  regime  creating  new  composite  imagery  with  a  new,  resultant  brightness  value 
(ATb).  a  study  by  Shenk  and  Curran  (1973)  used  tins  method  to  detect  cirrus  clouds  and  infer  their 
heighte.  Crosiar  et  al.  (1989)  abo  used  this  technique  to  automate  cloud  detection  over  oceans  during  the 
night  Allen  et  al.  (1990)  showed  that  tins  method  can  help  to  distinguish  clouds  fi^m  underlying  surfrce 
snow  cover  or  ice  cover.  Dilb  et  al.  (1996a  and  1996b)  used  the  channel  differencing  method  to  infer  and 
detect  cloud  phase,  low  clouds,  fog,  and  surface  features.  Again,  these  combinations  of  the  various  IR 
satellite  channeb  provide  additional  valuable  information  for  detecting,  cloud  and  near-surface  features 
with  respect  to  using  each  IR  channel  separately. 

For  thb  study,  the  chaimel  differencing  metiiod  follows  that  of  Ellrod  (1996a)  (described 
previously  in  Chapter  2)  and  was  limited  to  GOES-8  channeb  2  and  4.  These  channeb  were  chosen  in 
order  to  utiiiyft  some  of  their  unique  spectral  characteristics  during  the  daytime.  These  spectral  properties 
have  been  previously  described.  To  maximize  tiie  detected  reflective  component  of  tire  CH2  radiance,  flie 
sftfffllitft  image  times  were  selected  during  a  daytime  period  wWi  a  high  solar  zenttii  angle  and  satellite 
viewing  angle  in  the  vicinity  of  the  region.  Then,  to  simplify  the  interpretation  of  the  daytime  3.9nm 
(CH2)  image,  the  emitted  radiation  (or  tiiermal)  component  was  subtracted  fiom  it  by  using  the  10.7 
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micron  (CH4)  radiance  creating  Ae  desired  AT,  fw  use  in  Ae  study.  Tlius,  dominant  liquid  water  cloud 
features  should  be  evident  fiom  die  magnitude  of  die  reflected  solar  component  with  the  greatest  ATb 
representing  smaller  liquid  droplets  and  moderate  AT,  indicative  of  the  more  hazardous  larger  cloud 
droplets  at  cloud  top.  The  idea  is  to  use  die  unique  daytime  spectral  jMOperties  and  radiance  values  of  diese 
two  re.  channels  and  dieir  variation  with  cloud  microphysical  properties  togedier  to  combine  into  a  channel 
differenced  product  Then,  after  comparison  to  the  sensed  albedo  fiom  the  visible  channel  1,  an  indication 
of  cloud  microphysics  (cloud  droplet  size  and  phase)  for  potential  icing  conditions  fiom  interrogation  of 
cloud-top  regions  should  be  the  result 

4.4  GOES-8  Imager 
4.4.1  GOES  Background 

In  die  early  1970's,  NASA  research  and  development  began  die  Geostationary  Operational 
Environmental  (GOES)  program  widiin  the  National  Oceanic  and  Atmospheric  Administration 

(NOAA).  Several  spin  stabilized  satellites  were  buUt  and  in  October  1975  GOES-1  was  launched,  which 
introduced  a  new  era  of  satellite  service.  Overseas  the  European  communities  have  maintained  a  series  of 
geostationary  satellites  (Meteosat),  as  have  the  Jqianese  widi  die  GMS.  A  good  summary  of  the  European 
Meteosat  satellite's  capabilities  is  contained  in  Mason  and  Schmetz  (1992).  Unlike  polar-orbiting  satellites 
which  travel  around  the  Earth  in  a  low  altitude  orbit  several  times  a  day,  these  geostationary  satellites 
remain  nearly  fixed  at  a  point  about  35.8  km  above  die  Equator.  Imagery  for  these  early  GOES  satellites 
were  provided  by  the  Visible  and  Infiared  Spin  Scan  Radiometer  (VISSR).  GOES  significantly  advanced 
our  ability  to  observe  weather  systems  day  and  night  by  providing  regular  interval  visible  and  infirared 
imagery  of  the  Earth's  surfece  and  atmospheric  moisture,  thus  becoming  a  critical  part  of  National  Weather 
Service  (NWS)  operations.  Since  GOES-4  was  launched  in  1980  through  die  recently  retired  GOES-7, 
several  more  spectral  bands  had  been  added  and  die  introduction  of  die  VISSR  Atmospheric  Sounder 
(VAS)  enabled  active  sensing  of  atmospheric  temperature  and  moisture  as  well  as  imaging  capabilities. 
Although  these  new  instruments  represented  improvements  in  satellite  sensing  skill,  some  compromises 
were  made  (Menzel  and  Purdom,  1994).  First,  imaging  and  sounding  could  not  be  done  at  the  same  time. 
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Secondly,  past  GOES  satellite  wwnging  radiometers  would  ^in  continuousfy  (100  ipm)  while  perfonning 
its  sensor  TTie  spinning  GOES-VAS  viewed  the  Earth  roughly  5%  of  the  time,  so  it  was  not 

possible  to  attain  signal-to-noise  ratios  needed  for  either  high-quality  soundings  or  greater  spatial 
resolution  data.  TTiis  left  the  sounding  data  to  be  taken  only  in  special  experiments  and  research  during 
dwell  sounding  mode,  when  no  imaging  was  to  be  undertaken. 

Understanding  the  need  for  a  better  imaging  and  sounding  platform,  NOAA  began  in  1985  the 
development  of  the  next  generation  of  sateUites,  GOES  I-H  each  designed  for  a  five-year  lifetime. 
Finally,  in  April  1994,  GOES-8  (formerly  GOES-I  in  development)  was  successfully  launched  and  now 
resides  at  75“W  (Figure  4.7).  This  was  followed  nearly  one  year  later  by  GOES-9,  which  will  be 

moved  to  its  permanent  location  near  135“W  after  its  checkout  period.  The  most  significant  change  ftom 
previous  GOES  satellites  results  from  its  stabilization  technique.  GOES-8  is  three-axis  stabilized,  which 
enables  the  instruments  to  always  be  pointed  at  die  Ear*.  This  system  allows  better  multispectral  imaging 
accuracy  with  improved  spatial  resolution,  the  sounder  and  imager  (now  separate  instruments)  to  collect 
data  at  the  same  time,  more  stable  calibration,  and  mote  precise  image  flaming  registration.  Figure  4.8 
shows  a  sketch  of  the  GOES-8  imager.  The  new  earth-oriented  geostationary  spacecraft  design  required 
new  scanning,  navigation,  cahbration,  and  thermal  control  systems  (Menzel  and  Purdom,  1994).  For  more 
efficient  imaging  wMi  near-continuous  viewing  of  the  earth,  the  scaiming  system  moves  in  boustrophedon 
fashion  (i.c.  moving  (slewing)  back  and  forfli  west  to  east  then  east  to  west).  Also,  more  precise  pointing 
techniques  using  a  star  sensing  system  as  a  fixed  reference  allows  imager  and  sounder  pfacels  to  be  located 
within  4  km  at  nadir.  To  improve  calibration,  looks  to  space  and  reference  blackbodies  are  external  to  the 
instrument  unlike  the  previous  GOES-VAS.  Finally,  to  accommodate  the  large  daily  thermal  changes 
(varying  by  tens  of  degrees  Kelvin)  diat  occur  witii  staring  instruments,  better  monitoring  of  thermal 
gradients  and  focal  properties  has  been  introduced.  This  means  that  the  sensor  takes  no  longer  than  30 
miniitpg  between  looks  at  its  blackbody,  ensuring  that  the  calibration  equation  linear  coefficients  remain 
within  specified  noise  levels. 


68 


NORTH 


Figure  4. 7  Sketch  of  GOES-8  spacecraft  (from  Kidder  and  Yonder  Hoar  1 995) 
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Figure  4.8  Sketch  of  the  GOES~8  Imager  (from  Kidder  and  Yonder  Hoar  1 995) 

4.4.2  Imager  Characteristics  Compared  to  GOES-VAS 

The  GOES-8  imager  has  a  five  channel  multispectral  c^ability  consisting  of  a  visible  band 
(Channel  1:  0.52  -  0.72  pm),  three  infirared  window  bands  (Channel  2':  shortwave  IR  (SIR)— 3.78  -  4.03 
pm,  Channel  4:  diermal  IR  (LIR)-I02  -  1 1 2  pm,  Channel  5:  another  longwave  IR-1 1 .5  - 12.5  pm),  and 
a  water  v^)or  absorption  band  (Channel  3t  6.47  -  7.02  pm).  A  summary  of  the  GOES-8  imager's 
instrument  and  chatmel  characteristics  can  be  found  in  Tables  4.1a  and  4.1b  at  the  end  of  this  section. 
Table  4.2  describes  the  instrument  features  for  both  the  GOES-8  imager  and  sounder.  The  scanning  part  of 
the  radiometer  is  a  mirror  which  reflects  energy  into  a  telescope.  Both  visible  and  infiared  energy  are 
handled  by  Ae  same  optics,  however  different  detection  elements  are  used  for  different  spectral  bands. 
Visible  data  is  swept  out  8  lines  at  a  time  wiA  eight  silicon  photodiode  detectors  (Kidder  and  Yonder  Haar, 
1995).  The  IR  window  channels  are  swept  out  two  lines  at  a  time,  while  Ae  water  vapor  channel  is  swept 
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out  one  line  at  a  time.  Channels  3-5  use  a  metcury  cadmium  tellurium  (HgCdTe)  detector  and  Ae  SIR 
(Channel  2  uses  an  indium  antimonide  (InSb)  detector.  The  spectral  response  of  each  band  is  determined  by 
the  sensor  optics,  such  as  beam  splitters  or  filters,  and  by  Ae  detector  properties.  The  imager  provides  Ae 
visAle  at  nearly  1  km  ground  resolution  (comparable  to  GOES-VAS),  but  whh  10-bit  precision  and  a 
more  stable  linear  response  to  radiance.  The  detectors  sense  a  voltage  that  is  converted  into  a  digital 
brightness  count  value  (EUrod  and  Nelson,  1996b).  The  transmitted  data  stream  contains  calibration 
coefficients  for  users  vrishing  to  convert  Ae  Agital  count  values  to  raAances  (IR)  and  albedos  (VIS). 

A  comparison  of  Ae  characteristics  of  Ae  GOES-7  (VAS)  and  GOES-8  imagers  can  be  found  in 
Table  4.3a.  By  using  star  references  and  traditional  landmarks,  Ae  imagery  is  navigated  to  wiAin  2-4  km 
compared  to  3-10  km  wiA  VAS.  The  GOES-8  imager  provides  simultaneous  IR  imagery  in  four  channels 
(one  more  Aan  VAS).  The  ground  resolution  of  Ae  IR  channels  is  much  improved  over  Aat  possible  wiA 
VAS.  The  SIR  chaimel  has  roughly  four  times  Ae  resolutioit,  and  Ae  longwave  IR  channels  have  twice  Ae 
resolution,  while  Ae  water  v^ror  chatmel  has  nearly  Ae  same  resolution.  At  nadir,  Ae  ground  resolution  is 
4  km  for  Ae  IR  window  bands  and  8  km  for  Ae  water  vsqror  band.  The  radiometric  accuracy  and  onboard 
calibration  yield  brightoess  temperatures  wiAin  1  K  accuracy  and  0.3  K  relative  precision  at  most 
temperature  ranges  (Menzel  and  Purdom,  1994).  Additionally,  Ae  noise  levels  are  reduced  nearly  Aree 
times  over  Ae  GOES-VAS.  The  GOES-8  imager  also  slightly  reduced  Ae  amount  of  oversampling  of  its 
mfrared  channels,  because  of  its  better  signal-to-noise  ratios  Aan  Ae  VAS. 

For  real-time  weaAer  forecasting,  Ae  GOES-8  data  must  be  transmitted  nearly  instantaneously  to 
field  users.  A  new  ground  data  processing  system  for  Ae  next  generation  GOES  satellites  permits  efficient 
dissemination  of  imagery  data  and  products  to  users  (Menzel  and  Purdom,  1994).  This  is  done  using  a 
completely  redesigned  GOES  I-M  Variable  (GVAR)  format  for  direct-receive  users  (like  Ae  Cooperative 
Institute  for  Research  in  Ae  Atmosphere  -  CIRA/CSU)  and  existing  GOES-Tap  services  for  oAers.  The 
GVAR  format  lets  users  receive  real-time  calibrated  information  fittm  Ae  imager  and  Ae  sounder. 
Because  of  Ae  Afierence  in  sampling  fiequency  and  resolution  of  GOES  I-M  versus  Ae  previous  GOES- 
VAS,  direct  receive  imagery  fixtm  GOES-8  appears  stretched  in  Ae  east-west  direction  wiA  respect  to 


71 


GOES-7.  Hus  data  stream  is  also  being  archived  fw  future  use,  like  previous  GOES  data  starting  in  1978. 
User  access  to  the  GOES  archive  is  coordinated  through  die  NESDIS  SateUite  Data  Services  Division. 

Wifli  Ae  recent  successful  launchings  of  two  next-generation  satellites  (GOES-8/9)  and  dieir 
upgraded  sensor  packages,  researchers  are  rushing  to  apply  the  improved  remote  sensing  capabilities. 
Because  the  new  series  of  GOES  satelUtes  produces  imagery  images  as  often  as  every  15  minutes  widi 
much  better  resolution,  the  potential  for  using  these  products  operationally  has  greatly  increased.  This  is 
especially  true  for  die  SIR  Channel  2  (3.9  pm),  since  c^bilities  in  this  channel  on  past  satellites  have 
been  very  limited.  The  signal-to-noise  ratio  for  the  SIR  band  approaches  that  of  the  LIR  band  at  warm 
temperatures,  which  helps  in  die  quantitative  analysis  of  Eaidi  surfece  and  cloud-top  properties.  Better 
satellite  imagery  display  systems  with  looping  ciqiabilities  also  promise  to  extend  their  usefulness.  In 
addition,  the  scan  patterns  for  die  GOES-8  imager  is  much  more  flexible  than  those  on  the  GOES-VAS.  It 
is  even  possible  to  suspend  a  fuU-disk  scan  to  perform  a  rapid  scan  of  a  small  area,  and  then  resume  the 
larger  scan.  This  capability  lets  researchers  or  severe  weadier  forecasters  to  more  finely  interrogate  a  small 
geographic  region  widi  maximum  resolution  at  very  diort  time  intervals  (on  die  order  of  minutes 
depending  on  the  area  size). 

4.43  Spectral  Improvements  Over  GOES-7 

In  dieir  paper,  Menzel  and  Purdom  (1994)  did  a  thorough  comparison  of  the  spectral  packages  for 
the  imager  on  GOES-7  and  the  improved  imager  on  GOES-8.  Though  published  before  the  launch  of 
GOES-1  (8),  their  simulations  of  GOES-I  imagery  and  knowledge  of  the  different  channel  characteristics 
provide  a  good  source  expected  improvements  (listed  in  Table  4.3b).  These  channel  by  channel  advances 
are  summarized  below.  In  die  visible,  major  improvements  are  in  lO-bh  versus  6-bit  imagery  and 
increased  sampling  fiequency.  The  oversampling  should  allow  for  better  cloud-edge  detection,  while  die 
data  storage  improvement  will  provide  sixteen  times  more  available  (1024  versus  64)  brightness  levels  or 
gray-shades.  This  enables  a  wider  range  of  image  enhancements  that  could  upgrade  cloud-top  and  cloud- 
edge  feature  detection  for  use  in  determining  cloud-drift  winds,  better  quality  and  extended  use  of  visible 
imagery  in  low-light  situations,  potential  detection  of  haze  and  pollution,  and  improved  daytime  cloud 
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height  measurements  using  both  stereo  and  cloud  shadow  techniques.  For  the  SIR  window  channel,  major 
improvements  in  resolution  (2x4  km  vs.  4x16  km)  and  better  radiometric  sensitivity  at  warmer 
temperatures  wUl  greatly  increase  the  ability  to  detect  fog  at  night,  locate  water  clouds  over  snow  in  the 
daytime,  delineate  between  SLW  and  ice  clouds  during  die  daytime,  and  locate  thermally  hot  surface  areas. 
It  should  also  help  to  identify  die  hurricane  eye  when  it  is  obscured  by  this  cirrus  clouds  and  better 
determine  sea  surface  temperatures.  In  the  water  vapor  channel,  the  spatial  resolution  is  enhanced  by  twice 
that  of  the  VAS  and  signal-to-noise  is  improved  by  three  times.  This  should  provide  better  identification  of 
middle  and  upper  level  moisture  patterns  and  synoptic  weather  features  when  viewed  in  a  looping 
sequence.  In  die  LIR  channel,  a  near  fourfold  increase  in  spatial  resolution  will  help  to  provide  more 
refined  cloud-top  features,  severe  storm  identification,  and  heavy  rainfaU  estimations.  The  split  window 
(11.5  -  12.5  pm)  channel  will  experience  an  eight-fold  increase  in  resolution  and  covers  a  wider  spectral 
widdi  than  VAS.  Its  resulting  signal-to-noise  refinement,  same  as  the  LIR  channel,  should  give  more 
accurate  low-level  moisture  indications  and  enhanced  sea  sur&ce  temperature  retrievals. 
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Table  4.  la  GOES-8  imager  instruments '  characteristics  (from  Kidder  and  Yonder  Hoar  1 995) 


Farameter 

Value 

Optics 

Cassegrainian  telescope 

Earth  location  accuracy 

31.1-cm-diameter  primary  mirror 
Focal  lengdi  381.2  cm 
±30  /trad  (rms) 

Qiannel  coregistration 

±28  ^rad 

Scanning  rate 

Full  disk 

25  min 

3000  X  3000-km  area 

3.1  min 

500  X  500-km  area 

20  s 

Channels 

Five  (0.65,  3.9,  6.75,  10.7,  12.0  fim) 

Digitization 

10  bits 

Date  rate 

2.6208  Mbits  s’* 

In-orbit  calibration 

IR  (channels  2—5) 

Space  and  internal  290-K  blackbody 

Visible  (channel  1) 

None 

Size 

Sensor  module 

115  X  80  X  75  cm 

Electronics  module 

67  X  43  X  19  cm 

Power  supply 

29  X  20  X  16  cm 

Mass  (3-module  total) 

120  kg 

Power  consumption 

119  W,  daily  average 

Table  4.  lb  GOES-8  imager  channels'  characteristics  (from  Kidder  and  Yonder  Hoar  1 995) 


Channel 

1 

2 

3 

4 

5 

Wavelengdi 

(Atm) 

0.55-0.75 

3.8-4.0 

6.5-7.0 

10.2-11.2 

11.5-12.5 

Detector  type 

Silicon  pho¬ 
todiode 

InSb 

HgCdTe 

HgCdTe 

HgCdTe 

Noise  or 

NEAT 

NEAL 
(mW  m"^ 
sr”*  cm) 

±8  counts 

0.15  K  at  300  K 
3.5  K  at  230  K 
0.0058 

0.3  K  at  230  K 

0.044 

0.2  K  at  300  K  0.2  K  at  300  K 
0.4  K  at  230  K  0.4  K  at  230  K 
0.31  0.35 

IFOV  (Atrad) 

28 

112 

224 

112 

112 

Ground 

resolution 

(km) 

1.0 

4.0 

8.0 

4.0 

4.0 

Sampled 
resolution 
at  subpoint 
(E-W  X  N-S, 
km) 

0.57  X  1.0 

2.3  X  4.0 

2.3  X  8.0. 

2.3  X  4.0 

2.3  X  4.0 
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Table  4.2  GOES-8  imager  and  sounder  instrument  features  (from  Memel  and  Purdom  1 994) 


Feature 

Imager 

Sounder 

Optical  aperture 

^1.1  cm 

31.1  cm 

Type  optics 

Cassegrain 

Cassegrain 

Methods  of  scan 

Two  axes,  continuous 

Linear  E/W  64  /zrad  (2.3  km) 

Line  step  M/S  224  /irad  (8  km) 

Two  axes,  step  and  dwell 

E/W  280-/irad  steps 

N/S 1 120-//rad  steps 
(or2240-/irad  with  0.1- 
or  0.2-s  dwell) 

Spatial  resolution 

Visible  28  ;irad(1  km) 

IR  windows  112Airad(4km) 

H2O  band  224 /irad  (8  km) 

242/tfad(10km) 

Sampling 

Visible  1.75/IGFOV* 

IR  windows  1 .75/IGFOV 
Hpband3.5/IGFOV 

Four  IGFOVs  sampled 
at  the  same  time 

Sampling  rate 

20®  r' 

1 83.3  /isec  per  pixel  (IR) 

45.8  /isec  per  pixel  (vis) 

40  soundings  s*' 

0.1 , 0.2,  or  0.4  s  per  sample 

Spectral  band 
coregistration 

±28 //rad 

Within  22 //rad  of  IR 

10.7-/tfn  window 

Data  output 

10-bit  quantization 

13-bitquantization 

Data  rate 

2.6208  Mbs-' 

40kbs-' 

Time  between 
space  looks 

2.2s 

(nominally  for  large  frame) 

9.2  or  36.6  8 

(nominally  for  small  frame) 

2  min 

Time  between 

blackbody 

calibrations 

10-30  min 

20  min 

'Instantaneous  geometric  field  of  view 
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Table  4.3a  GOES-7  and  GOES-8  imager  characteristics  (from  Menzel  and Purdom  1994) 


Wavelength 

Unn) 

IGFOV(km) 

E/WxN/S 

8SR(lan) 

E/WxN/S 

Noise 

GOES-7 

055-0.75 

0.75x0.86 

0.75x0.86 

6-bit  data  ±  2  counts  3  (T 

3.84-4.06 

13.8x13.8 

3.0x13.8 

0.25  K  ®  300  K,  6.00  K  @  230  K 

6.40-7.08 

13.8x13.8 

3.0x13.8 

1.00  K  ®  230  K 

10.4-12.1 

6.9  X  6.9 

3.0X6.9 

0.10  K  0  300  K,  0.20  K  0  230  K 

12.5-12.8 

13.8x13.8ir 

3.0x13.8 

0.40  K  O  300  K,  0.80  K  ®  230  K 

GOES-1 

0.52-0.72 

I.Oxl.0 

0.57x1.0 

1 0-bit  data  ±  8  counts  3  (T 

3.7&-4,03 

4.0  X  4.0 

2.3  X  4.0 

0.15  K  O  300  K.  3.50  K  0  230  K 

6.47-7.02 

8.0  X  8.0 

2.3X6.0 

0.30K0230K 

10.2-11.2 

4.0  X  4.0 

2.3  X  4.0 

0.20  K  9  300  K.  0.40  K  0  230  K 

115-12.5 

4.0  X  4.0 

2.3  X  4.0 

0.20  K  0  300  K.  0.40  K  ®  230  K 

Table  4.3b  Expected  improvements  in  GOES-8  imager  products  (from  Menzel  and  Purdom  1 994) 


^ 

•  More  details  in  imagery  (4-km  IR  resolution,  oversampted 
visible,  better  signal  to  noise,  higher  bit  depth). 

•  Improved  composite  imagery  (five  spectral  bands). 

•  Low-light  visible  imagery  (1 0-bit  visible  data) . 

•  Better  synchronization  with  other  observations  (separate 
imager). 

•  Better  cloud-drift  winds  (4*km  resolution,  better  edge 
distinction). 

•  improved  water  vapor  motion  (winds  in  dear  regions). 

•  Enhanced  severe  storm  forecasting  (timely  rapid  imaging, 
derived  products  at  4  km,  10-bit  visible  data). 

•  Timely  fog  dete^on  at  night  (continuous  3.9-Mm  imaging). 
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Chapter  5  Icing  Case  Studies  Research 
5.1  Overview 

Tnuiitionally,  aircraft  icing  forecast  products  have  specified  threat  regions,  heights  and  bases,  and 
occasionally  icing  type.  These  forecasts  have  been  accomplished  using  a  variety  of  different  metiiods, 
several  of  which  were  discussed  in  Ch^ter  2.  They  also  exhibit  a  large  variability  in  forecast  accuracy. 
The  information  used  to  prepare  these  forecasts  has  mainly  included  visible  and  infiared  satellite  imagery, 
ouq)ut  fi-om  numerical  weather  prediction  models  and  cloud  advection  models,  precipitation  forecasts,  and 
pilot  reports.  This  study  exploits  the  individual  characteristics  of  three  GOES-8  imager  channels  (VIS, 
SIR,  and  LIR)  along  widi  varying  cloud-top  microphysical  properties  and  their  change  with  each  spectral 
band  to  develop  a  scheme  that  will  infer  a  potential  hazard  for  aircraft  icing.  Much  of  the  background  has 
been  laid  out  in  the  previous  chspters  describing  likely  fisctors  that  could  lead  to  a  technique  to  better 
isolate  the  aircraft  icing  environment  and  improve  horizontal  and  vertical  spatial  forecasting  of  this  region. 
In  consideration  of  the  operational  forecaster,  the  me&ods  have  been  purposely  kept  as  simple  as  possible 
in  order  to  fiicilitate  its  usage  in  their  environment 

To  conduct  a  study  of  this  type,  it  proved  necessary  to  retrieve  the  verification  data  for  the  study 
first.  This  was  only  logical,  since  the  analysis  of  satellite  imagery  would  be  moot  without  a  way  to 
determine  the  accuracy  of  a  predicted  icing  wivirorunent  or  icing  forecast  through  the  only  known 
quantification  of  aircraft  icing  on  a  large  scale  basis,  namely  pilot  reports  (PIREPs).  Also,  to  use  the 
PIREPs  effectively,  it  was  essential  to  understand  the  proper  statistical  or  qualitative  use  of  PIREPs  for  a 
research  study,  because  as  described  in  Ch^ter  2,  diey  are  not  tire  most  reliable  or  easiest  way  to  confirm  a 
scientific  study  (due  to  their  subjectivity).  The  bulk  of  the  data  collection  process  time  involved  trying  to 
find  a  source  of  PIREPs  in  a  form  which  could  be  readily  sorted  by  type,  intensity,  altitude,  and  latitude  / 
longitude,  and  conveniently  displayed  or  plotted.  After  trying  unsuccessfully  with  the  Aviation  Weather 
Center  (AWC;  formerly  the  National  Aviation  Weather  Advisory  Unit  (NAWAU))  in  Kansas  City  and 
through  contacts  in  die  Air  Force  in  die  Air  Weadier  Service  to  get  the  needed  PIREP  data,  a  search  of  die 
Internet  provided  an  excellent  source  of  this  PIREP  data,  and  accompanying  decoding  and  plotting 
software  (FORTRAN,  C,  and  NCAR  Graphics  computer  code)  from  die  National  Center  for  Atmospheric 
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Reseaieh  (NCAR)  Research  and  AppUcations  Program  (RAP).  Uhanks  mainfy  to  Greg  Thompson,  CSV 
grad.  Dr.  Barbara  Brown,  and  others  at  NCAR/RAP:  CSV  grad  Ben  Bernstein,  Dr.  Marcia  Politovich] 

Once  die  desired  PIREPs  had  been  decoded  and  plotted,  the  optimal  case  study  periods  were 
determined  through  the  use  of  prevailing  synoptic  conditions  (significant  atmospheric  moisture  and 
sufficiently  cold  atmospheric  temperatures)  and  noting  periods  of  a  large  number  of  PIREPs  collections. 
These  fectors  helped  to  narrow  the  likely  study  periods  to  the  transitional  seasons  with  a  focus  on  late 
winter  and  early  springtime  in  particular.  With  a  narrowed  scope  for  the  case  studies,  my  attention  turned 
back  to  the  satellite  imagery  collection  and  finally  the  analysis.  Because  of  changes  occurring  with  the 
CIRA/CSU  GOES-8  imagery  collection  S3rstem,  a  few  data  archive  losses,  satellite  movement  from  its 
checkout  position,  and  image  display  format  changes  from  the  fall  of  1994  into  the  spring  of  1995,  the 
cases  selected  for  fliis  study  were  limited  to  tire  period  after  1  April  1995.  After  this  date,  the  available 
archive  of  raw  GOES-8  data  was  much  more  stable  and  convenient  to  retrieve  thanks  to  some  software 
developed  by  Duane  Whitcomb  of  CIRA/CSU.  The  SubSect  software  package  allows  the  user  to  retrieve 
near  real-time  GOES-8  imager  (and  sounder)  data  for  any  geogr^hic  area,  as  determined  by  the  user, 
scanned  by  the  satellite.  The  defined,  sectorized  images  can  then  be  viewed  using  another  piece  of 
software  written  in  IDL  (either  Imgbar36,  Imgbar41,  or  GVIEW  standing  for  GVAR  View).  Upon  finding 
that  the  imagery  was  available  for  flie  desired  case  shufy  dates,  die  GOES-8  imager  data  for 

channels  1,  2,  and  4  was  restored  fiom  the  archive  tapes.  Then,  die  analysis  began  using  techniques 
described  below. 

5.2  Data  Collection 

5.2.1  PIREPs 

5.2.1.1  Collection  and  Display 

All  of  the  PIREPs  came  from  die  NCAR  storage  t^s  and  were  collected  via  an  anonymous  ftp  to 
a  guest  user  directory  into  which  they  were  transferred.  They  came  in  an  encoded  format  that  required 
software  to  decode  and  further  software  manipulation  to  create  an  icing  gr^hic  plot  against  a  background 
m^  of  the  United  States  (as  in  Figure  5.1).  This  software  was  also  made  available  to  me  via  anonymous 
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^  to  NCAR.  The  massive  number  crunching  of  flie  PIREPs*  software  was  accomplished  with  much  help 
from  my  research  group’s  computer  system  expert,  Rick  Taft,  on  my  remote  link  from  Galileo  to  Tornado. 
All  raw,  decoded,  and  plotted  PIREPs  are  separated  into  hour  segments  diroughout  die  case  study  period 
and  stored  on  Tornado  (no  doubt  taking  up  much  space).  The  NCAR  Gr^hics  formatted  plots  are  also 
viewed  on  Tornado  using  the  "idt"  command.  Each  plot  can  be  tailored  to  plot  specific  intensities  by 
changing  the  configuration  file  or  by  selecting  coUection  periods  vaiying  in  length  of  time  for  the  plot 
using  multiple  PIREP  file  hours  to  build  a  larger  file. 

5.2.1.2  Selection  for  Verification 

Every  PIREP  within  a  3  hr  PIREP  file  period,  before  or  after,  die  analyzed  satellite  image’s  time 
were  used  in  the  verification  of  die  mediods  listed  below  applied  for  icing  conditions.  This  window  of  time 
was  chosen  so  that  an  adequate  number  of  PIREPs  could  be  obtained  to  verify  each  time  period. 
AdditionaUy,  the  6  hr  window  takes  into  account  the  reporting  timing  errors  that  are  inherent  in  the  PIREP 
database.  Studies  have  shown  diat  by  setting  a  certain  temporal  window  for  PIREP  collection  around  die 
desired  analysis  time,  the  statistical  benefits  of  the  study  verification  far  outweigh  the  possible  errors  that 
may  be  introduced  (Politovich  and  Olson,  1991;  Brown  et  al.,  1993). 

5.2.2  Case  Selection 
5.2.2.1  Casel  (10  Apr  95) 

Case  1  covered  an  area  approximately  2400  km  by  1200  km  over  the  Northern  Rockies,  Northern 
and  Central  Plains,  and  western  Great  Lakes  of  the  central  United  States.  This  area  (grid  region  shown  in 
Figure  52)  covered  portions  of  Montana,  Wyoming,  Colorado,  North  and  South  Dakota,  Nebraska, 
Kansas,  northern  Oklahoma,  Minnesota,  Iowa,  Missouri,  Wisconsin,  lUinois,  Michigan,  Indiana,  Kentucky, 
and  western  Ohio.  The  date  and  time  (10/1745Z  April  1995;  Julian  Day  100)  were  mainly  chosen,  because 
of  a  known  abundance  of  PIREPs,  on  which  to  base  the  study  verification,  were  collected  during  the  day  in 
the  vicinity  of  a  dynamic  spring-time  storm  system.  The  necessary  elements  for  icing  to  occur  on  a  large 
scale  were  present  as  abundant  warm.  Gulf  of  Mexico  moisture  streamed  northward  to  meet  the  cold, 
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continental  polar  air  mass  plunging  southward  out  of  Canada  (see  Figures  53  and  5.4).  The  air  masses 
converged  on  the  Central  Plains  creating  a  wide  area  of  cloud  coverage.  For  both  cases,  as  would  be 
expected,  the  convective  (cumuliform)  clouds  dominated  the  cloud  shield  along  and  ahead  of  the  frontal 
zones  (both  cold  and  warm  fronts),  and  ncariy  uniform  stratiform  cloud  decks  traded  the  surfece  low  to  the 

nordi  and  west 

At  10/12Z  Apr  95,  a  strong,  cold  upper-level  low  pressure  system  (see  Figure  5.5)  was  continumg 
to  deepen  over  northern  New  Mexico  and  starting  to  move  north-eastward  into  the  Central  Plains.  The 
return  flow  of  Gulf  of  Mexico  moisture  was  enhanced  by  the  strong  pressure  gradient  between  the 
deepening  surface  low  pressure  center  in  northern  Texas  and  a  region  of  high  pressure  parked  over  the 
southeast  U.S.  (see  Figure  5.6).  A  high  pressure  system  (continental  polar  air  mass)  over  Ontario  and  the 
northern  Great  Lakes  was  helping  to  feed  a  broad  area  of  cold  air  into  die  Central  Plains  and  Great  Lakes 
regions.  Precipitation  in  the  form  of  snow  showers  occurred  most  of  the  day  over  eastern  Colorado  and 
north-eastward  flirough  centra!  Minnesota,  while  rain  showers  and  thunderstorms  were  prevalent  from 
eastern  Texas  through  eastern  Kansas,  western  Missouri,  and  into  northern  Illinois  (see  Figure  5.7).  Two 
of  the  most  important  ingredients  for  aircraft  icing  conditions  are  the  upper-level  temperatures  and 
moisture.  These  elements  are  clearly  evident  from  the  upper-level  analyses  of  10/12Z  at  850mb  (roughly 
5000  ft  (050)  altitude  above  MSL),  700mb  (roughly  10000  ft  (100)  above  MSL),  and  500mb  (roughly 
18000  ft  (180)  above  MSL),  as  temperatures  for  the  region  of  interest  ranged  from  +15‘*C  to  -28‘’C  and 
dewpoint  depressions  were  low  (between  0  and  5)  particularly  for  the  lower  levels  (see  Figures  5.8  and 
5.9).  By  10/1 8Z,  near  die  icing  analysis  time,  much  of  the  Northern  Plains  had  surface  temperatures  of 
+3»C  or  less.  Thus,  the  probability  was  high  that  the  potential  icing  conditions  could  have  extended  very 
close  to  the  ground,  which  is  a  potentially  hazardous  situation  for  smaller  aircraft. 

5.2.23  Case  2  (11  Apr  95) 

Case  2  included  much  of  the  same  area  in  the  central  United  States  though  slightiy  larger  with  the 
addition  of  portions  of  northern  New  Mexico,  northern  Arkansas,  and  Tennessee.  It  occurs  one  day  after 
Case  1  on  11/1745Z  April  95  (Julian  Day  101).  The  synoptic  situation  is  much  the  same,  yet  shifted 
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Figure  5.3  Example  visible  GOES-8  image  from  10/1 815Z  Apr  95 


CONUS  -  CH2 
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Figure  5.4  Example  SIR  GOES-8  image  from  10/1815  Apr  95 


Figure  5. 7  Weather  radar  summary  for  I0/2135Z  Apr  95 
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Figure  5.9  Upper  air  weather  chart  (700mb)  for  10/I200Z  Apr  95 
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slightly  eastward.  By  1 1/12Z,  the  upper-level  low  pressure  center  (see  Figure  5.10)  over  northwest  Kansas 
had  nearly  cut-off  torn  the  steering  flow  to  its  north,  allowing  it  to  drift  much  more  slowly  and  maintain 
the  icing  conditions  for  an  extended  period  of  time.  It  turns  out  that  this  synoptic  system  didn't  fully 
release  its  grip  over  the  continental  U.S.  until  14  April.  At  the  surfiice  (see  Figure  5.1 1),  the  maturing 
fi-ontal  system  and  occluding  low  pressure  center  over  Iowa  was  helping  to  spread  a  mixed  bag  of 
precipitation  types  over  the  central  U.S.  SnowfeU  continued  to  accumulate  in  a  band  from  eastern 
Colorado  through  Minnesota,  while  showers  and  tiiunderstorms  continued  from  Wisconsin  and  Illinois 
soutiiward  to  the  Gulf  coast  (see  Figure  5.12).  The  extensive  coverage  of  precipitation  and  snowfell 
amounts  can  be  seen  with  this  storm  system  (see  Figures  5.13  -  5.16).  From  the  1 1/12Z  Apr  95  upper-level 
analyses  (see  Figures  5.17  and  5.18),  it  is  evident  tiiat  some  of  the  necessary  parameters  for  icing 
environments  are  still  present  over  tiie  region  of  study.  Air  temperatures  range  from  +15°C  to  -27®C  with 
most  of  the  dewpoint  depressions  remaining  low,  representing  a  high  relative  humidity  for  levels  below 
15000  ft  MSL. 

5.2.3  Satellite  Imagery 
5.23.1  Collection 

Raw  10-bit  GOES-8  data  for  the  study  period  is  archived  on  digital  (8-mm  exabyte)  tapes, 
recorded  and  stored  by  CIRA.  The  imagery  data  is  condensed  into  8-bits  from  the  10-bit  raw  data  on  the 
tapes  to  facilitate  the  display  constraints  of  the  monitor.  Case  study  dates  were  selected  and  the 
corresponding  tapes  were  used  to  restore  the  raw  data  into  my  directory  on  the  CSU  Atmospheric  Science 
Building's  CIGS  Room  Windows  NT  (named  IRIS)  operating  computer  system.  With  the  help  of  the 
SubSect  and  SeetDef  image  processing  software,  flje  selected  raw  imagery  data  was  reduced  to  specified 
regions.  Primarily  using  the  GVIEW  function  of  tiie  IDL  display  software  written  by  Kelly  Dean  of 
CIRA/CSU  tiie  imagtss  were  displayed,  enhanced,  and  edited  on  a  1024  X  1024  pixel  screen. 
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Figure  5.11  Surface  weather  chart  for  ll/1800ZApr  95 
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Figure  5.14  24  hr  precipitation  summary  for  12/I200Z  Apr  95 
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Figure  5.16  Observed  snow  caver  summary  for  12/I200Z  Apr  95 


Figure  5.17  Upper  air  weather  chart  (850mb)for  ll/1200ZApr  95 


Figure  5.18  Upper  air  weather  chart  (700mb)  for  10/1200Z  Apr  95 


S23:i  Display 

GVIEW  enables  the  user  to  interrogate  toe  displayed  imagery  by  toe  computer  mouse  movement 
of  an  on-screen  cursor  (+)  to  a  desired  point  on  toe  image.  A  mouse  click  on  toe  point  creates  a  "pixel 
dump"  to  a  small  window  on  toe  screen  which  displays  toe  valuable  information  needed  for  this  study. 
Each  pixel  dump  gives  toe  X-Y  monitor  image  pixel  coordinates,  pixel  count  value,  brightness  temperature 
(Tb)  in  degrees  Kelvin  (for  CH2,  4)  or  albedo  value  (for  CHI)  ranging  from  zero  to  one,  and  the  Earth 
latitude  and  longitude  of  the  selected  image  pixel.  Since  there  is  not  a  quick  and  easy  way  at  this  time  to 
dump  this  data  directly  into  a  data  ffle,  care  was  taken  to  write  down  each  pixel  dump  value  from  a 
predetermined  sector  grid  into  an  empty  spreadsheet  grid  for  later  manual  entering.  For  future  reference, 
one  should  be  able  to  write  some  IDL  code  that  would  do  this,  but  it  was  beyond  the  reach  of  my 
capabilities  and  tone  limits  at  present.  AdditionaUy,  there  diould  be  a  similar  way  to  display  a  grid  of  pixel 
dump  values  in  a  user  determined  sector  size  to  easily  view  all  pixel  values  at  once  for  a  region  on  toe 
image. 


5.2.3J  Selection 

For  Cases  1  and  2,  pixel  dump  values  (brightness  temperature  for  IR  and  albedo  for  visible)  for 
fBch  of  the  three  channels  were  selected  according  to  toe  grid  region  (shown  in  Figure  5.2)  with  each  value 
assumed  representative  of  toe  entire  sector  (roughly  10x20  pixels).  While  this  method  is  quite  tedious  for 
large  areas  and  fine  grids,  it  represents  toe  simplest  way  to  retrieve  toe  required  data  values.  The  Case  1 
and  Case  2  grids  were  determined  by  toe  large  cloud  field  as  seen  from  toe  CHI  image  and  a  rough  idea  of 
atmospheric  temperatures,  since  for  icing  one  needs  a  cloud  and  a  certain  temperature  regime.  The  fact 
that  a  grid  is  used  instead  of  toe  exact  image  pixel  associated  with  each  PIREP  causes  some  problems  over 
regions  where  toe  cloud  cover  is  not  uniform  in  levels  or  texture,  as  would  be  toe  case  for  embedded 
cumulus  within  a  stratiform  region,  edges  of  cloud  shields,  or  for  areas  of  widely  broken  /  scattered  cloud 
top  elements  at  different  levels.  This  would  lead  to  somewhat  erroneous  satellite-sensed  cloud-top  values 
and  thus,  incorrectly  assigned  grid  values. 
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5^4  Rawinsonde  (Sounding)  Data 

The  rawinsonde  or  upper-air  sounding  data  used  to  support  diis  study  was  collected  witii  the  help 
of  Steve  Finley  (CSU/Atmos  Science.  Dept.  Computer  System  Coordinator),  who  maintains  archived  upper 
air  data  in  GEMPAK  format  A  C-programming  language  script  was  run  to  recover  the  requu-ed  upper-an 
soundings  on  the  desired  dates  for  stations  across  die  norfli-central  United  States  with  complete  coverage  of 
the  desired  grid  region  for  case  1.  The  data  was  collected  and  stored,  and  viewed  (using  Ghostview)  on 
separate  PC’s  (Aurora  and  Horizon)  in  the  Atmospheric  Science  Dept.  Weatiier  Ub.  These  in-situ 
soundings  give  the  most  accurate  knowledge  of  the  vertical  profile  of  ambient  air  and  dewpoint 
temperatures,  and  winds.  Unfortunately,  die  wide  horizontal  spacing  of  these  stations  (-  400  km)  through 
this  area  does  not  lead  to  a  great  resolution  for  the  region.  Their  representativeness  of  the  atmosphere 
around  the  station  extends  ^proximately  150-200  km  radially  or  more  so  downstream  of  the  station. 

The  sounding  data  in  5x>me  cases  can  give  a  good  idea  of  the  air  mass  and  environment  in  which 
the  icing  events  occurred  and  icing  altitudes  can  be  related  to  temperatures  using  U.S.  Standard 
Atmosphere  profiles.  However,  they  should  be  used  wMi  some  caution  in  the  vicinity  of  fi-ontal  zones 
where  air  masses  change  rapidly  due  to  the  turbulent  nature  of  the  convection  which  usually  takes  place 
there.  For  comparison  purposes  to  PIREPs,  diey  would  only  be  useful  if  within  these  areas  of  die  nearest 
rawinsonde  location.  Since  die  soundings  take  place  only  routinely  twice  a  day  (12Z,  OOZ),  some  care 
must  be  used  in  interpolating  between  the  soundings  for  comparisons  that  take  place  in  the  middle  of  the 
day.  The  hard  copies  of  the  soundings  are  shown  in  a  standard  GEMPAK  Skew-T,  Log-P  format,  which 
also  calculates  some  atmospheric  stability  indices  for  the  sounding. 

5^.5  Upper  Air,  Surface,  and  Forecast  Weather  Charts 

Weather  charts  of  surfiuse  and  upper-air  features  were  collected  fixim  the  archive  stored  across 
from  the  Weather  Lab.  Paper  copies  were  made  of  die  surface  and  standard  upper-air  level  analyses 
(SSOmb,  700mb,  500mb,  SOOmb)  that  coincided  as  close  as  possible  to  the  desired  times.  These  charts  can 
help  to  horizontal  temperature  and  moisture  regimes  across  the  grid  region.  The  surface  chart  was 

also  used  to  confirm  die  CH4  temperature  retrievals  of  surface  targets  with  known  land  and  water  surface 
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temperatures.  Radar  summary  charts  were  coUected  to  confinn  die  presence  of  thunderstorm  activity, 
precipitation,  and  cloud  base  and  tops.  Some  NGM  forecast  model  charts  were  collected  to  see  if  they 
were  accurately  handling  Ae  storm  systems  development  The  Eta  model  forecasts  are  a  basis  for  some  of 
the  i<^ing  forecasting  algorithms  in  use  today  at  the  Aviation  Wea&er  Center. 

53  Icing  Study 

53.1  Methodology 

For  this  study,  several  assumptions  need  to  be  made  and  an  understanding  of  the  potential 
crudeness  of  the  verification  items  (PIREPs)  must  be  weighed  against  the  benefits  of  their  use. 
Additionally,  some  limiting  factors,  such  as  computer  software  compatibility  and  display  capabilities, 
specific  PIREP  and  satellite  data  collection,  satellite  imagery  interpretation,  and  computer  programming 
and  time  limitations  must  be  taken  into  consideration  during  a  review  of  the  research  process. 

5.3.1.1  Temperature  Regime 

An  in-situ  aircraft  icing  study  by  Pobanz  et  al.  (1994)  determined  that  cloud-top  temperature 
(CTT)  ranges  of -5  to  -16*0  were  found  to  characterize  SLW  at  cloud-tops  during  known  icing  encounters 
with  large  drop  regions  by  a  research  aircraft.  While  cloud-top  temperatures  deduced  fi'om  infrared 
satellite  imagery  may  contain  some  error,  CH4  IR  brighmess  temperatures  can  still  be  used  in  a  relative 
sense  to  create  boundaries  or  threshold  techniques  to  infer  icing  conditions  for  a  cloud  system  with  the  help 
of  PIREPs  comparisons  and  verification.  For  Cases  1  and  2,  the  CTT  thresholds  for  the  suspected  icing 
zone  (IZ)  were  set  between  +3»C  and  -25“C  (or  276K  to  248K).  These  limits  are  set  to  take  into 
consideration  the  dynamic  temperature  increase  of  the  airframe,  inversion  layers  where  cloud  tops  may 
exist  warmer  than  the  interior,  and  moderately  deep  or  growing  cumuliform  clouds  with  slightly  colder 
cloud  tops.  During  die  study,  it  became  evident  dial  many  icing  reports  occurred  within  a  colder  CTT 
zone,  and  it  was  dubbed  the  convective  zone  (CZ).  The  CZ  ranged  fr'om  -25®C  to  -63®C  (or  248K  to 
210K). 
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53.1  J  Supercooled  Liquid  Water  (SLW)  Locations 

Detection  of  supercooled  liquid  water  (SLW)  cloud  droplets  is  essential  in  dre  detection  of  icing 
environments.  Popa  Fotino  et  al.  (1986)  found  from  statistical  analysis  that  using  ground  based  remote 
sensing  they  could  correlate  75  percent  of  die  confhmed  icing  cases  widi  SLW  measurements 

alone.  Ideally,  diese  mediods  will  be  able  to  distinguish  regions  of  SLW  or  more  importantly  large  SLW 
droplets  from  satellite  sensing  of  their  cloud  tops.  AddWonaUy,  Pobanz  et  al.  (1994)  found  that  under 
proper  conditions  large  SLW  droplets  tended  to  develop  and  congregate  in  layers  near  the  tops  of  some 
stratiform  clouds. 

533  Multlspectral  Channel  Dlfierenclng  (ATb  )  Study 

53.2.1  Past  Studies 

Brightness  temperature  differencing  mediods  (using  ATj  =  3.9nm  (SIR)  minus  lO.Tpm  (LIR)) 
have  been  previously  considered  by  Liljas  (1986),  Ellrod  (1996a),  and  others  previously  mentioned.  Liljas, 
using  the  SIR  and  LIR  channels  from  a  polar-orbiting  satellite,  encountered  some  interference  or  noise 
problems  that  curtaUed  a  thorough  study  of  this  differencing  method  for  daytime  usage.  Ellrod  has  just 
recently  put  fortii  the  first  paper  trying  to  use  quantified  ATg  results  to  detect  aircraft  icing  regions.  His 
study  (more  detail  in  Chapter  2)  concentrated  on  large,  stratiform  cloud  areas  during  the  winter  and  spring 
seasons,  used  visible  channel  imagery  to  gain  additional  qualitative  information,  and  also  relied  on  PIREPs 
as  a  rough  verification  tool.  Ellrod  determined  tiiat  there  seems  to  be  some  skill  in  correlating  icing 
PIREPs  to  potential  icing  regions. 

5.3.23  Application  to  Cases  1  and  2 

The  first  metiiod  used  for  the  April  cases  involves  the  use  of  GOES-8  imager  charmels  1,2,  and  4 
(VIS,  SIR,  LIR).  It  follows  closely  to  the  method  suggested  by  Ellrod  (1996a),  but  with  some  slight 
changes.  Much  of  tiie  descrqrtion  of  the  channel  differencing  (ATg  )  technique  is  described  in  Chapter  4. 
Digital  CH2  and  CH4  IR  brighmess  temperatures  (Tb)  and  CHI  visible  albedo  values  were  collected  for 
each  sector  within  a  region  (shown  in  Figure  53)  of  vast  cloud  cover,  some  of  which  was  convective 
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(cumulifonn)  and  some  stratiform.  The  PIREPs  (125  for  Case  1;  164  for  Case  2)  taken  roughly  within  3 
hours  of  die  satellite  imagery  time  were  matched  to  each  respective  sector  on  die  grid  and  corresponding 
values  of  Tb  and  albedo  were  assigned  to  each  PIREPs  by  location.  These  values  were  inserted  into  a 
spreadsheet  for  graphical  analysis  (reference  FIG  in  Ch^  6-  CH2-4.xls).  The  LIR  channek*  brighmess 
temperatures  were  then  subtracted  from  those  of  the  SIR  channel  (CH2  -  CH4).  The  result  gives  a  measure 
of  die  extra  reflected  radiance  component  sensed  from  die  cloud  tops  by  CH2,  since  during  the  daytime, 
the  emitted  component  of  the  CH2  radiance  is  nearly  (assumed)  equal  to  diat  of  CH4  (of  which  emission  is 
the  only  significant  component  of  the  sensed  radiance).  This  k  important,  because  as  shown  in  Chapter  2, 
the  variance  in  cloud-top  reflectivity  k  largely  a  function  of  the  cloud  particle  size  and  phase  (liquid  or 
ice).  For  a  suspected  icing  hazard,  die  cloud  must  have  SLW  droplets  present  at  cloud  top  or  within  the 
cloud,  and  larger  SLW  droplets  contribute  to  a  greater  icing  hazard.  Using  diese  differenced  IR  Tb  values 
and  the  vkible  albedo  values  for  the  PIREPs  together  in  a  scatter  plot  (see  Figures  6.1,  6.2,  6.8,  and  6.9), 
some  interBsting  and  some  seemingly  obvious  findings  can  be  noted,  which  are  discussed  in  the  next 

chapter. 

Since  the  cloud-top  CH4  Tb  values  can  be  roughly  assumed  to  equal  the  cloud-top  temperature 
(CTT),  this  information  can  be  useful  to  estimate  the  height  of  the  cloud  top  to  compare  to  the  icing  PIREP 
altitude  and  whether  die  cloud  top  temperature  folk  into  a  known  icing  temperature  zone  (IZ)  (using  +3®C 
to  -25®C  or  276K  -  248K).  Using  the  rawinsonde  data  for  stations  closest  to  each  PIREP,  the  icing 
altitudes  can  also  be  compared  to  determine  the  ^proximate  temperature  for  the  levels  at  which  the  icing 
occurred  and  sounding  estimated  CTTs  can  be  used  to  verify  the  precision  of  the  CH4  CTT  retrievals. 
Additionally,  upper-air  and  surface  weadier  charts,  and  radar  summaries  are  used  to  further  confirm  the 
overall  synoptic  conditions. 

53:23  Strengths  and  Weaknesses  of  Method 

While  thk  particular  method  can  be  used  only  during  the  daytime  when  there  is  a  Tb  difference 
between  CH2  and  CH4,  the  majority  of  the  PIREPs  occur  during  the  daytime  when  conunercial  and 
business  air  traffic  k  at  a  maximum.  Thk  metiiod  has  several  strengths  over  other  icing  prediction 
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mediods.  Most  importantly,  when  the  digitally  differenced  image  can  be  automated  for  display,  it  can  give 
a  real-time  indication  of  icing  dveat  areas  to  a  trained  operational  forecaster  fiom  visual  inspection  over 
large  geognqrhical  areas.  The  greater  resolution,  multispectrai  imt^ety  available  fiom  die  GOES-8/9  series 
of  diow  potential  to  improve  the  assessment  of  icing  hazards  over  continental  as  well  as  over 

regions  where  little  in-situ  measurements  are  available  on  vdiich  to  base  model  forecasts.  Though 
this  technique  also  has  several  weaknesses.  Unlike  the  icing  algorithm  forecast  methods  mentioned  in 
Ch^ter  2,  diis  method  would  have  to  be  automatically  combined  or  overlaid  with  other  data  to  infer  cloud 
depdis  and  heights  which  could  then  be  used  to  make  an  icing  forecast  Additionally,  diis  method  falters 
when  trying  to  sense  icing  conditions  for  most  convective  or  very  cold  ice  cloud  tops. 


5  J  J  Multispectrai  RGB  Color  Enhancement  Study 

Separate  channels  with  false  color  and  RGB  enhanced  images  were  also  investigated.  The  RGB 
mediod  each  a  satellite  image  channel  to  each  of  die  colors  (Red,  Green,  Blue).  The  resultant 

image  is  a  rainbow  of  colors  wMi  varying  hues  resulting  fiom  die  assignment  of  each  channel  to  a  range  of 
hues  within  each  primary  color.  RGB  or  other  false  color  methods  have  not  been  widely  used  in  the  past 
and  only  recentiy  have  more  studies  concentrated  on  the  potential  scientific  ^plications  for  these  mediods 
with  hope  of  developing  cloud  or  land  feature  identification  techmques. 

With  Kelly  Dean's  he^,  a  program  called  Ticolor"  running  in  IDL  was  created  to  select  individual 
channel  GOES-8  raw  imagery  and  assign  them  in  any  order  to  RGB  and  display  the  image  on  the  terminal 
screen  using  die  GVIEW  function.  Hiis  method  was  not  designed  to  determine  icing  regions 
quantitatively,  so  the  pixel  dump  value  was  not  a  concern  in  this  case,  due  mainly  to  the  difficulty  in 
^Ptwmining  and  Understanding  die  resultant  combined  value's  significance.  The  IR  channel  (2  and  4) 
imagery  were  also  inverted  so  as  to  compare  fevorably  with  die  CHI  image  (i.e.  Clouds  would  appear  white 
or  light  colored  for  all  images).  Using  die  techniques  described  in  Ch^ter  4,  regarding  Red-Green-Blue 
(RGB)  imng*^  enhancement,  two  different  combinations  of  channels  seemed  to  hold  some  promise  for 
<i.»»Arring  potential  icing  environments.  They  included  GOES-8  charmels  1, 2,  and  4.  The  addition  of  die 
visible  channel  to  the  IR  channels  helps  to  give  the  combined  image  more  textural  features.  The  two  RGB 
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combinations,  which  will  be  described  in  further  detail  in  Chapter  6,  are  indicated  by  the  three  channel 
numbers  as  assigned  to  the  sequence  of  red,  green,  and  blue.  They  are  refened  to  as  241  (CH2-red,  CH4- 
green,  CHl-blue)  and  41 1  (CH4-red,  CHI-green  and  blue). 


Chapter  6  ~  Resnits  and  Condnsions  of  the  Idng  Studies 

NOTE:  All  Figures  and  Tables  are  at  Ihe  ^  of  this  duqiter  (ordered  as  Aty  qipear  in  die  text). 

<1.1  Discussion  of  Results 

6.1.1  AT,  Method  for  Case  1-10  April  1995 

From  die  10  Apr  95  (Julian  Day  100)  case,  die  clourfy  region  over  die  Nordiem  plains  through  die 
western  Great  Traces  region  has  been  analyzed  using  die  multispectral  channel  differencing  (AT,)  metiiod 
mentioned  in  the-previous  chapter.  For  the  125  PIREP  set,  the  icing  intensities  covered  the  fuU  range  from 
a  Trace  (1)  to  Severe  (7),  although  over  half  of  the  icing  PIREPs  (46  of  86)  were  of  significant  intensity 
(Lgt-Mdt  or  greater,  denoted  as  (4+)  icing  PIREPs)  (see  Tables  6.1a<).  Overall,  the  cloud  shield  within 
the  area  of  interest  was  mostly  amvective  in  nature  with  only  smaU  repons  of  stratiform  cloud  tops  evident 
over  die  western  Plains  states. 

From  the  case  1  scatter  plot  of  CH2-CH4  (AT,)  vs.  CHI  (Figures  6.1a-<l),  we  find  a  clustering 
(81%)  of  all  icing  (1-7)  PIREPs  within  arange  of  cloud-top  AT,  between  20-44®C  and  albedo  values  of  42- 
60  and  the  same  is  true  (85%)  with  the  greater  intensity  (4+)  icing  PIREPs.  Most  of  the  outlying  PIREPs 
outside  of  the  range  can  be  explained  by  questionable  grid  values  along  cloud  edges  due  to  the  relatively 
large  gradient  of  the  values  across  a  short  distance  in  die  vicinity  of  the  measurement  or  in  regions  of 
broken  clouds  which  let  some  surfiice  radiation  through  between  the  clouds.  The  oudying  icing  reports 
could  also  resuh  from  the  feet  that  die  clouds  are  transient  and  are  not  always  propo-ly  represented  properly 
(using  this  gridded  scheme)  over  die  entire  six  hour  collectirm  period  for  the  PIREPs.  In  this  case,  the  4+ 
intensity  reports  generally  corresponded  to  greater  AT,  and  albedo  values.  Since  higher  albedo  values 
would  indir^*^  a  greater  visible  reflective  component,  characteristic  of  water  droplets  versus  ice  crystals, 
and  greater  positive  AT,  between  CH2  and  CH4  would  tend  to  result  from  warmer  and  liquid  (smaller) 
droplet  cloud  tops  widi  a  greater  reflective  component,  diis  analysis  would  provide  a  good  indication  of 
and  agreement  for  potmtial  icing  conditions  within  these  particular  cloud  regions.  However,  when  diese 
AT,  and  albedo  tango  tiiresholds  are  applied  to  the  spreadsheet  grid  of  pixel  dump  values  (shown  in  Tables 
62a-d)  by  coloring  the  grid  sectors  which  are  wMiin  tiie  ranges,  the  result  is  only  a  gridded  reflection  of 
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the  combined  cumulifoim  and  stratifonn  cloud  shield.  Ibese  ranges  for  which  the  PIREPs  clustered  just 
proved  that  diey  occurred  in  clouds,  which  is  a  given,  wMiout  respect  to  eiflrer  type  of  cloud. 

When  the  cumuliform  cloud  tops  (convective  zone  (CZ)  -  CH4  brightness  temperatures  less  than 
25*C)  were  screened  out  leaving  only  the  stratiform  cloud  tops  (see  Figures  62a-b)  (icing  zone  (IZ) 
CH4  brightness  temperatures  from  +3"C  to  -25“C)  grid  sectors  (shown  in  Tables  6.2c-d  /  shaded  in  the 
color  pink),  a  better  indication  of  icing  conditions  resulted.  While  die  shaded  regions  in  Tables  6.2c-d  for 
the  CH2-CH4  (ATg)  grid  values  represented  stratiform  regions  better,  it  did  also  highlight  the  cloud-top 
temperature  (CTT)  regions  along  the  cloud  shield  edges  of  convective  or  cirrus  clouds  which  are  not 
expected  to  be  used  in  die  icing  detection  technique.  The  AT,  values  fiir  die  IZ  were  in  the  higher  part  of 
die  range  (30-44)  which  is  representative  of  a  greater  reflected  component  from  CH2  for  liquid  droplets. 
But,  again  these  range  of  values  were  not  characteristic  of  icing  regions  only.  There  seemed  to  be  nearly  as 
many  instances  of  these  values  in  die  unshaded  CZ  portion  of  die  spreadsheet,  leading  to  the  result  diat  die 
channel  difference  mediod  is  not  a  particularly  stand-alone  technique  for  icing  detection  in  all  type  of 
clouds. 

Highlighting  the  region  of  the  IZ  using  only  CH4  Tb  on  the  grid  (see  Figure  6.3)  gives  a  better 
visual  representation  of  die  stratiform  areas  against  die  visible  CHI  image  for  that  time.  The  cloud  types 
are  also  qualitatively  verified  fiem  the  individual  GOES-8  channel  imagery  (see  Figures  6.4a-d).  Then, 
transferring  die  stratiform  cloud  area  to  die  composite  PIREP  plot  for  die  6  hr  period  encompassing  the 
analysis  time  (15-21Z)  (see  Figure  6.5)  gives  an  idea  of  how  the  icing  PIREPs  compared  to  the  stratiform 
regions  alone.  It  is  obvious  fimm  inspection  of  die  icing  PIREP  and  cloud-top  CH4  temperature  (CTT)  data 
that  Hiiring  Case  1  most  of  die  icing  PIREPs  coincided  widi  very  cold  CTTs,  which  were  mainly  in  the 
vicinity  of  convective  cloud  plumes  and  thunderstorm  activity.  The  greater  intensity  icing  PIREPs  of 
magnitude  4+  were  mainly  in  a  convective,  CTT  temperature  zone  (CZ)  between  248K  and  205K  (-25  C  to 
-68'’C)  (see  Figure  62b).  The  PIREPs  within  die  stratiform  region  were  sparse,  possibly  due  to  a  lack  of 
flight  through  die  area,  but  die  majority  of  the  PIREPs  in  die  area  were  of  the  significant  4+  intensity 
(particularly  over  western  Nebraska  and  Kansas).  A  rawinsonde  analysis  fmm  North  Platte,  NE  (LBF) 
(Figure  6.6)  shows  an  atmosphere  that  would  be  conducive  to  icing  conditions  at  altitudes  between  die 
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surface  and  15000  ft  AGL.  The  same  holds  true  for  die  Dodge  C%,  KS  (DDC)  (Figure  6.7)  sounding 
which  would  icing  conditions  far  tiie  levels  between  tiie  surface  and  9000  ft  AGL.  A  nice 

clustering  of  a  few  moderate  intensity  (4+)  icing  PIREPs  were  reported  in  tiiese  areas  and  at  tiiese  altitudes, 
which  also  means  that  tiiis  area  might  be  a  good  candidate  for  this  satellite  analysis  technique. 

The  introduction  of  tiie  negative  (none  or  nuU)  icing  remits  (Figure  6.1a)  to  the  analysis  did 
nothing  more  but  confirm  that  the  majority  of  the  null  PIREPs  were  beneath  very  cold  clouds  or  over  clear 
areas  where  tire  satellite  was  sensing  surfttce  temperatures  only.  This  is  often  the  case  when  the  CH2-CH4 
Tb  difference  (ATb)  is  smaller,  because  most  land  surfaces  are  not  very  reflective,  so  that  the  radiance 
sensed  in  CH2  is  nearly  the  same  as  that  in  CH4.  Also,  many  of  the  null  PIREPs  occurred  within  the 
stratiform  areas,  however  they  were  for  altitudes  above  the  stratiform  cloud  deck. 

6.1.2  AT,  Method  for  Case  2  - 11  April  1995 

From  tiie  1 1  Apr  95  case,  a  similar  cloutty  region  to  Case  1  over  tiie  Northern  Plains  through  tiie 
western  Great  Lakes  region  has  been  analyzed  using  tiie  multispectral  channel  differencing  (AT,)  method 
mentioned  in  the  previous  chapter.  For  the  164  PIREP  set,  the  icing  intensities  covered  the  full  range  from 
a  Trace  (1)  to  Severe  (7),  and  only  40%  of  tiie  icing  PIREPs  (44  of  109)  were  of  significant  intensity  (Lgt- 
Mdt  or  greater,  denot«vt  as  (4+)  icing  PIREPs)  (see  Tables  6.3a*d).  Overall,  the  cloud  shield  within  the 
area  of  interest  was  much  more  stratiform  in  nature  (wMi  respect  to  Case  1)  resulting  in  almost  equal  cloud 
coverage  between  the  convective  and  stratiform  regions.  Areas  wifli  regions  of  stratiform  cloud  tops  were 
evident  over  the  western  and  central  Plains  states,  and  over  northern  Wisconsin. 

From  the  scatter  plot  of  CH2-CH4  (AT,)  vs.  CHI  (Figures  6.8a-d),  we  find  a  clustering  (87%)  of 
all  icing  (1-7)  PIREPs  witiiin  a  range  of  AT,  between  10-44‘‘C  (or  somewhat  broader  (darker)  than  Case  1) 
and  albedo  values  of  40-60  (similar  to  Case  1)  and  the  same  is  true  with  the  greater  intensity  (4+)  icing 
PIREPs  (86%).  As  with  Case  1,  most  of  tiie  outlying  PIREPs  outside  of  tiie  range  can  be  explained  by 
questionable  grid  values  along  cloud  edges  due  to  the  relatively  large  gradient  of  the  values  across  a  short 
distance  in  tiie  vicinity  of  the  measurement  or  in  regions  of  broken  clouds  which  let  some  surface  radiation 
tiirough  between  the  clouds.  In  this  case,  the  4+  intensity  reports  also  generally  corresponded  to  greater 
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AT*  and  albedo  values.  However,  when  Aese  AT*  and  albedo  range  thresholds  are  again  ^lied  to  the 
spreadsheet  grid  of  pixel  dump  values  (shown  in  Tables  6.4a-d)  by  coloring  die  grid  sectors  which  are 
within  die  ranges,  the  resuh  is  die  same  as  in  Case  1  with  die  entire  cloud  shield  being  outlined.  These 
ranges  for  which  die  PIREPs  clustered  just  proved  that  diey  occurred  in  clouds,  which  is  a  given,  without 
respect  to  eidier  type  of  cloud. 

When  the  cumuliform  cloud  tops  (convective  CTT  zone  (CZ))  were  screened  out  leaving  only  the 
stratiform  cloud  tops  (see  Figures  6.9a-b)  (icing  CTT  zone  (IZ)  shaded  in  the  color  pink)  grid  sectors 
(shown  in  Tables  6.4c-d),  a  better  indication  of  icing  conditions  resulted  similarly  to  Case  1.  The  shaded 
regions  in  Tables  6.4c-d  for  the  CH2-CH4  (AT*)  grid  values  again  represented  stratiform  regions  better. 
Unlike  Case  1,  the  AT*  values  for  die  IZ  were  in  die  evenly  distributed  within  the  range  (10-44).  But, 
again  diese  range  of  values  were  not  characteristic  of  icing  regions  only,  and  diere  seemed  to  be  nearly  as 
many  instances  of  these  values  in  the  unshaded  CZ  portion  of  the  grid.  This  leads  to  a  similar  result  for 
Case  1,  that  die  channel  difference  method  is  not  a  particularly  stand-alone  technique  for  icing  detection. 

Highlighting  the  region  of  the  IZ  using  only  CH4  T*  on  the  grid  (see  Figure  6.10)  again  gives  a 
better  visual  representation  of  the  stratiform  areas  against  the  visible  CHI  image  for  that  tune. 
Comparatively,  this  case  was  much  less  convective  overall  dian  in  Case  1  as  can  be  seen  fiom  the  larger 

stratiform  cloud  regions  in  the  individual  GOES-8  channel  imagery  (see  Figures  6.1  la-c).  Transferring 

the  stratiform  cloud  area  to  the  composite  PIREP  plot  for  the  6  hr  period  encompassing  the  analysis  time 
^j5.2lZ)  (see  Figure  6.12)  gives  an  idea  of  how  die  PIREPs  compared  to  the  stratiform  regions  alone.  It  is 
obvious  from  inspection  of  the  icing  PIREP  and  cloud-top  CH4  temperature  (CTI)  data  that  also  during 
this  period  most  of  the  icing  PIREPs  coincided  with  very  cold  CTTs.  The  greater  intensity  icing  PIREPs 
(IP)  of  TnagnitiiHa  4+  were  again  mainly  in  a  convective  CTT  zone  between  248K  and  205K  (-25^  to  - 
68®C)  (see  Figure  6.9b)  over  the  mid  Mississippi  River  Valley  area.  The  PIREPs  within  die  stratiform 
region  were  more  abundant  in  Case  2,  but  there  was  no  clear  majority  of  PIREPs  of  the  significant  4+ 
intensity  (though  most  occurred  over  eastern  Nebraska,  Kansas  and  Arkansas).  In  fact,  the  amount  of  4+ 
PIREPs  overall  decreased  from  Case  1  (46  to  44),  while  die  overall  number  icing  PIREPs  increased  from 
86  to  109.  This  result  may  indicate  that  the  more  significant  icing  cases  over  the  entire  cloud  region  is  due 
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mainly  to  convective  cloud  encounters  (as  seen  in  Case  1),  however  more  but  less  intense  icing  encounters 
ocCTir  wMiin  stratiform  cloud  areas.  A  rawinsonde  analysis  from  Omaha,  NE  (OAX)  (Figure  6.14),  for 
Topeka,  KS  (TOP)  (Figure  6.15),  and  frmn  other  stations  near  tire  stratiform  cloud  regions  (Figures  6.13 
and  6.15-17)  shows  an  atmosphere  that  would  be  conducive  to  icing  conditions  at  altitudes  between  tire 
surfrwe  and  12000  ft  AGL.  Again,  some  clustering  of  moderate  intensity  (4+)  icing  PIREPs  were  reported 
in  these  areas  and  at  these  altitudes,  which  again  confirms  that  this  area  might  be  a  good  candidate  for  this 
satellite  analysis  technique. 

6.13  RGB  Method  for  10  Apr  95 

When  the  RGB-411  combination  technique  is  used  for  imagery  taken  at  10/1815Z  Apr  95  (see 
Figure  6.18),  several  cloud  features  can  be  readily  distinguished  by  visual  interrogation  when  checked 
against  individual  charmel  imagery  for  confirmation.  The  lowest  clouds  appear  a  dark  green-blue  or  cyan 
color.  Middle  clouds  transitioning  from  liquid  to  ice  phase  clouds  tend  to  correlate  with  a  light  shade  of 
cyan  and  the  lightest  pinkish  hues.  Then,  the  deepest  (convective  and  likely  ice-topped)  clouds  with  the 
coldest  tops  appear  white  or  nearly  white.  Some  upwelling  surface  radiance  seems  to  be  bleeding  through 
the  northern  fiinges  of  the  thin,  cirrus  clouds  fiom  Norfli  Dakota  tiirough  the  Great  Lakes,  making  them 
appear  red  or  pink  instead  of  a  whiter  shade.  Using  tiiis  combination  seems  particularly  usefiil  in 
distinguishing  low,  middle  and  high  clouds  with  the  convective  thunderstorm  tops  very  evident  as  the 
brightest  whites.  As  far  as  icing  detection,  the  majority  of  the  4+  icing  PIREPs  (over  western  Kansas, 
eastern  Nebraska  and  South  Dakota,  soutiiem  Minnesota  and  nortiiem  Illinois;  see  Figure  6.5)  seem  to 
correspond  to  the  darkest  pink  shades  and  lightest  blue  diades  with  no  consistent  hue  being  dominant 
throughout  the  region.  Also,  there  are  not  enough  PIREPs  in  other  regions  of  the  image  (western  South 
Dakota  and  Nebraska  and  with  tiiese  colors  to  make  a  distinct  or  consistent  probability  for  icing  conditions. 
For  use  as  a  forecasting  tool,  this  combination  may  be  helpful  to  outline  areas  as  a  first  guess  for  icing 
regions. 

A  qualitative  inspection  of  the  RGB-241  imagery  from  the  same  date  and  time  assists  in  further 
resolving  the  highest  clouds  (bright  white,  light  pinks,  and  yellows)  from  tiie  middle  and  lower  cloud  decks 
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(see  Figure  6.19).  Again  the  lowest  clouds  appear  dark  blue  with  middle  clouds  in  a  color  range  from  dark 
cyan  to  lightw  cyan  for  hi^er  middle  clouds,  and  die  yellow  hue  clouds  represent  diin  cirrus.  Then,  die 
pmlfish  or  lavender  hues  seem  to  signal  the  change  from  the  liquid  to  ice  phase  clouds,  since  their  presence 
over  Colorado,  Wyoming,  and  South  Dakota  correlate  weU  with  a  conspicuous  lack  of  positive  icing 
PIREPs.  Yet,  die  pinkish  tinged  cloud  tops  over  eastern  Nebraska  are  in  an  area  of  4+  intensity  icing 
PIREPs.  This  may  be  a  resuh  of  its  close  proximity  to  die  deeper  convection  occurring  just  to  the  south 

(bright  white  clouds). 

6.2  Conclusions 

Whereas  Ellrod  (1996a)  confined  himself  to  only  applying  diis  technique  to  nearly  uniform 
stratiform  (mainly  "warm")  cloud  regions,  diis  study  first  attacked  a  bigger  problem,  namely  die  entire 
synoptic  system  widi  both  convective  (cumuliform)  and  stratiform  regions  at  different  levels,  to  see  how  it 
would  do  at  detecting  icing  conditions  for  all  types  of  clouds,  which.  Evidence  from  diis  study  agrees  diat 
it  may  be  wise  to  concentrate  on  stratiform  regions  as  opposed  to  cumuliform  cloud  areas,  since  the 
satellite-based  sensor  has  a  hard  time  in  die  visible  and  IR  "seeing"  into  the  cloud  and  inferring  inner  cloud 
microphysics  or  phase  dirough  verticaUy  extensive  and  non-uniform  convective  clouds.  This  may  account 
for  the  general  abundance  of  icing  reports  in  the  convective  portions  of  the  cloud  shield  while  CTTs  were 
much  to  cold  to  indicate  a  support  for  SLW  collection.  Conclusive  results  may  have  been  better  had  die 
study  taken  place  in  February  or  March  when  a  climatological  maximum  of  icing  PIREPs  occurs  (Brown, 
1995b)  or  if  more  aircraft  had  chosen  to  fly  dirough  some  of  the  probable  areas  with  icing  conditions. 

In  trying  to  apply  the  cloud-top  sensed  ATb  technique  to  detection  of  aircraft  icing  within  all 
clouds  in  the  areas  of  interest,  die  results  from  die  two  cases  and  techniques  used  in  this  study  were  not 
definitive  enough  yet  to  be  a  helpful  technique  for  determining  icing  potential.  Convective  CTTs  tended 
to  produce  erroneous  assumptions  of  cloud  phase  and  SLW  content  when  only  sensing  their  cloud  tops. 
The  ATb  technique  seems  to  work  best  only  when  a  specific  area  of  stratiform  clouds  can  first  be  specified. 
When  it  comes  down  to  the  efficiency  of  this  technique  within  the  forecaster's  workload,  this  method 
unless  highly  automated,  would  not  seem  to  provide  any  added  value  to  the  icing  detection  problem.  It 
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may  be  that  the  grid  method  in  which  I  used  is  6r  loo  coaise  in  relating  specific  image  pixeb  to  the 
specific  location  of  each  PIREP,  and  better  results  may  come  from  a  highly  selective  examination  of  the 
exact  PIREP  and  pixel  co-location  and  corresponding  pixel  values.  This  ^)pears  to  be  the  way  Ellrod 
(1996a)  obtained  his  positive  results. 

This  ATb  method  seems  destined  for  use  on  "warm"  horizontally  consistent  (stratiform)  cloud 
tops,  instead  of  "cold"  (cumulifoim)  cloud  tops.  The  nature  of  deep  convective  clouds  is  such  that  their 
tops  are  often  glaciated  (consisting  mainly  of  ice  crystals),  however  the  stronger  updrafts  within  the  cells 
also  tend  to  accumulate  supercooled  liquid  water  (SLW)  in  regions  within  the  cloud  and  sometimes  near 
cloud  top.  Thus,  a  varying  array  of  satellite-sensed  IR  temperature  fliresholds  would  be  required  to  try  to 
pinpoint  icing  regions  within  deep  cumuliform  clouds,  but  probably  without  much  success. 

Satellite-based  microwave  retrieval  methods  hold  some  promise  for  Aese  types  of  clouds, 
however  their  returns  are  affected  by  precipitation-sized  particles,  which  are  also  quite  likely  within  diese 
cloud  types.  Promising  results  from  Lawyer  (1995)  (further  discussed  in  Chapter  2)  using  satellite 
microwave  cloud  liquid  water  (CLW)  retrieval  techniques  that  can  penetrate  thin  ice  cloud  tops  have 
shown  that  potential  icing  conditions  (verified  by  significant  intensity  icing  PIREPs)  can  be  inferred  using 
CLW  measurements  (CLW  >  0.5  kg  m*^)  in  conjunction  with  a  range  of  CTTs  (warmer  than  -25‘'C)  firom 
LIR  imagery.  His  CTT  ranges  seem  to  agree  wifli  die  icing  zone  CTTs  used  in  this  study,  however  this 
study  did  not  reach  any  quantitative  conclusions  on  liquid  water  content  of  the  clouds  in  the  areas  of 
interest.  Uwyer  found  that  very  little  CLW  existed  when  CTTs  were-colder  than  -25»C,  which  would 
imply  increasing  numbers  of  ice  crystals  scavenging  die  available  SLW.  Also,  die  CLW  amounts  were 
more  variable  at  CTTs  warmer  than  -25®C,  showing  the  unstable  nature  of  the  SLW  within  the  cloud 
(possibly  an  indication  of  the  clouds  life  cycle).  Lawyer's  work  would  also  further  indicate  that  microwave 
CLW  retrievals  in  combination  widi  other  sensing  techniques  (multispectral  satellite,  ground-based  remote 
sensing,  or  numerical  weather  prediction  model  output)  would  be  necessary  to  create  improved  icing 
forecasting  tools  for  the  future 

Most  aircraft  pilots  have  been  educated  that  the  existence  of  a  deep  cumulifoim  cloud 
(cumulonimbus  and  "towering"  cumulus  congestus)  implies  the  hazard  for  icing  conditions  among  other 
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tilings,  and  tiicy  tend  to  steer  clear  of  tiiese  areas.  Hus  demonstrates  anotiier  reason  that  tiiere  are  not 
sufficient  PIREPs  available  to  verify  the  existence  of  icing  regions  witiiin  tiiese  clouds  any  way.  Yet, 
aircraft  pUots  do  not  seem  averse  to  flying  tiirou^  tiiick,  non-precipitating  regions  of  stratiform  clouds 
even  if  some  ir-ing  has  been  forecasted,  as  was  seen  fiom  the  Cole  and  Sand  (1991)  study  of  aircraft  icing 
accidents.  In  bofli  cases,  tiie  majority  of  the  icing  PIREPs  occurred  within  areas  outside  of  (with  CTTs 
colder  than)  the  typical  icing  CTT  zone  (IZ)  (73%  within  the  CZ  and  27%  in  the  IZ  in  Case  1),  however 
many  cignifirant  intensity  icing  PIREPs  did  occur  witiiin  tiie  IZ  in  Case  2  (44%  within  the  IZ)  when  the 
cloud  shield  became  more  stratiform.  Ihe  icing  PIREP  verification  was  feirly  inconclusive  for  stratiform 
regions  due  to  a  relative  lack  of  icing  reports  in  some  portions  of  the  stratiform  cloud  deck.  It  is  important 
to  note  that  the  icing  PIREP  intensities  are  used  witiiout  regard  to  the  type  of  aircraft  which  might  have 
reported  it,  as  intensities  will  vary  on  aircraft  type  and  their  individual  anti-icing  c^abilities  (noted  in 
Ch^ter  3). 

However,  tiiere  would  be  an  underlying  problem  for  the  operational  aviation  weatiier  forecaster 
who  would  use  the  ATb  satellite  techniques  from  this  study.  The  feet  that  icing  only  in  stratiform  clouds 
would  be  the  main  ^plication  for  this  ATb  technique  leaves  a  whole  realm  of  convective  clouds  for  which 
there  is  yet  no  clear-cut  indicator  of  icing  conditions  in  cumuliform  clouds.  This  is  especially  important 
since  these  cloud  types  make  up  the  majority  of  tiie  clouds  near  surface  frontal  systems,  which  are  the  most 
common  locations  for  icing  conditions.  Even  satellite  microwave  cloud  liquid  water  retrieval  methods 
have  problems  with  the  majority  of  the  convective  (precipitating)  clouds.  So,  the  future  emphasis  of  icing 
conditions  detection  and  prediction  research  should  likely  focus  on  convective  clouds  with  the  most  likely 
techniques  coming  fixim  mesoscale  numerical  weather  prediction  models  or  precipitation  radars 
that  can  someday  closely  simulate  or  detect  individual  cloud  system  life  cycles.  Model  output  of  inner 
cloud  structure,  when  combined  with  infonnation  from  satellite  and/or  ground-based  remote  detection 
systems,  should  help  to  pinpoint  tiie  atmospheric  layers  for  which  icing  conditions  are  most  probable.  This 
would  fnavimiyf!  tiie  efficiency  of  the  icing  forecast  product,  while  enabling  safe  flight  through  more 
regions  of  the  atmosphere. 
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The  multispectial  RGB  color  enhancement  meAod  does  a  good  job  at  distinguishing  clouds  of 
different  heights  widi  the  greatest  color  and  shade  contrast  between  die  high  and  low  clouds.  Some 
infonnation  on  cloud  phase  and  SLW  clouds  seems  evident  from  inspection  of  the  middle  cloud  deck  color 
hues,  but  the  resuh  is  not  yet  reliably  quantitative.  This  infonnation  could  be  used  qualitatively  to  help  in 
initial  icing  assessment  for  larger-scale  cloud  regions.  For  die  prospects  of  the  multispectral  RGB 
enhancement  method  for  detecting  icing  conditions,  diese  color  enhancements  provide  an  easier 
identification  of  cloud  heights  and  can  be  used  in  distinguishing  stratiform  clouds  (relatively  darker  with  a 
more  uniform  shade)  from  cumulifoim  clouds  (mainly  white  or  very  light  color  shades).  While  there  was 
no  definitive  proof  found  that  these  methods  could  identify  potential  icing  regions,  these  do  seem  a 
valuable  tool  in  separating  clouds  into  their  respective  height  categories  or  types.  Though,  for  an 
operational  forecaster,  to  use  these  RGB  techniques  might  be  too  great  of  a  task  to  do  effectively  and 
consistently  with  out  much  practice  or  experience  at  it  Unfortunately,  while  the  human  eye  is  effective  at 
Hicringiiighing  these  subtle  hue  differences  in  diese  images,  for  proper  accuracy  and  time  efficiency,  it 
would  likely  need  to  be  automated  using  a  computer’s  help.  As  of  yet  it  has  not  been  shown  that  a 
computer  can  effectively  assimilate  color  changes  effectively  enough  to  develop  a  forecast  product  to  take 
advantage  of  diese  RGB  metiiods.  However,  neural  networic  (computer  learning  algorithms)  methods  are 
being  researched  presently  at  CSU/CIRA  and  odier  locations,  which  would  help  in  the  automated  computer 
determination  and  processing  of  pixel  color  hues  and  could  take  advantage  of  diese  methods. 

63  Suggestions  for  Future  Work 

Detection  of  potential  icing  environments  using  infrared  and  visible  satellite  imagery  remains  a 
challftnge  for  cloud  regions  made  up  of  convective  clouds  or  a  mix  of  both  convective  and  stratiform  types. 
Aircraft  iring  is  such  a  complex  phenomenon  based  on  many  highly  and  non-linearly,  temporally  and 
spatially  variant  atmospheric  and  mechanical  factors,  many  of  which  have  been  described  herein.  As  long 
as  the  aircraft  icing  problem  exists,  researchers  will  be  trying  new  scientific  techniques  to  try  to  unlock  die 
many  remaining  mysteries  of  its  formation,  detection,  and  prediction.  With  further  improvements  in  the 
next  generation  of  GOES  channel  imagery  and  better  analysis  and  display  system  capabilities,  the  wide 
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variety  of  ipiaga  channel  combinations  can  be  fuitiier  explored  to  he^)  in  tiie  detection  of  aircraft 

icing  environments  using  the  channel  differencing,  RGB  enhancement,  and  oflier  mefliods.  In  tiie  future, 
improved  icing  environment  detection  and  prediction  wffl  most  likely  result  fiom  a  combination  of  a 
variety  of  remote  sensing  techniques  both  active  and  passive,  such  as  ground>based  radiometers,  advanced 
precipitation  radars,  profilers,  in-shu  aircraft  measurements,  and  satellite-based  radiometers.  Also,  a  better 
understanding  of  the  microphysics  of  stratiform  and  cumuliform  clouds,  as  well  as  precipitation  formation 
processes,  must  be  reached  before  the  full  compliment  of  conditions  of  an  icing  region  can  be  successfully 
determined.  Many  more  icing  cases  with  better  detailed  PIREP  data  sets  or  other  verification  methods 
need  to  be  examined  before  a  definitive  detection  and  prediction  scheme  will  result. 

In  the  future,  as  numerical  models  become  more  sophisticated  in  terms  of  tiieir  microphysical  and 
oflier  parameterizations,  and  additional,  new  data  sources  (including  satellite  techniques)  are  brought 
together  to  form  improved  icing  algorithms,  icing  forecasts  should  continue  to  be  more  precise/accurate, 
while  reducing  the  overforecasting  of  temporal  and  spatial  coverage  (false  alarms).  Also,  future  studies  of 
regional  biases  in  the  performance  of  the  algorithms  and  those  directed  towards  developing  methods  to 
forecast  icing  severity/intensity  will  benefit  from  an  improved  suite  of  atmospheric  information  on  icing. 

Improvements  in  tiie  temporal  and  spatial  coverage  of  PIREPs  would  provide  a  much  improved 
basis  on  which  to  judge  the  effectiveness  of  future  aircraft  icing  forecasts.  At  the  moment,  a  more 
systematic  collection  of  icing  verification  would  likely  require  tiie  assistance  of  commercial  and  commuter 
airlines.  Until  FAA  rules  require  pilots  to  give  more  PIREPs  or  a  way  is  found  to  better  mechanically 
sense  the  ice  accretion  formation  and  automatically  record/transmit  more  frequent  in-flight  weather 
environment  information  (temperature,  relative  humidity,  cloud  droplet  sizes,  etc.),  icing  forecasting  and 
verification  will  be  limited  to  present  PIREP  verification  mefliods  and  smaller  scale  experiments  for 
evaluating  forecasts.  But,  of  course,  if  a  commercial  pilot  is  asked  to  fly  into  threatening  cloudy 
regions  or  suspected  aviation  hazard  areas  that  would  imply  possibly  severe  weather,  one  can  expect  that  it 
wouldn't  happen  voluntarily! 

Ultimately,  «gnificnnt  advances  in  aircraft  icing  forecasts  depend  on  better  observations  of  icing. 
Direct,  automated  measurements  of  airframe  ice  are  not  likely,  but  indirect  techniques  using  data  from  on- 
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board  sensors  of  temperature  and  humidity  might  be  feasible.  Several  major  airlines  are  currently  using 
equipment  for  collecting  temperature  and  wind  measurements.  This  data  is  relayed 
automatically  by  a  system  called  ACARS  (ARINC  Communications,  Addressing,  and  Reporting  System, 
operated  by  Aeronautical  Radio,  Inc.  (ARINC)).  They  cover  a  substantial  portion  of  the  continental  U.S. 
and  are  being  used  experimentally  for  model  initialization.  Although  they  do  not  yet  include  humidity 
measurements,  die  equipment  for  providing  these  measurements  is  presently  being  tested. 
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Table  6.1a  PlREPsfor  Case  1  (10  Apr  95  1500-2100Z) 
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Table  6.  lb  PIREPs  for  Case  1  (10  Ap-  95  1500-21002) 
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CH2-CH4  IRaTb  (deg  C) 
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CH2-CH4IRATs(degC)  CH2  -  CH4 IR  aTb  (deg  C) 


CH2  -  CH4 IR  Tb  Diff  vs  CHI  Visible  Albedo 


0  J 


20  30  40  50  60  70  80 

CHI  Visible  Albedo  (scaled  0  - 100) 


Figure  6.1c  Scatter  plot  of  null  and  all  icing  PIREPs  (Case  1) 


CH2  -  CH4  IR  Tb  Diff  vs  CHI  Visible  Albedo 


Figure  6.  Id  Scatter  plot  of  all,  moderate  or  greater  (4+),  and  extreme  (5+)  icing  PIREPs  (Case  1) 


6 


117 


118 


119 


Table  6.2d  Grid  of  pixel  dump  values  for  Case  1  from  1745Z  imagery 
I  Channel  Differenced  Pixel  Dump  T.  Value  (E-W  for  pbcst  number  out  of  1 
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il  number] 


IRCTT(degK)  CH4 IR  CTT  (deg  K) 


CH4  IR  Cloud-Top  Temp  vs  Icing  Intensity 


Figure  6.2b  Scatter  plot  of  all  PIREPs  within  the  Convective  Zone  (CZ)  (Case  1) 
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Figure  6.3  Grid  showing  shaded  sectors  representing  stratiform  cloud  tops  within  IZ  (Case  1 ) 


GOES-8  CHI  10/1745ZApr95 
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Figure  6.4a  GOES-8  channel  1  (visible)  image  of  U.S.  for  10/1745Z  Apr  95  (Case  1) 


GOES-8  CH2  10/1745Z  Apr  95 
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Figure  6.4c  GOES-8  channel  ^ 


GOES-8  CH4  10/1745Z  Apr  95 
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Figure  6.4d  GOES-8  channel  4  (inverted  LIR)  image  of  U.S.  for  10/1745Z  Apr  95  (Case  1) 


Figure  6.5  Composite  plot  of  all  PIREPs  from  1500-2100Z  on  10  Apr  95  (Case  1) 
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Figure  6.6  Skew-T,  log-P  diagram  for  Nor^  Platte,  NE  (LBF)  at  I0/1200ZApr  95  (Case  1) 
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Figure  6.7  Skew-T,  log-P  diagram  for  Dodge  City,  KS  (DDC)  at  10/l200ZApr  95  (Case  1) 
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Table  6.3a  PIREPsforCase2(lIApr95  1500-2100Z) 
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Table  6.3b  PIREPsforCase2(llApr95  1500-2100Z) 
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TtAle  6.3c  PlREPs  for  Case  2  (11  Apr  95  1500-2100Z) 
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CH2-CH4  IRATstdeg  C) 


CH2  -  CH4 IR  Tb  Diff  vs  CHI  Visible  Albedo 
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Figure  6. 8a  Scatter  plot  of  null  (0  intensity  -  no  icing)  PIREPs  (Case  2) 


Figure  6.8b  Scatter  plot  of  all  icing  (1-7  intensity)  PIREPs  (Case  2) 
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CH2  -  CH4  IR  Tb  Diff  vs  CHI  Visible  Albedo 
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Figure  6.8c  Scatter  plot  of  null  and  all  icing  PIREPs  (Case  2) 


CH2  -  CH4  IR  Tb  Diff  vs  CHI  Visible  Albedo 
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Figure  6.8d  Scatter  plot  of  all,  moderate  or  greater  (4+),  and  extreme  (5+)  icing  PIREPs  (Case  2) 
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TaUe  6.4b  Grid  of  pixel  dump  values  for  Case  2  from  174SZ  imagery 
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Table  6.4c  Grid  of  pixel  dump  values  for  Case  2  from  1745Z  imagery 
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Table  6.4d  Grid  of  pixel  dump  values  for  Case  2  from  1745Z  imagery 
BrUS  I  I  Chwind  Differenced  Wxel  Dump  T«  Value  (E-W  for  pixel  numt 
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Icing  Intensity  (0  •  7  /  low  ->  high) 


Figure  6.9a  Scatter  plot  of  all  PIREPs  (0-7  intensity)  within  the  Icing  Zone  (IZ)  (Case  2) 


CH4 IR  Cloud-Top  Temp  vs  Icing  Intensity 


Icing  Intensity  (0  -  7  /  low  ->  high) 


Figure  6,9b  Scatter  plot  of  all  PIREPs  within  the  Convective  Zone  (CZ)  (Case  2) 


aarw 

•!  . 

1 ;  '  f  If 

li/x 

t* 

t  J 

\ 

■Hianri 

■ 

El 

E 

'tA' 

'’  •  ..r  „ 

OBI 

BB 

B 

■IBBBI 

■ 

ai'*!' 

X 

JiM 

'^3 

8 

B 

BB 

B 

■IBBBI 

E 

7. 

r7*»>^ 

Ml 

riiBB 

B 

BBB 

B 

■IBBBI 

n 

E 

’^i ' 

BB 

B 

B 

BB 

o 

u 

Aui 

-•  *•> 

^iiaesB 

BB 

B 

B 

BB 

B 

■■■■1 

S 

/■'? 

.IZ311 

B 

BB 

BBSS 

BH 

B 

mm 

B 

BB 

mmi 

B 

BBB 

■■■■1 

Ilf 

^  r: 

IDBBB 

BB 

BiJ 

□ 

B 

BB 

n 

■■■■1 

■ 

p. 

T^-  V:  *-.  •.’? 

mm 

B 

O 

BB 

BBBBBI 

■ 

BB 

lZBRB 

Gi: 

BB 

B 

BBBBSBBBI 

BBHHBKaHHaBnaaaBaa 

B 

■smi 

B 

B 

BBSai 

BIIB 

B 

BB 

BB 

BB 

BB 

B 

EBHBI 

V  - 

i  1 

B 

[:.lB 

BB 

BIIB 

B 

BB 

HB 

B 

B 

BBBBBEllll 

B 

li 

i^lB 

BB 

BIIB 

BBB 

IIB 

B 

B 

BE 

;v. 

B 

1  0»l 

is  :l 

i  f.  A 

IjBBB 

BIIB 

B 

BBlfBB 

B 

□B 

B 

Billll?.l 

B 

ii^rjBBB 

BIIB 

B 

BB 

IIB 

B 

B 

BB 

B 

Bll^^l 

BB 

SB 

BB 

BIIB 

BBB 

«B 

B 

B 

BB 

B 

■Esl^XiI 

B 

B 

BlIBBSIIB 

B 

BBWBB 

n 

BB 

B 

ElilSISI 

BBBlIBBBKBBBfillEll' 

Bll2:iMllEli^O 

IB 

B 

Bll 

BB 

BBB 

B 

[2M^LS' 

HiiiiMMaHMaMaiqittAM. 'MMiirK  /?\ 

BBBWBBBBH 

B 

^Kr 

•JK  1* 

BBBHBBBLB 

E 

TTf-^ 

iM|n 

-ss 

B 

BBIBBflEl! 

,~:2Sr 

a  “ 

S 

1 

«BSBBBL:i~; 

mm 

(mICl 

31 

hi 

•- 

BB 

□BBEie'. 

'•rs 

iir*T 

1 

L  '  ** 

BB 

mz 

BBE 

T 

IiSi22 

Ei 

‘iy.l 

^  1 

■■“■.  1 

ii: 


■nilBHi  C^naHBHiWiSOIII^O 'm 

- If  ji/  v!r-Tr5SB*«BBiBiBif’~«rr'*»’rtf - 


«&££![  r 
m 
m 


iHB! 


ilia; 


■  (> 


MmnriA-m 
fe  ‘ ;  l«J*«i.-x-.?55S3ailSi33Ziilci'® ^  : 
:}»S  '-;5K:W30a3SS33SELllt;?  : 

. ....:_  .:::n3s:sE“-ini*iE;‘ 

irK^r'-'iES  iS^ESIiniBrSli  Jil. 


^  k  li.  .'J^q 

jr.-.  riC'J*a 
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Figure  6.1  la  GOES-8  channel  1  (visible)  image  of  U.S.  for  1 1/1745Z  Apr  95  (Case  2) 


GOES-8  CH2  1 1/1745Z  Apr  95 
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Figure  6. 


GOES-8  CH4  11/1745Z  Apr  95 
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Figure  6.1  Ic  GOES-8  channel  4  (inverted  LIR)  image  of  U.S.  for  1 1/1745Z  Apr  95  (Case  2) 
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TRC  Figure  6.12  Composite  plot  of  all  PlRElH^fMifi'f 500-2  lOOZ  on  1 1  2) 
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Figure  6.13  SkewT.  log-P  diagram  for  North  Platte.  NE  (LBF)  at  11/12002 Apr  95  (Case  2) 
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Figure  6.14  Skew-T.  log-P  dit^am  for  Omaha,  NE  (OAJQ  at  1  l/1200ZApr  95  (Case  2) 
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Figure  6.15a  Skew-T,  log-P  diap’om  for  Topeka,  KS  (TOP)  at  ll/1200ZApr  95  (Case  2) 
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Figure  6. 15b  Skew-T.  log-P  diagram  far  Denver,  CO  (DNR)  at  1 1/1200Z  Apr  95  (Case  2) 
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MULTI  -411 
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Figure  6.18  Multispectral  RGB-41 1  imagery  for  10/1815Z  Apr  95 


MULTI  -  241 
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Figure  6.19  Multispectral  RGB-241  imagery  for  10/1815Z  Apr  95 
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